马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear all
clc
A=1;m=0.38;B=0.98;C=1;fr=120;fz=3000;fk=80;
w0=2*pi*fr;
w1=2*pi*fz;
w2=2*pi*fk;
t=1/16384;
t1=0:t:0.0625;
x1=A*(1+m*cos(w0.*t1)).*cos(w1.*t1+B.*sin(w0*t1));
x2=sin(w2*t1);
x3=C*randn(size(t1));
x=x1+x2+x3;
[cory,lag]=xcorr(x,'unbiased');
subplot(5,1,1);
plot(lag*t,cory);
axis([0 0.06 -4 4]);
xlabel('时间 t/s');
ylabel('振幅 cm');
[c,l]=wavedec(cory,3,'db4');
a3=wrcoef('a',c,l,'db4',3);
d3=wrcoef('d',c,l,'db4',3);
d2=wrcoef('d',c,l,'db4',2);
d1=wrcoef('d',c,l,'db4',1);
subplot(5,1,2);
plot(d1,'LineWidth',2);
ylabel('d1');
subplot(5,1,3);
plot(d2,'LineWidth',2);
ylabel('d2');
subplot(5,1,4);
plot(d3,'LineWidth',2);
ylabel('d3');
xlabel('时间 t/s');
subplot(5,1,5);
plot(a3,'LineWidth',2);
ylabel('a3');
xlabel('时间 t/s');
==============================================
这个程序是要用小波分析对原信号进行三层的分解,但是执行了之后,原信号的横坐标是(0,0.06)
而在分解之后的信号d1,d2,d3的横坐标却是[0,2100左右],...
是设置参数的时候错了吗..如果我想把所有波形的横坐标都统一为原信号的横坐标,即(0,0.06)应该怎么样修改程序????
搜索更多相关主题的 |