%去均值 for p=1:89 ave=mean(X(:,p)); X(:,p)=X(:,p)-ave; for o=1:160 XX((p-1)*160+o)=X(o,p); end end figure;plot(1:14240,XX,'k'); axis([5000 5400 ,-inf,inf]);xlabel('时间(s)');ylabel('时间波形');title('展开信号'); hold on; b1=[0.008233 -0.004879 0.007632 -0.004879 0.008233]; a1=[1 -3.6868 5.8926 -5.0085 2.2518 -0.4271]; [H,w]=freqz(b1,a1,512); for r=1:89 sf_1(:,r)=filter(b1,a1,X(:,r)); end plot(1:14240,sf,'b');xlabel('时间(s)');ylabel('幅值');title('椭圆滤波');axis([5000 5400 ,-inf,inf]) %改变横轴范围 hold on; |
最佳答案
2 条回复