资源描述:
《数字信号处理 matlab参考资料new》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、MATLAB相关内容一信号与系统一、目的1、熟悉连续信号理想采样前后的频谱变化,加深对时域采样定理的理解。2、熟悉时域二、内容1、理想采样信号序列(1)首先产生信号x(n),0<=n<=50n=0:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;%T为采样率w0=50*sqrt(2.0)*pi;%ω符号在MatLab中不能输入,用w代替x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π,信号乘可采用“.*”cl
2、oseall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title('理想采样信号序列');(2)绘制信号x(n)的幅度谱和相位谱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('理想采
3、样信号序列的相位谱');(3)改变参数为:A=1,a=0.4,w0=2.0734,T=1A=1;a=0.4;w0=2.0734;T=1;%设置信号有关的参数和采样率Tx=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;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);
4、%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('理想采样信号序列的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('理想采样信号序列的相位谱');2、单位脉冲序列在MatLab中,这一函数可以用zeros函数实现:n=1:50;%定义序列的长度是50x=zeros(1,50);%注意:MATLAB中数组下标从1开始x(1)=1;closeall;subplot(3,1,1);stem(x);title
5、('单位冲击信号序列');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('单位冲击信号的相位谱');3、矩形序列n=1:50;x=sign(sign(10-n)+1);closeall;subplot(3,1,1);stem(x);tit
6、le('矩形信号序列');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('矩形序列相位谱');4、特定冲击串n=1:50;%定义序列的长度是50x=zeros(1,50);%注意:MATLAB中数组下标从1开始x(1)=1;x(2)=2.5;x(
7、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('特定冲击串序列相位谱');MATLAB相关内容二卷积计算及定理
8、一、目的利用卷积方法观察分析信号、系统的频谱特性二、内容1、卷积计算在MATLAB中,提供了卷积函数conv,即y=conv(x,h),调用十分方便。n=1:50;%定义序列的长度是50hb=zeros(1,50);%注意:MATLAB中数组下标从1开始hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;closeall;s