信号与系统实验一、二VIP免费

chapter1实验内容：1、画出以下连续时间信号的波形1-0)f(t)=cos(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=cos(2*pi*t);plot(t,fa);1-1)f(t)=sin(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=sin(2*t*pi);plot(t,fa);2-0)f(t)=Sa(t/π)代码如下：pi=3.14159;t=0:0.01:8;fa=sinc(t/pi);plot(t,fa);3-0)f(t)=2[u(t3)-u(t5)]代码如下：t=-1:0.01:10;ft=2*((t>=3)-(t>=5));plot(t,ft);axis([-1,10,0,3]);4-1)f(t)=et代码如下：t=0:0.01:10;ft=exp(t);plot(t,ft);4-2)f(t)=e-tu(t)代码如下：t=0:0.01:10;f1=(t>=0);f2=exp(-t);plot(t,f1.*f2);5-0）f(t)=2ej(π/4)t，画出实部、虚部、模和相角的波形代码如下：t=0:0.01:10;ft=2*exp(j*(pi/4)*t);h=real(ft);%实部g=imag(ft);%虚部r=abs(ft);%模a=angle(ft);%相角subplot(2,2,1),plot(t,h),title('实部')subplot(2,2,3),plot(t,g),title('虚部')subplot(2,2,2),plot(t,r),title('模')subplot(2,2,4),plot(t,a),title('相角')7）f(t)=u(t)代码如下：t=-1:0.01:5ft=(t>=0);plot(t,ft);axis([-1,5,0,1.5]);8）f(t)=δ(t)代码如下：t=-1:0.01:5;ft=(t>=0)-(t>=0.1);plot(t,ft);axis([-1,1,0,1.1]);9)f9为周期矩形信号，其幅度从-1到1，占空比为75%代码如下：pi=3.14159;t=-10:0.01/pi:10;ft=square(t,75);plot(t,ft);2、信号本身运算画出f1(t)为宽度是4，高为1，斜度为0.5的三角脉冲，然后画出f1(-t)，f1(2t)，f1(2-2t)的波形以及f1(t)的微分和积分波形。代码如下：t=-5:0.001/pi:5;ft=sawtooth(t*pi/2,0.5);plot(t,ft),title('三角脉冲')3、多个信号运算(1)信号f1(t)和f2(t)波形如图所示，画出f1(t)+f2(t)和f1(t)-f2(t)波形。代码如下：求f1(t)+f2(t)的波形：t=-4:0.01/pi:4;f1=(((t>=-2)-(t>=0)).*(t+1))+(((t>=0)-(t>=2)).*(t-1));f2=((((t>=-2)-(t>=0)).*((-t-2)*0.5)))+((((t>=0)-(t>=2)).*((-t+2)*0.5)));plot(t,f1+f2);求f1(t)-f2(t)的波形如下：t=-4:0.01/pi:4;f1=(((t>=-2)-(t>=0)).*(t+1))+(((t>=0)-(t>=2)).*(t-1));f2=((((t>=-2)-(t>=0)).*((-t-2)*0.5)))+((((t>=0)-(t>=2)).*((-t+2)*0.5)));plot(t,f1-f2);实验心得：MATLAB提供了大量的生成基本信号的函数，例如指数信号、正余弦信号。表示连续时间信号有两种方法，一是数值法，二是符号法。数值法是定义某一时间范围和取样时间间隔，然后调用该函数计算这些点的函数值，得到两组数值矢量，可用绘图语句画出其波形；符号法是利用MATLAB的符号运算功能，需定义符号变量和符号函数，运算结果是符号表达的解析式，也可用绘图语句画出其波形图。信号基本运算是乘法、加法、尺度、反转、平移、微分、积分，实现方法有数值法和符号法。实验方法归纳与总结：Matlab是一个很好的信号表示工具，利用matlab中的三角函数、sinc(t)函数、阶跃函数、三角函数sawtooth()，real()、imag()、abs()、angle()，以及各种运算符号，可以实现上述的波形，清晰明了。Chapter2实验内容：1、求解下列微分方程的零状态响应，并画出系统的响应波形。d2y(t)/dt2+2dy(t)/dt+100y(t)=10f(t),f(t)=sin(2πt)u(t)代码如下：pi=3.14159;t=0:0.01:5;sys=tf([10],[1,2,100]);f=sin(2*pi*t).*(t>=0);y=lsim(sys,f,t);plot(t,y),gridon;2、求下列系统的单位冲激响应和单位阶跃响应，并画出响应波形。d2y(t)/dt2+3dy(t)/dt+2y(t)=8f(t)单位冲击响应代码如下：t=0:0.001:5;sys=tf([8],[1,3,2]);f=(t==0);y=impulse(sys,t);plot(t,y);gridon;title('单位冲击响应');单位阶跃响应代码如下：t=0:0.001:5;sys=tf([8],[1,3,2]);f=(t>=0);y=impulse(sys,t);y=step(sys,t);plot(t,y);gridon;title('单位阶跃响应');3、利用符号计算法求下列微分方程的显式解，并画出响应的波形。d2y(t)/dt2+3dy(t)/dt+2y(t)=df(t)/dt+2f(t)f(t)=t2,y(0)=y'(0)=1求系统的自由响应yh(t)，强迫响应yp(t)，零输入响应yzi(t)，零状态响应yzs(t)以及系统的完全响应y(t)。完全响应y(t),代码如下：y=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=1','Dy(0)=1');ezplot(simplify(y));自由响应yh(t)，代码如下：y=dsolve('D2y+3*Dy+2*y=0','y(0)=1','Dy(0)=1');e...