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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】关于matlab循环
71/1返回列表
查看: 948  |  回复: 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-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的回帖

brqhl_ing

银虫 (小有名气)
引用回帖:
Originally posted by lijinfeng042 at 2010-10-08 09:23:18:

不理解你的意思 你是要实时的改变相位还是求出结果统一加?

我的意思是分别定义出两个相位（瞬时），然后做相位差。画出相位差的图。hilbert变换后相位u一直都在-pi到pi之间变化。但相位应该是随时间增长的 p1和p3做出来的都是从-pi增加到pi然后又突然落到-pi这么个状态。现在想消除这个落差。把他们弄到一条线上。你可以先画出p1p3的图看一下
赞一下 回复此楼
5楼2010-10-08 09:54:37
已阅   回复此楼   关注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的回帖

brqhl_ing

银虫 (小有名气)
引用回帖:
Originally posted by lijinfeng042 at 2010-10-09 18:12:40:
不是很理解 相位 呵呵
    if    p3(i)          dy2=dy2+2*pi;
     end
         p4(i)=p4(i)+dy2;
      
  end
通过这个来?

是通过这个来的啊。matlab能做做一些大的循环？
赞一下 回复此楼
7楼2010-10-11 08:42:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 brqhl_ing 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 生物学调剂 +3 Surekei 2026-03-21 3/150 2026-03-21 18:31 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
信息提示
关闭
请填处理意见
关闭