请问,滤波后如何出图呢?是否用需用反傅里叶变化?freq_sig1=fft(quake_data(:,4)); lfreq=length(freq_sig1); nyq=50; half_lfreq=floor(lfreq/2); freq=[1:1:(half_lfreq)]/(half_lfreq)*nyq posfreq=abs(freq_sig1); figureplot(freq,posfreq(1:h
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/26 22:23:55
请问,滤波后如何出图呢?是否用需用反傅里叶变化?freq_sig1=fft(quake_data(:,4)); lfreq=length(freq_sig1); nyq=50; half_lfreq=floor(lfreq/2); freq=[1:1:(half_lfreq)]/(half_lfreq)*nyq posfreq=abs(freq_sig1); figureplot(freq,posfreq(1:h
请问,滤波后如何出图呢?是否用需用反傅里叶变化?
freq_sig1=fft(quake_data(:,4));
lfreq=length(freq_sig1);
nyq=50;
half_lfreq=floor(lfreq/2);
freq=[1:1:(half_lfreq)]/(half_lfreq)*nyq
posfreq=abs(freq_sig1);
figure
plot(freq,posfreq(1:half_lfreq));
[b_low,a_low]=butter(8,0.2);
low=filter(b_low,a_low,quake_data(:,3));
[b_high,a_high]=butter(8,0.8,'high')
[b_band.a_band]=butter(8,[0.2 0.4])
freqz(b_low,a_low)
这是我的代码,quake_data里面一共有三个信号波,但是用上面代码滤波后出图结果都一样,麻烦请看看代码应该如何修改.十分感谢了~
请问,滤波后如何出图呢?是否用需用反傅里叶变化?freq_sig1=fft(quake_data(:,4)); lfreq=length(freq_sig1); nyq=50; half_lfreq=floor(lfreq/2); freq=[1:1:(half_lfreq)]/(half_lfreq)*nyq posfreq=abs(freq_sig1); figureplot(freq,posfreq(1:h
如果我没记错,freqz函数是会画图的.所以你最后得到的图就是freqz(b_low,a_low).