matlab计算e有效值近似列表,MATLAB_实验5

实验5 MATLAB在信号与系统的应用

【实验目的】

1、了解并熟悉MATLAB实现常用连续时间信号的产生方法；连续系统的时域和频域分析方法。

2、掌握MATLAB程序的编程方法。

3、了解MATLAB信号处理工具箱函数的调用方法。

【实验内容】

1、对书上例题的总结

例1 绘制了单位冲激函数、单位阶跃函数与复指数函数的图形，用到的关键语句

单位冲激函数：x1=zeros(1,length(t)); %对所有信号初始化，t为时间

x1(n1)=1/dt; %给出t1时刻脉冲信号，dt为时间间隔

单位阶跃函数：x2 = [zeros(1,n1),ones(1,length(t)-n1)];%产生阶跃信

号,n1为t1对应的样本序号

例2 求解了LTI系统的零输入响应，用到了经典法求解，关键语句如下：

p=roots(a); %求系统的极点

V=rot90(vander(p));c= V\Y0'; % 用范德蒙特矩阵求对应于各极点的系数

for k=1:n y= y+c(k)*exp(p(k)*t);end

例3 求解n阶LTI系统的冲激响应，用到了residue函数求解极点和留数，还用到了for语句叠加各根分量。

例4 主要是卷积的计算。用到的函数有conv函数计算卷积，input函数等待用户输入数值，fliplr(h)函数将h左右翻转。

例5 求LTI系统的零状态响应，和求解零输入同样用到留数residue函数来求解冲激响应函数，并且绘制冲激响应函数的图形，然后将其与输入信号的卷积求出，即输出函数，最后绘制输出函数的图形。

例7 将周期方波分解为多次正弦波之和，用程序演示谐波合成情况。用到的新函数是line([0,pi+0.5],[pi/4,pi/4]) % 加上方波幅度线及标注

halft=ceil(length(t)/2);

mesh(t(1:halft),[1:10],y(:,1:halft)) %绘制三维图形