发布时间 : 星期五 文章长江大学信号与系统matlab实验答案更新完毕开始阅读
实验1 信号变换与系统非时变性质的波形绘制
? 用MATLAB画出习题1-8的波形。 ? 用MATLAB画出习题1-10的波形。 Eg 1.8
代码如下：
function [y]=zdyt(t) %定义函数zdyt
y=-2/3*(t-3).*(heaviside(-t+3)-heaviside(-t)); end
t0=-10;t1=4;dt=0.02; t=t0:dt:t1; f=zdyt(t); y=zdyt(t+3); x=zdyt(2*t-2); g=zdyt(2-2*t); h=zdyt(-0.5*t-1);
fe=0.5*(zdyt(t)+zdyt(-t)); fo=0.5*(zdyt(t)-zdyt(-t)); subplot(7,1,1),plot(t,f); title('信号波形的变化') ylabel('f(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,2),plot(t,y); ylabel('y(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,3),plot(t,x); ylabel('x(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,4),plot(t,g); ylabel('g(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,5),plot(t,h); ylabel('h(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,6),plot(t,fe); ylabel('fe(t)')
grid;line([t0 t1],[0 0]); subplot(7,1,7),plot(t,fo); ylabel('fo(t)')
grid;line([t0 t1],[0 0]); xlabel('Time(sec)')
结果：
信号波形的变化210-10210-10210-10210-10210-1010.50-1010-1-10f(t)-8-6-4-2024y(t)-8-6-4-2024x(t)-8-6-4-2024g(t)-8-6-4-2024h(t)-8-6-4-2024fe(t)-8-6-4-2024fo(t)-8-6-4-2Time(sec)024
Eg1.10
代码如下：
function [u]=f(t) %定义函数f(t) u= heaviside(t)-heaviside(t-2); end
function [u] =y(t) %定义函数y(t)
u=2*(t.*heaviside(t)-2*(t-1).*heaviside(t-1)+(t-2).*heaviside(t-2)); end
t0=-2;t1=5;dt=0.01; t=t0:dt:t1; f1=f(t); y1=y(t); f2=f(t)-f(t-2); y2=y(t)-y(t-2); f3=f(t)-f(t+1); y3=y(t)-y(t+1);
subplot(3,2,1),plot(t,f1); title('激励——响应波形图') ylabel('f1(t)')
grid;line([t0 t1],[0 0]);
subplot(3,2,2),plot(t,y1); ylabel('y1(t)')
grid;line([t0 t1],[0 0]); subplot(3,2,3),plot(t,f2); ylabel('f2(t)')
grid;line([t0 t1],[0 0]); subplot(3,2,4),plot(t,y2); ylabel('y2(t)')
grid;line([t0 t1],[0 0]); subplot(3,2,5),plot(t,f3); ylabel('f3(t)')
grid;line([t0 t1],[0 0]); subplot(3,2,6),plot(t,y3); ylabel('y3(t)')
grid;line([t0 t1],[0 0]); xlabel('Time(sec)')
结果：
激励——响应波形图10.50-210-1-210-1-2y3(t)f3(t)y2(t)f2(t)y1(t)f1(t)210-220-2-220-2-20246024602460246024602Time(sec)46
实验2 微分方程的符号计算和波形绘制
上机内容
? 用MATLAB计算习题2-1，并画出系统响应的波形。
注：由于MATLAB版本的不同，求齐次通解和非齐次通解中的常数C可能不同，
这时可令C1=1；C2=1；C3=1；C4=1；yht=subs(yht); yt=subs(yt); 则特解可求出yp=yt-yht
? 用MATLAB计算习题2-3，并画出系统响应的波形。 Eg2.1
代码如下：
C5=1;C6=1;C8=1;C9=1;
y=dsolve('D2y+3*Dy+2*y=4*exp(-3*t)','y(0)=3,Dy(0)=4') yht=dsolve('D2y+3*Dy+2*y=0')
yt=dsolve('D2y+3*Dy+2*y=4*exp(-3*t)') yp=yt-yht yh=y-yp
yzi=dsolve('D2y+3*Dy+2*y=0','y(0)=3,Dy(0)=4')
yzs=dsolve('D2y+3*Dy+2*y=4*exp(-3*t)','y(0)=0,Dy(0)=0') t=linspace(0,4,300);figure(1)
y_n=subs(y);yh_n=subs(yh);yp_n=subs(yp);
plot(t,y_n,t,yh_n,'m:',t,yp_n,'r-.','linewidth',2) xlabel('Time(sec)'),title('全响应，自由响应，强迫响应') legend('全响应','自由响应','强迫响应',0) figure(2)
yzi_n=subs(yzi);yzs_n=subs(yzs);
plot(t,y_n,t,yzi_n,'m:',t,yzs_n,'r-.','linewidth',2) xlabel('Time(sec)'),title('全响应，零输入响应，零状态响应') legend('全响应','零输入响应','零状态响应',0)
结果： >> eg2_1 y =
12*exp(-t) - 11*exp(-2*t) + 2*exp(-3*t)
yht =
C5*exp(-t) + C6*exp(-2*t) yt =
2*exp(-3*t) + C8*exp(-t) + C9*exp(-2*t) yp =
2*exp(-3*t) - C5*exp(-t) - C6*exp(-2*t) + C8*exp(-t) + C9*exp(-2*t) yh =
12*exp(-t) - 11*exp(-2*t) + C5*exp(-t) + C6*exp(-2*t) - C8*exp(-t) - C9*exp(-2*t)
yzi =
10*exp(-t) - 7*exp(-2*t)