2021电工杯数学建模A题思路（更新）--最强思路

Posted 2021-06-01 微信公众号：您好啊数模君

思路程序均免费提供给各位参赛者，思路原创，看下你们是不是亏了

 

​

​

 

load('Data.mat')

Fs=50;%频率根据上一问中不同波形填入，这里仅给个案例，这个频率是波的频率不是采样频率

H=4000;%采样个数

x=1:4000;%取数

y=Ia(x);%取数

y0=abs(fft(y,H));%快速傅里叶变换的幅值

y1=fftshift(y0);%幅值修正得到真实幅值

f=x*Fs/H;%将横坐标转化，显示为频率f= n*(fs/N)

figure

subplot(121)

plot(x,y)

subplot(122)%直接FFT结果

plot(f(1:H/2),y0(1:H/2))

xlabel('Frequency');

ylabel('Amplitude');

Z=[f(1:H/2)',y0(1:H/2)];%储存频率及谐波

%以幅值为主进行排序

[~,a]=sort(y0(1:H/2),'descend');%降序

Z=Z(a,:);

Z=[a,Z];

load('Data.mat')
Fs=50;%频率根据上一问中不同波形填入，这里仅给个案例，这个频率是波的频率不是采样频率
H=4000;%采样个数
x=1:4000;%取数
y=Ia(x);%取数
y0=abs(fft(y,H));%快速傅里叶变换的幅值
y1=fftshift(y0);%幅值修正得到真实幅值
f=x*Fs/H;%将横坐标转化，显示为频率f= n*(fs/N)
figure
subplot(121)
plot(x,y)
subplot(122)%直接FFT结果
plot(f(1:H/2),y0(1:H/2))
xlabel('Frequency'); 
ylabel('Amplitude');
Z=[f(1:H/2)',y0(1:H/2)];%储存频率及谐波