微分出的图不太理想，帮忙看下。
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'); 微分求什么呢？感兴趣问问哈!
页:
[1]