Question

matlab数字信号处理的采样用程序怎么实现？

Answer 1

1、理想采样信号序列\x0d\x0a（1）首先产生信号x(n),0<=n<=50\x0d\x0an=0:50; %定义序列的长度是50\x0d\x0aA=444.128; %设置信号有关的参数\x0d\x0aa=50*sqrt(2.0)*pi;\x0d\x0aT=0.001; %采样率\x0d\x0aw0=50*sqrt(2.0)*pi; %ω符号在MatLab 中不能输入，用w 代替\x0d\x0ax=A*exp(-a*n*T).*sin(w0*n*T); %pi 是MATLAB 定义的π，信号乘可采用“.*”\x0d\x0aclose all %清除已经绘制的x(n)图形\x0d\x0asubplot(3,1,1);stem(x); %绘制x(n)的图形\x0d\x0atitle(‘理想采样信号序列’); %设置结果图形的标题\x0d\x0a（2）绘制信号x(n)的幅度谱和相位谱\x0d\x0ak=-25:25;\x0d\x0aW=(pi/12.5)*k;\x0d\x0aX=x*(exp(-j*pi/12.5)).^(n’*k);\x0d\x0amagX=abs(X); %绘制x(n)的幅度谱\x0d\x0asubplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);\x0d\x0aangX=angle(X); %绘制x(n)的相位谱\x0d\x0asubplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)\x0d\x0a（3）改变参数为： 1, 0.4, 2.0734, 1 0 A = α = Ω = T =\x0d\x0an=0:50; %定义序列的长度是50\x0d\x0aA=1; a=0.4; w0=2.0734; T=1; %设置信号有关的参数和采样率T\x0d\x0ax=A*exp(-a*n*T).*sin(w0*n*T); %pi 是MATLAB 定义的π，信号乘可采用“.*”\x0d\x0aclose all %清除已经绘制的x(n)图形\x0d\x0asubplot(3,1,1);stem(x); %绘制x(n)的图形\x0d\x0atitle(‘理想采样信号序列’);\x0d\x0ak=-25:25;\x0d\x0aW=(pi/12.5)*k;\x0d\x0aX=x*(exp(-j*pi/12.5)).^(n’*k);\x0d\x0amagX=abs(X); %绘制x(n)的幅度谱\x0d\x0asubplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);\x0d\x0aangX=angle(X); %绘制x(n)的相位谱\x0d\x0asubplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)\x0d\x0a2、单位脉冲序列\x0d\x0a在 MatLab 中，这一函数可以用zeros 函数实现：\x0d\x0an=1:50; %定义序列的长度是50\x0d\x0ax=zeros(1,50); %注意：MATLAB 中数组下标从1 开始\x0d\x0ax(1)=1;close all;\x0d\x0asubplot(3,1,1);stem(x);title(‘单位冲击信号序列’);\x0d\x0ak=-25:25;\x0d\x0aX=x*(exp(-j*pi/12.5)).^(n’*k);\x0d\x0amagX=abs(X); %绘制x(n)的幅度谱\x0d\x0asubplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);\x0d\x0aangX=angle(X); %绘制x(n)的相位谱\x0d\x0asubplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)\x0d\x0a3、矩形序列\x0d\x0an=1:5; 0x=sign(sign(10-n)+1);\x0d\x0aclose all; subplot(3,1,1); stem(x);title(‘单位冲击信号序列’);\x0d\x0ak=-25:25; X=x*(exp(-j*pi/25)).^(n’*k);\x0d\x0amagX=abs(X); %绘制x(n)的幅度谱\x0d\x0asubplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);\x0d\x0aangX=angle(X); %绘制x(n)的相位谱\x0d\x0asubplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)\x0d\x0a4、特定冲击串\x0d\x0ax(n) =δ (n) + 2.5δ (n ?1) + 2.5δ (n ? 2) +δ (n ? 3)\x0d\x0an=1:50; %定义序列的长度是50\x0d\x0ax=zeros(1,50); %注意：MATLAB 中数组下标从1 开始\x0d\x0ax(1)=1; x(2)=2.5; x(3)=2.5; x(4)=1;\x0d\x0aclose all; subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);