- 积分
- 0
- 注册时间
- 2009-2-23
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2020-5-22 17:00:41
|
显示全部楼层
来自 广东深圳
m=length(Scope1); %scope1,时域响应,1列时间,2列幅值 见附件N=max(m); %【提取 采样点数】
t=Scope1(:,1);
y=Scope1(:,2);
N=2^(log(m)/log(2));
T=1/N; %采样周期
kT=(0:N-1)*T; %固定周期采样
yT=interp1(t,y,kT);
s=fft(yT,N);%【傅里叶变换】
w=(0:N-1)/(N*T);
%【时域】
figure(1);
plot(kT,yT);
title('时域图');
xlabel('时间/s');
ylabel('幅值');
%【频谱图】
figure(2);
plot(w(1:N/2),abs(s(1:N/2)))%由于对称只取一半频谱;
title('频谱图');
xlabel('频率/Hz');
ylabel('幅值');
axis([0 300 0 max(abs(s))]) |
|