Matlab_DSP教程.ppt

function y=xy(x) if x=0 y=0 elseif x=1 y=1 elseif x=2 y=2*x else y=2*x+5 end （2）输入xy并保存这个M文件。 （3）在Mat lab命令窗口输入：xy(-1)，就有其值；输入：xy(0.7)，就可求出对应的值；依次可求得所有值。 2、循环语句 （1）for 循环变量=初值：增量：终值 循环体 end (2) while 条件 循环体 end (1)打开M文件，输入： s=0 for n=1:100 s=s+1/n^2 end （2）输入sg并保存这个M文件。 （3）在Mat lab命令窗口输入：sg ， 就有其值。 八、曲线的绘图 绘图命令fplot ----- 是专门用于绘制一元函数曲线。格式： fplot (‘fun’ ，[a , b]） 表示绘制区间[a ,b]上函数 y=fun 的图形。 Example5：绘出下列函数图形，说明其奇偶性。 fplot(x^4/2+x^2-1,[-10,10]) （b）按回车键Enter。 （c）出现图形窗口，上有其图形。 （d）同样的方法输入： （a）输入： fplot(x+sin(x),[-5,5]) （e）同样的方法，按回车键，出现其对应图形。 fplot(x^2*exp(-x^2),[-10,10]) 第四节 数学实验 1、MATLAB求导 2、MATLAB求驻点 3、矩阵的特征值与特征向量 一、MATLAB求导 1、导数运算命令 diff (f,x) diff (f,x,n) 2、运行格式 syms x 命令 syms x y = ….. d=diff (y,x,n) % 求 y的n阶导数 pretty(d) syms x y = 3*x-exp(x)/2+1 dy=diff (y,x,1) 或dy=diff(y) % 求 y的1阶导数 pretty(dy) syms x y1 y2 y1=2^(sin(x));y2=log(log(log(x))); dy1=diff(y1),dy2=diff(y2) dy1 = 2^sin(x)*cos(x)*log(2) dy2 = 1/x/log(x)/log(log(x)) pretty(dy1),pretty(dy2) Example2 求下列函数导数 syms x diff((cos(x))^2*log(x),2) ans = 2*sin(x)^2*log(x)-4*cos(x)/x*sin(x)-2*cos(x)^2*log(x)-cos(x)^2/x^2 pretty(ans) 二、MATLAB求驻点 求 y=f(x) 驻点命令格式为： x= solve(‘f(x)’) clear x=solve(x^2-x-6) x = [ -2] [ 3] 3、矩阵的特征值与特征向量 MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有3种： (1)E=eig(A) 求矩阵A的全部特征值，构成向量E。 (2)[V,D]=eig(A) 求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。 (3)[V,D]=eig(A,nobalance) 与第2种格式类似，但第2种格式中先对A作相似变换后求矩阵A的特征值和特征向量，而格式3直接求矩阵A的特征值和特征向量。 例 用3种不同的格式求A的特征值和特征向量。 命令如下： A=[1,2,2;1,-1,1;4,-12,1]; E=eig(A) [V,D]=eig(A) [V,D]=eig(A,nobalance) 第五节 信号处理应用 1、信号波形的MATLAB产生 2、用Simulink建立仿真系统 3、应用实例 三种方法 1.使用funtool 图形化函数计算器 可以自定义自变量x变化范围,缺省为 [-2pi,2pi], a为自由参数. 1.5sin(x)、5exp(-x)、sin(x)/x、sqrt(a*x) 2.用信号处理工具箱(Signal Processing toolbox) 3.向量表示法 t= -10:0.5:10 ;f=sin(t).