matlab信号与系统代码整理.doc

连续时间系统 1 离散时间系统 2 拉普拉斯变换 4 Z变换 5 傅里叶 7 连续时间系统 %%%%%%%%%%向量法%%%%%%%%%%%%%%%% t1=-2:0.01:5; f1=4*sin(2*pi*t1-pi/4); figure(1) subplot(2,2,1),plot(t1,f1),grid on %%%%%%%%%符号运算法%%%%%%%%%%%% syms t f1=sym(4*sin(2*pi*t-pi/4)); figure(2) subplot(2,2,1),ezplot(f1,[-2 5])跟plot相比，ezplot不用指定t，自动生成。 axis([-5,5,-0.1,1])控制坐标轴的范围xx,yy； 求一个函数的各种响应 Y’’(t)+4y’(t)+2y(t)=f”(t)+3f(t) %P187 第一题 %(2) clear all; a1=[1 4 2]; b1=[1 0 3]; [A1,B1,C1,D1]=tf2ss(b1,a1); t1=0:0.01:10; x1=exp(-t1).*Heaviside(t1); rc1=[2 1];（起始条件） figure(1) subplot(3,1,1),initial(A1,B1,C1,D1,rc1,t1);title(零输入响应) subplot(3,1,2),lsim(A1,B1,C1,D1,x1,t1);title(零状态响应) subplot(3,1,3),lsim(A1,B1,C1,D1,x1,t1,rc1);title(全响应) Y=lsim(A1,B1,C1,D1,x1,t1,rc1);title(全响应)则是输出数值解 subplot(2,1,1),impulse(b1,a1,t1:t:t2可加),grid on,title(冲激响应) subplot(2,1,2),step(b1,a1,t1:t:t2可加),grid on,title(阶跃响应) 卷积 %第九题 P189 clear all; %(1) t1=-1:0.01:3; f1=Heaviside(t1)-Heaviside(t1-2); %定义信号 t2=0:0.01:4; f2=(Heaviside(t2-1)-Heaviside(t2-3)); %定义信号 figure(1) gggfconv(f1,f2,t1,t2); %计算卷积积分并绘出时域波形 离散时间系统 画图 x=[-1,2,3,3,5,-4]; n=[-2,-1,0,1,2,3]; figure(1) stem(n,x)离散序列画图函数 各种序列操作函数 %ex_2 clear all; x=[1 2 3 4 5 6 7 6 5 4 3 2 1]; n=[-4 -3 -2 -1 0 1 2 3 4 5 6 7 8]; % [x1_1,n1_1]=sigshift(2*x,n,5);%2x(n-5) [x1_2,n1_2]=sigshift(-3*x,n,-4);%-3x(n+4) [x1,n1]=sigadd(x1_1,n1_1,x1_2,n1_2);%2x(n-5)-3x(n+4) subplot(3,1,1),stem(n1,x1) % [x2_1,n2_1]=sigshift(x,n,-3);%x(n+3) [x2_2,n2_2]=sigfold(x2_1,n2_1);%x(-n+3) [x2_3,n2_3]=sigshift(x,n,2);%x(n-2) [x2_4,n2_4]=sigmult(x,n,x2_3,n2_3);%x(n)x(n-2) [x2,n2]=sigadd(x2_4,n2_4,x2_2,n2_2);%x(n)x(n-2)+x(-n+3) subplot(3,1,2),stem(n2,x2) % x3_2=zeros(1,length(x)); n3_2=n; for j=1:5 [x3_1,n3_1]=sigshift(x,n,j) [x3_2,n3_2]=sigadd(x3_2,n3_2,x3_1,n3_1) end x3=x3_2.*n3_2 n3=n3_2 subplot(3,1,3),stem(n3,x3) 复指数序列各因数画图 %P157 第8题 %(1) clear all; n1=-10:10; f1=exp(j*n1*pi) f1_m=abs(f1); f1_a=angle(f1); f1_r=real(f1); f1_i=imag(f1); figure(1) subplot(2,2,1),ste