微分出的图不太理想，帮忙看下。

cammer534 · 发表于 2008-10-7 14:52

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

fs=10^(-9);
w=0:1/fs:128000000000;
ew=sin((w*0.5*10^(-9)-50))./(w-10^11+eps)+sin((w*0.5*10^(-9)+50))./(w+10^11+eps);
et=fftshift(fft(ew,8192));
et=(abs(et)/fs);
k=0:8191;
figure
plot(fs/8192*k,et);
xlabel('TIME(SECOND)');ylabel('e(t)');
figure
plot(w,ew);xlabel('FREQUENCY(Hz)');ylabel('e(w)');
w=0:1/fs:128000000000-1;
dew=diff(ew)/fs;
figure
plot(w,dew)
dot=fftshift(fft(dew,8192));
dot=(abs(dot)/fs);
k=0:8191;
figure
plot(fs/8192*k,dot);xlabel('TIME(SECOND)');ylabel('dot');

vampire101 · 发表于 2008-10-7 19:22

微分求什么呢？感兴趣问问哈!

账号		自动登录	找回密码
密码			我要加入

[FFT] 微分出的图不太理想，帮忙看下。

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。