应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】关于matlab循环
51/1返回列表
查看: 1122  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

brqhl_ing

银虫 (小有名气)

[交流] 【求助】关于matlab循环

目前问题是这样的：
阶微分方程组得到一组数据（2500个时间单位） 现取其中一组例如x的序列进行hilbert变换得到数据（定义相位） 现在想把得到的数据有一个判断：当后一个时刻的值比前一个时刻的值大的时候，我们定义的相位+2pi   请问这个循环怎么写！！
回复此楼

» 猜你喜欢
1楼 2010-10-06 19:58:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab
★
ctgu_zheng(金币+1):欢迎继续交流讨论 2010-10-09 23:42:37
不是很理解 相位 呵呵
    if    p3(i)          dy2=dy2+2*pi;
     end
         p4(i)=p4(i)+dy2;
      
  end
通过这个来?
赞一下(1人) 回复此楼
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
6楼2010-10-09 18:12:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

lijinfeng042

木虫 (小有名气)

Matlab
★ ★
ctgu_zheng(金币+1):鼓励继续交流应助。。。 2010-10-07 20:25:16
hiqun(金币+1):感谢应助 2010-10-08 11:53:49
取diff >0的用2pi取代 加回原数组
赞一下(2人) 回复此楼
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
2楼2010-10-07 12:30:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brqhl_ing

银虫 (小有名气)
引用回帖:
Originally posted by lijinfeng042 at 2010-10-07 12:30:00:
取diff >0的用2pi取代 加回原数组

下面是我写的程序，如何修改?
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)          dy1=dy1+2*pi;
     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)          dy2=dy2+2*pi;
     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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab
★
hiqun(金币+1):感谢应助 2010-10-08 09:36:13
引用回帖:
Originally posted by brqhl_ing at 2010-10-07 19:50:01:

下面是我写的程序，如何修改?
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(m ...

不理解你的意思 你是要实时的改变相位还是求出结果统一加?
赞一下(1人) 回复此楼
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
4楼2010-10-08 09:23:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 13:08 by winsaint
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 5/250 2026-08-28 12:41 by Tide man
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
信息提示
关闭
请填处理意见
关闭