matlab练习有程序超全

0102030405060-2000200理想采样信号序列010203040506005001000理想采样信号序列的幅度谱0102030405060-505理想采样信号序列的相位谱0102030405060-101理想采样信号序列0102030405060012理想采样信号序列的幅度谱0102030405060-505理想采样信号序列的相位谱0510152025303540455000.51单位冲击信号序列0102030405060012单位冲击信号的幅度谱0102030405060-505单位冲击信号的相位谱0510152025303540455000.51单位冲击信号序列01020304050600510单位冲击信号的幅度谱0102030405060-505单位冲击信号的相位谱特定冲击串卷积05101520253035404550024单位冲击信号序列01020304050600510单位冲击信号的幅度谱0102030405060-505单位冲击信号的相位谱05101520253035404550024系统hb[n]05101520253035404550-2000200输入信号x[n]0102030405060708090100-100001000输出信号y[n]020406005001000输入信号的幅度谱0204060-505输入信号的相位谱02040600510系统响应的幅度谱0204060-505系统响应的相位谱0501000500010000输出信号的幅度谱050100-505输出信号的相位谱024681012141601230246810121416024024681012141602460246810121416-0.500.5102468101214160123405101501234024681012141605101520程序一fora=1:3kk(a)=4-a;endfora=4:7kk(a)=a-3;endfora=8:15kk(a)=0;end05101501234024681012141605101520closeall;subplot(2,1,1);stem(kk);subplot(2,1,2);stem(abs(fft(kk,16)))程序二n=0:15;%定义序列的长度是15a=0.1;f=0.0625;x=exp(-a*n).*sin(2*pi*f*n);closeall;subplot(2,1,1);stem(x);subplot(2,1,2);stem(abs(fft(x)))程序三n=0:15;%定义序列的长度是15p=8;q=2;x=exp(-1*(n-p).^2/q);closeall;subplot(3,1,1);stem(abs(fft(x)))p=8;q=4;x=exp(-1*(n-p).^2/q);subplot(3,1,2);stem(abs(fft(x)))p=8;q=8;x=exp(-1*(n-p).^2/q);subplot(3,1,3);stem(abs(fft(x)))程序四n=1:50;%定义序列的长度是50hb=zeros(1,50);%注意：MATLAB中数组下标从1开始hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;closeall;subplot(3,1,1);stem(hb);title('系统hb[n]');m=1:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;%采样率w0=50*sqrt(2.0)*pi;x=A*exp(-a*m*T).*sin(w0*m*T);%pi是MATLAB定义的π，信号乘可采用“.*”subplot(3,1,2);stem(x);title('输入信号x[n]');y=conv(x,hb);subplot(3,1,3);stem(y);title('输出信号y[n]');k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,2,1);stem(magX);title('输入信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,2,2);stem(angX);title('输入信号的相位谱')Hb=hb*(exp(-j*pi/12.5)).^(n'*k);magHb=abs(Hb);%绘制hb(n)的幅度谱subplot(3,2,3);stem(magHb);title('系统响应的幅度谱');angHb=angle(Hb);%绘制hb(n)的相位谱subplot(3,2,4);stem(angHb);title('系统响应的相位谱')n=1:99;k=1:99;Y=y*(exp(-j*pi/12.5)).^(n'*k);magY=abs(Y);%绘制y(n)的幅度谱subplot(3,2,5);stem(magY);title('输出信号的幅度谱');angY=angle(Y);%绘制y(n)的相位谱subplot(3,2,6);stem(angY);title('输出信号的相位谱')程序五n=1:50;%定义序列的长度是50x=zeros(1,50);%注意：MATLAB中数组下标从1开始x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;closeall;subplot(3,1,1);stem(x);title('单位冲击信号序列');k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱')程序六n=1:50x=sign(sign(10-n)+1);closeall;subplot(3,1,1);stem(x);title('单位冲击信号序列');k=-25:25;X=x*(exp(-j*pi/25)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱')程序7n=1:50;%定义序列的长度是50x=zeros(1,50);%注意：MATLAB中数组下标从1开始x(1)=1;closeall;subplot(3,1,1);stem(x);title('单位冲击信号序列');X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱')程序8n=0:50;%定义序列的长度是50A=1;%设置信号有关的参数a=0.4;T=1;%采样率w0=2.0734;x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title('理想采样信号序列');k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('理想采样信号序列的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('理想采样信号序列的相位谱')程序9n=0:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;%采样率w0=50*sqrt(2.0)*pi;x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title('理想采样信号序列');k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('理想采样信号序列的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('理想采样信号序列的相位谱')