| 查看: 1119 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】关于matlab循环
|
|||
|
目前问题是这样的： 阶微分方程组得到一组数据（2500个时间单位） 现取其中一组例如x的序列进行hilbert变换得到数据（定义相位） 现在想把得到的数据有一个判断：当后一个时刻的值比前一个时刻的值大的时候，我们定义的相位+2pi 请问这个循环怎么写！！ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有16人回复
-
系统查不到
已经有10人回复
-
基金不中，共勉
已经有3人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
-
我不理解！
已经有13人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有9人回复
-
2026国自然函评费到账
已经有25人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有28人回复
-
能否退出参与的面上项目解除限项
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】三维矩阵运算
已经有8人回复
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
6楼2010-10-09 18:12:40
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
2楼2010-10-07 12:30:00
|
下面是我写的程序，如何修改? function rosslertest2() clc; clear all; tic; [T,x]=ode45(@Rossler,[0,2000],[0.2;0.3;0.6;0.4;0.7;0.1]); y1=x(:,1); m1=hilbert(y1); N=length(m1); n1=imag(m1); l1=real(m1); p1=atan2(n1,l1); p2=p1; dy1=0; for i=2:N; if p1(i) end p2(i)=p2(i)+dy1; end y2=x(:,4); m2=hilbert(y2); M=length(m2); n2=imag(m2); l2=real(m2); p3=atan2(n2,l2); p4=p3; dy2=0; for i=2:M; if p3(i) end p4(i)=p4(i)+dy2; end p5=abs(p4-p2); plot(p5) toc; function dx = Rossler(t,x) dx1=-(1+0.015)*x(2)-x(3)+0.01*(x(4)-x(1)); dx2=(1+0.015)*x(1)+0.15*x(2); dx3=0.2+x(3)*(x(1)-10); dx4=-(1-0.015)*x(5)-x(6)+0.01*(x(1)-x(4)); dx5=(1-0.015)*x(4)+0.15*x(5); dx6=0.2+x(6)*(x(4)-10); dx=[dx1;dx2;dx3;dx4;dx5;dx6]; |
3楼2010-10-07 19:50:01
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
4楼2010-10-08 09:23:18
投票:
