24小时热门版块排行榜    

查看: 329  |  回复: 3

shuishui1988

铜虫 (小有名气)

[求助] Matlab 程序求助

本人初学Matlab,最近在研究小波变换,下面的程序得到的小波变换的波形非常不好,不知道问题出在哪里,希望高手解答,金币奉上~
clear;
close all;
clc

m=10^(-17);  
k=10;      
E0=1*1.60217646*10^(-19);
%r=2*sqrt(k*m);
r=10^(-3);%6.32*10^(-4);
beta=30*pi/180;
vm=400;
ax=4;
ay=7;
jiange=0.002;
tend=75;
x0=[0 0];
[t x]=ode45(@ode1DMoS2,[0:jiange:tend],x0,[],m,k,E0,r,vm,ax,ay,jiange,beta);
xx=x(:,1);
vx=x(:,2);
counter=1;
for i=1:length(xx)

          b(counter)=xx(i);
          f(1,counter)=k*(vm*t(i)*jiange-b(counter))*0.1;
          f(2,counter)=vm*t(i)*jiange*0.1;
          counter=counter+1;
end
figure;
plot(f(2,,f(1,);
figure;
[c,l]=wavedec(f(1,,5,'db9');
a4=wrcoef('a',c,l,'db9',5);
a3=wrcoef('d',c,l,'db9',4);
a2=wrcoef('d',c,l,'db9',3);
a1=wrcoef('d',c,l,'db9',2);
d1=wrcoef('d',c,l,'db9',1);
figure(2);
subplot(5,1,1);
plot(a4,'k-','LineWidth',2);
ylabel('a4');
subplot(5,1,2);
plot(a3,'k-','LineWidth',2);
ylabel('a3');
subplot(5,1,3);
plot(a2,'k-','LineWidth',2);
ylabel('2');
subplot(5,1,4);
plot(a1,'k-','LineWidth',2);
ylabel('a1');
subplot(5,1,5);
plot(d1,'k-','LineWidth',2);
ylabel('d1');
xlabel('时间 t/s');
subplot(6,1,1);plot(f(1,);

Matlab 程序求助
小波变换.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuishui1988

铜虫 (小有名气)

程序如下:
clear;
close all;
clc

m=10^(-17);  
k=10;      
E0=1*1.60217646*10^(-19);
%r=2*sqrt(k*m);
r=10^(-3);%6.32*10^(-4);
beta=20*pi/180;
vm=400;
ax=5.7;
ay=9.88;
jiange=0.002;
tend=75;
x0=[0 0];
[t x]=ode45(@ode1DMoS2,[0:jiange:tend],x0,[],m,k,E0,r,vm,ax,ay,jiange,beta);
xx=x(:,1);
vx=x(:,2);
counter=1;
for i=1:length(xx)

          b(counter)=xx(i);
          f(1,counter)=k*(vm*t(i)*jiange-b(counter))*0.1;
          f(2,counter)=vm*t(i)*jiange*0.1;
          counter=counter+1;
end
figure;
plot(f(2,,f(1,);
figure;
[c,l]=wavedec(f(1,,5,'db9');
a4=wrcoef('a',c,l,'db9',5);
a3=wrcoef('d',c,l,'db9',4);
a2=wrcoef('d',c,l,'db9',3);
a1=wrcoef('d',c,l,'db9',2);
d1=wrcoef('d',c,l,'db9',1);
figure(2);
subplot(5,1,1);
plot(a4,'k-','LineWidth',2);
ylabel('a4');
subplot(5,1,2);
plot(a3,'k-','LineWidth',2);
ylabel('a3');
subplot(5,1,3);
plot(a2,'k-','LineWidth',2);
ylabel('2');
subplot(5,1,4);
plot(a1,'k-','LineWidth',2);
ylabel('a1');
subplot(5,1,5);
plot(d1,'k-','LineWidth',2);
ylabel('d1');
xlabel('时间 t/s');
subplot(6,1,1);plot(f(1,);
%figure;plot(t,f(1,);
%figure;plot(t,b);
%figure;plot(t,f(2,);
%figure;plot(t);
2楼2015-05-29 14:27:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuishui1988

铜虫 (小有名气)

笑脸实际是程序“:)”,希望高手帮帮忙啊!
3楼2015-05-29 14:29:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuishui1988

铜虫 (小有名气)

求助各位,帮帮忙啊~
4楼2015-06-01 09:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shuishui1988 的主题更新
信息提示
请填处理意见