Matlab小程序

微信公开课小程序热门讨论「Matlab小程序」最新回复-刚好这里有一个现成的%自己编程剔除数据异常值clearoffclcX=[13574638471235717467853];ls=length(X);K=3;figure(1);plot(X);fori=1:3;A1(i)=abs(X(i)(X(i+3)+X(i+2)+X(i+1))/3);Y1=[abs(X(i+3)X(i+2));abs(X(i+3)X(i+1));abs(X(i+1)X(i+2))]Y(i)=max(Y1)if(A1(i)/Y(i)K)X(i)=(X(i+1)+X(i+2)+X(i+3))/3;elseendendfori=ls2:ls;A2(i)=abs(X(i)(X(i1)+X(i2)+X(i3))/3);Z1=[abs(X(i3)X(i2));abs(X(i3)X(i1));abs(X(i1)X(i2))]Z(i)=max(Z1)if(A2(i)/Z(i)K)X(i)=(X(i3)+X(i2)+X(i1))/3;elseendendfori=4:ls3;A1(i)=abs(X(i)(X(i3)+X(i2)+X(i1))/3);A2(i)=abs(X(i)(X(i+1)+X(i+2)+X(i+3))/3);Y1=[abs(X(i3)X(i2));abs(X(i2)X(i1));abs(X(i3)X(i2))]Y(i)=max(Y1)Z1=[abs(X(i+3)X(i+2));abs(X(i+3)X(i+1));abs(X(i+1)X(i+2))]Z(i)=max(Z1)if(A1(i)/Y(i)K&A2(i)/Z(i)K)X(i)=(X(i3)+X(i2)+X(i1)+X(i+1)+X(i+2)+X(i+3))/6;elseendendfigure(2);plot(X);再来一个简单的clcclearall;loadXK=1;N=length(X);x=X;figure(1);plot(x);fori=4:N3;if(x(i)x(i1)K&x(i)x(i2)K&x(i)x(i3)K)x(i)=(x(i3)+x(i2)+x(i1)+x(i+1)+x(i+2)+x(i+3))/6;endendfigure(2);plot(x);
clearoffclcbeta=0.1102*(83-8.7);dw=(0.7-0.5)*pi;N=(83-8)/(2.286*dw);N=ceil(N);ifmod(N,2)==0N=N+1endwn=kaiser(N,beta)figure(1);stem(wn);gridon;alfa=(N)/2;wc=(0.7+0.5)*pi/2;n=0:N-1;hd=sin(wc*(n-alfa))./(pi*(n-alfa));hn=hd.*wn;figure(2);stem(hn);gridon;[M,P]=freqz(hn,1,512,whole);ma=abs(M);mal=20*log10(ma);figure(3)plot(P/pi,mal)gridonan=angle(M);figure(4)plot(P/pi,an);gridon 更多有关「Matlab小程序」的疑问请扫码关注微信公开课+小程序！了解更多