| 查看: 611 | 回复: 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, ); 小波变换.jpg |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有85人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB中程序看不懂的简单问题！谢谢！
已经有3人回复
- 高分求将我的MATLAB程序改写为FORTRAN程序（1000BB）
已经有13人回复
- Matlab中的fmincon功能求助
已经有16人回复
- 求在matlab中进行分类的程序
已经有5人回复
- TSVM的matlab程序
已经有7人回复
- matlab新手求助
已经有6人回复
- 有个matlab程序，但是运行不下来，求高手
已经有23人回复
- matlab求助？？？
已经有6人回复
- 模式搜索MATLAB程序
已经有9人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 求助Matlab解方程问题
已经有6人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
- 看不懂matlab程序该怎么办？
已经有7人回复
- 求助matlab编程
已经有13人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】请问matlab编程语言如何打包成为独立的应用程序啊？谢谢
已经有8人回复
- 【求助】请教把java程序转化成matlab程序
已经有9人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【求助】matlab 启动的问题
已经有7人回复
- 【原创】多目标规划matlab程序实现---(原创)
已经有120人回复
shuishui1988
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1234.9
- 散金: 50
- 红花: 5
- 帖子: 128
- 在线: 79.2小时
- 虫号: 2817868
- 注册: 2013-11-21
- 性别: MM
- 专业: 半导体物理
|
程序如下： 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
shuishui1988
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1234.9
- 散金: 50
- 红花: 5
- 帖子: 128
- 在线: 79.2小时
- 虫号: 2817868
- 注册: 2013-11-21
- 性别: MM
- 专业: 半导体物理
3楼2015-05-29 14:29:09
shuishui1988
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1234.9
- 散金: 50
- 红花: 5
- 帖子: 128
- 在线: 79.2小时
- 虫号: 2817868
- 注册: 2013-11-21
- 性别: MM
- 专业: 半导体物理
4楼2015-06-01 09:54:35
投票:
30