clc;

source=simout1;
filter=simout2;

%%??原始數據
subplot（311）;
plot（source.Timesource.Data‘b‘filter.Timefilter.Data‘r‘）;
title（‘數據‘）;
legend（‘原始數據‘‘濾波之后數據‘）;
grid?on;

%%?fft?
Ts=1000;
x=source.Data;
N=source.Length;
n=1:N;
n_2=1:N/2;
f=n*Ts/N;
y=fft（xN）;
mag=abs（y）;
subplot（312）;
plot（f（n_2）?mag（n_2））;
title（‘原始數據fft‘）;
axis（[0?Ts/2?0?3*1e6]）;
grid?on;

%%?fft?
Ts=1000;
x=filter.Data;
N=filter.Length;
n=1:N;
n_2=1:N/2;
f=n*Ts/N;
y=fft（xN）;
mag=abs（y）;
subplot（313）;
plot（f（n_2）?mag（n_2））;
title（‘濾波之后數據fft‘）;
axis（[0?Ts/2?0?3*1e6]）;
grid?on;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????631??2018-11-27?23:07??fft_analysis.m
?????文件??????352792??2018-11-27?23:03??simout1.mat
?????文件??????793538??2018-11-27?23:03??simout2.mat