陷波器程序如下,麻烦帮忙看看
clc;
clear all;
close all;
x = 0 : 0.05 : 10;
s= 210 * exp(-0.25 * x) .* cos(2 * pi * 1.5 * x + pi) + 110 * exp(0.25 * x) .* cos(2 * pi * 1.0 * x + pi/2 );
%load data;
%s=cc;
N=length(s); % 信号长度
fs=1000; % 采样频率
n=1:N;
n2=1:N/2;
tt=(n-1)/fs; % 时间刻度
ff=(n2-1)*fs/N; % 频率刻度
X=fft(s); % 谱分析
figure(1)
subplot 311;
plot(tt,s,'k');
title('原始数据'); xlabel('时间/s'); ylabel('幅值');
X=X/max(abs(X));
subplot 312;
plot(ff,abs(X(n2)),'k');
axis tight; title('谱分析');
xlabel('频率/Hz'); ylabel('幅值&#