以下是一段MATLAB代虚并码,嫌拦可以画出周期为4,脉冲宽度为1的矩形脉冲信号的频谱图:
T = 4 % 周期为4
tau = 1 % 脉冲宽度为1
% 生成矩差者迹形脉冲信号
t = linspace(0, T, 1000) % 时间范围为0到T,分1000个点
x = rectpuls(mod(t, T), tau)
% 计算信号的傅里叶变换
X = fft(x)
% 计算频率范围
N = length(x)
f = (0:N-1) / N / (t(2)-t(1))
% 绘制频谱图
plot(f, abs(fftshift(X)) / N)
title('频谱图')
xlabel('频率')
ylabel('振幅')
运行这段代码可以得到矩形脉冲信号的频谱图。其中,x轴表示频率,y轴表示振幅。可以看到,在频率为0处有一个峰,其余频率处振幅都很小。这是因为矩形脉冲信号的频谱是一个sinc函数,它在频率为0处有一个峰,其余频率处振幅都很小。
可以参考下面的代码:
t=-10*pi:0.01*pi:10*pi%这个由采样频率而定,这个是200的采样频率。
x=sin(2*pi*50*t)+0.5*cos(2*pi*27*t)
f=fft(x,2001)
y=abs(f)%取模,整合洞物幅频和相频
plot(y,'r-')
扩展资料:
Matlab常用函数和命令纳禅液
normcdf 正态分布累计概率密度函数
normest 估计矩阵2范数
norminv 正态分布逆累计概率密度函数
normpdf 正态分布概率密度函数
semilogx: x轴为对数刻度,y轴为线性刻度
semilogy: x轴为线性刻度,袭祥y轴为对数刻度
odephas2 ODE 输出函数的二维相平面图
odephas3 ODE 输出函数的三维相空间图
参考资料来源:百度百科-MATLAB
1、假设该信号的域被借出并且向量t是N维向量,则该信号离散化的采样周期是Ts = 1 / fs = lent /(N-1),其中fs是采样频率渗搏。
2、从上面的离散傅立叶公式,我们可以知道在使用fft函数之后,我们仍然得到一个N维向量。
3、频域的长度定义型蠢为lenf,其满足lent =(N-1)* Ts,solid lenf =((N-1)*(N-1)/ N)/ lent = fs *(N-1) )/ N.当N很大时,lenf = fs。
4、到目前为止在频域中获得的矢量是f =(0:N-1)/(N-1)* lenf =(0:N-1)* fs / N.根据奈奎斯特采样定理,采样频率高于信号中最高频率的两倍,因此不会失真。
5、因此,频谱中的(fs / 2,fs)是无用的频谱,因丛租祥此域矢量和频域信号的范围矢量应该减半。即:f = f(1,N / 2)y = abs(fft(x))y = y(1. N / 2)好吧,我们可以绘制信号的时域和频域图像。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)