如何用matlab画出周期矩形脉冲信号及频谱图
谢谢哈,不过我试了,好像不太对哦,图像不对,而且有错误,没有频谱图,麻烦您再帮下忙哦~ 展开
这样:
x1=ones(1,2);
x2=[x1,zeros(1,6)];
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk;
magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k,magXk);
xlabel('k');ylabel('Xtilde(k)');
title('DFS of SQ.wave :L=2,N=8')%该信号频谱图
扩展资料:
注意事项
1、周期性矩形波(方波)信号:在MATLAB中用square函数来表示,其调用形式为
y=square(t,DUTY)
其作用类似于sin(t),用以产生一个时长为t、幅值为±1的周期性方波信号,其中的DUTY参数表示占空比,即在信号的一个周期中正值所占的百分比。例如频率为30Hz的周期性方波信号的 MATLAB 参考程序如下:
t=-2*pi/100:pi/1024:2*pi/100;
y=square(2*pi*30*t,50);
plot(t,y);
grid
ylim([-1.5 1.5])
2、matlab产生方波脉冲和周期性方波信号
矩形脉冲信号:在MATLAB中用rectpuls函数来表示,其调用形式为
y=rectpuls(t,width)
用以产生一个幅值为1,宽度为width,相对于t=0点左右对称的矩形波信号,该函数的横坐标范围由向量t决定,是以t=0为中心向左右各展开width/2的范围,width的默认值为1。
例:画宽为2的矩形脉冲信号的MATLAB源程序如下:
width=2;
t=-2:0.001:3;
ft=rectpuls(t,width);
plot(t,ft);
grid on;
ylim([-0.5 1.5])
2023-05-25 广告
x2=[x1,zeros(1,6)];
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk;
magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k,magXk);
xlabel('k');ylabel('Xtilde(k)');
title('DFS of SQ.wave :L=2,N=8')%该信号频谱图
不好意思写掉了一点
x1=ones(1,2);
x2=[x1,zeros(1,6)];
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
k1=[-N/2:N/2];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk; magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k1,magXk);
xlabel('k1');ylabel('Xtilde(k)');
title('DFS of SQ.wave :L=2,N=8')%该信号频谱图
你再试一试看看,我运行了是正确的
频谱图就用fft()命令做傅里叶变换,之后再处理一下就行了
具体不多解释,给个程序样例:
(为了省点事,这里采用的是rectpuls函数输出一个矩形脉冲)
t=-20:20;
w=10;
y=rectpuls(t,w);
%矩形脉冲信号
yy=fft(y);
%傅里叶变换
n=size(y);
n=n(2);
%宽度
fy=abs(fft(y))/n*2;
%得到频谱
subplot(2,1,1)
plot(y);
subplot(2,1,2)
plot(fy);
输出见图。
x=10*x2;%所求的周期脉冲信号
N=8;%长度为8
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;Xk=x*WNnk;magXk=abs([Xk(N/2+1:N),Xk(1:N/2+1)]);
subplot(2,1,1);stem(n,x);%画出周期脉冲信号
subplot(2,1,2);stem(k,magXk);
xlabel('k');ylabel('Xtilde(k)');
title('DFS
of
SQ.wave
:L=2,N=8')%该信号频谱图