sa(t)信号对应的Matlab程序怎么编写?

sa(t)信号对应的Matlab程序怎么编写?,第1张

程序如下:

t=(-5*pi:100/pi:5*pi)

x=sinc(t/pi)

plot(t,x)

xlabel('t')

ylabel('x(t)')

title('Sa(t)')

gridon

抽样信号也被称为抽样函数或Sa(t)函数,是指sint与t之比构成的函数。

MATLAB是美国MathWorks公司出品的商业数学软件,用于数据分析、无线通信、深度学习、图像处理与计算机型帆视觉、信号处理、量化金融与风险管理、机器人,控卜誉雹制系统等领域虚链。

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(‘单位冲击信号序列数厅’)

Y=fft(y,1024)

pyy=Y.*conj(Y)/1024

这两句讲的是对随机信号y的功率谱估计,可以参考数字信号处理中功率谱估计直接法部分,这个是定义式,要说由来要看 胡广书著《数字信号处理理论、算法与实现》的10.4节、11.2.1节,其他相关函数、平稳随机信号的功率谱的知识也可以看看便于理解

f=1000*(0:512)/1024

这句讲的是构建真实的频率轴,这个频率轴烂磨灶是0Hz ~ 512*(1000/1024)Hz,要理解这个就要理解DFT变换的定义的推导过程,简单的说对时域进行游派以Ts为间隔的采样相当于对频谱进行fs=1/Ts的周期延拓,反之对频域进行以fs为间隔的采样相当于对时域波形进行Ts=1/fs的周期延拓,这个程序中Ts=0.001,对应fs=1000,也就是频谱一个周期长度为1000Hz,由于进行的是1024点的fft就相当于频率采样间隔为1000/1024,再来,程序中只想画出正频段部分因此就只去了0~512这几个点,DFT知识可以参考胡广书著《数字信号处理饥扮理论、算法与实现》3.5节

求采纳为满意回答。


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/12317358.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-24
下一篇 2023-05-24

发表评论

登录后才能评论

评论列表(0条)

保存