3.1matlab应用基础.ppt

习题4参考答案 （2）t*u(t)的导数。 f=sym((t)*(heaviside(t))) f = t*heaviside(t) df=diff(f) df = heaviside(t) + t*dirac(t) 习题4参考答案 （3）(t-1)*u(t-1)的导数 f=sym((t-1)*(heaviside(t-1))) f = heaviside(t - 1)*(t - 1) df=diff(f) df = heaviside(t - 1) + dirac(t - 1)*(t - 1) （4）exp(2*j*t)的微分和积分 f=sym(exp(2*j*t)) f = exp(2*j*t) df=diff(f) df = 2*j*exp(2*j*t) intf=int(f) intf = exp(2*j*t)/(2*j) 习题5答案:δ(n) 习题5答案:u(n) * MATLAB 矩阵实验室 * MATLAB 矩阵实验室 * 注意：.* ./ 表示两个同维数组中的对应元素相乘除，或 其中之一为标量。 A.^B 表示对矩阵A的各元素与矩阵B的对应元素的乘方运算。A 与B必须为同维矩阵。 * 注意：plot命令不能在plot后给出座标轴的取值范围（与ezplot)对比，是因为取值范围已经由自变量在前面定义了取值范围。 * 注意：plot命令不能在plot后给出座标轴的取值范围（与ezplot)对比，是因为取值范围已经由自变量在前面定义了取值范围。 * 注意：z是指函数值，而非自变量的取值。 * 注：plot(f1,f2)表示f1为横轴，f2为纵轴 * 注意： plot(t,rectpuls(t-2,3))中的rectpuls(t-2,3）其中“ t-2” 表示中心点为t＝2，“3”表示宽度为3的矩形脉冲 * axis([a,b,c,d]):a与b表示x轴的显示范围，c与d表示y轴的显示范围。 * sinc(x)： sinc(x）＝ sinc(pi*x）/pi*x grid on:表示显示栅格 * (4)单位阶跃信号 2)调用单位阶跃函数Heaviside 调用单位阶跃函数Heaviside可方便的表示出单位阶跃信号，但是不能直接用ezplot来画图，因为函数ezplot只能画出即存在与Symbolic Math toolbox中，又存在于总MATLAB工具箱中的函数，而Heaviside只存在于Symbolic Math toolbox中。 (4)单位阶跃信号 解决方法： 在自己的工作目录work下创建Heaviside.m文件，该文件如下: Function f=Heaviside(t) f=(t0); 正确定义并保存该函数后，就可调用该函数了。 注意：调用该函数只能使用EZPLOT来画图，Heaviside是符号函数。 例：用matlab画出u(t+3)-u(t)的波形 f=sym(heaviside(t+3)-heaviside(t)) f = heaviside(t+3)-heaviside(t) ezplot(f,[-5,4]) 注意： 用ezplot时，坐标范围只能限定自变量的范围，而不能限定函数值的范围，如ezplot(f,[-5,4])后的范围里只有两项，表示所画出图形的自变量的范围。 (5) 取样函数： 用sinc(t)命令 N=1000; t=-10:20/N:10; y=sinc(t/pi); plot(t,y);grid on 用符号运算功能 y=sin(t)/t; ezplot(y) (6)指数信号 画出指数信号e-0.3t的波形； t=0:pi/50:4*pi; y=exp(-0.3*t); plot(t,y,r); title(exp(-0.3*t)); xlabel(t); ylabel(y); 指数信号e-0.3t的波形 画出指数信号2t的波形 t=0:0.01:3; y=2.^t; plot(t,y,r); title(2^t); xlabel(t); ylabel(y); 指数信号2t的波形 (7)随机信号 rand产生均匀分布的白嘈声，randn产生高斯分布的白嘈声 t=0:0.01:1; y=randn(1,length(t)); plot(t,y);grid on; 2.离散信号 （1）离散时间正弦信号 产生一个幅度为2,频率为1/12Hz,相位为 pi/6的正弦信号 w=2*pi/12; n=-10:10; x=2*sin(w*n+pi/6); stem(n,x,fill); grid on; (2) 离