1、打开搜旅MATLAB,新建脚本。
2、在编辑器里输入代码。
3、Tp是观察时间,我们设慎漏穗置为64微妙,即64/1000;接着输入的代码是产生M长采样序列x(n)。
4、把实验内容中给定的值设置好;再把给定的模拟信号式子输入进去。
5、subplot(3,2,1)的意思是:三行两列,在第一个位置作图。
6、代码输入完毕后,点击运行。
7、得宽卜到结果,对信号进行采样。
1、理想采样信号序列(1)首先产生信号x(n),0<=n<=50
n=0:50%定义序列的长度是50
A=444.128%设置信号有关的参数
a=50*sqrt(2.0)*pi
T=0.001%采样率
w0=50*sqrt(2.0)*pi%ω符号在MatLab 中不能输入,用w 代替
x=A*exp(-a*n*T).*sin(w0*n*T)%pi 是MATLAB 定义的π,信号乘可采用“.*”
close all %清除已经绘制的x(n)图形
subplot(3,1,1)stem(x)%绘制x(n)的图形
title(‘理想采样信号序列’)%设置结果图形的标题
(2)绘制信号x(n)的幅度谱和相位谱
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 (‘理想采样信号序列的相位谱’)
(3)改变参数为: 1, 0.4, 2.0734, 1 0 A = α = Ω = T =
n=0:50%定义序列的长度是50
A=1a=0.4w0=2.0734T=1%设置信号有关顷搜的参数和采样率T
x=A*exp(-a*n*T).*sin(w0*n*T)%pi 是MATLAB 定义的π,信号乘可采用“.*”
close all %清除已经绘制的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 (‘理想采样信号序列的相位谱’)
2、单位脉冲序列
在 MatLab 中,这一函数可以用zeros 函数实现:
n=1:50%定义序列的长度是50
x=zeros(1,50)%注意:MATLAB 中数组下标从1 开始
x(1)=1close all
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 (‘单位冲击信号的相位谱’)
3、矩形序列
n=1:50x=sign(sign(10-n)+1)
close allsubplot(3,1,1)stem(x)title(‘单位冲击信号序列’)
k=-25:25X=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、特定冲击串
x(n) =δ (n) + 2.5δ (n ?1) + 2.5δ (n ? 2) +δ (n ? 3)
n=1:50%定义序列的长度是50
x=zeros(1,50)%注意:MATLAB 中数组下标从1 开始
x(1)=1x(2)=2.5x(3)=2.5x(4)=1
close allsubplot(3,1,1)stem(x)title(‘单位冲击信号序列数厅’)
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)