马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- clc;
- % load L20X.txt;
- % y = L20X;
- load M1;
- load L1;
- figure(1)
- subplot(2,1,1),plot(M1);
- subplot(2,1,2),plot(L1);
- N = size(M1);
- figure(2)
- % n1 = find(y<=min(y));
- % plot(abs(fft(y)));
- % y1 =ones(N(1,1),1);
- % a = y(1);
- % y1 = a.*y1;
- % y = y - y1;
- b = fir1(100,0.2,0.3,blackman(101));
- m1 = filter(b,1,M1);
- l1 = filter(b,1,L1);
- subplot(2,1,1),plot(m1);
- subplot(2,1,2),plot(l1);
- % x = x + y1;
- % n2 = find(x<=min(x));
- % n = n2-n1;
- % figure(3);
- % % plot(abs(fft(x)));
- % plot(x);
- % save y.txt y -ASCII
- p_a = find(l1>=max(l1))
- p_b = find(l1<=min(l1))
- figure(3)
- plot([l1(31:(p_b-p_a+1)*2); l1(31)],[m1(31:(p_b-p_a+1)*2); m1(31)])
- hold on
复制代码
[ 本帖最后由 suffer 于 2006-10-9 20:34 编辑 ] |