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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】关于matlab循环
71/1返回列表
查看: 940  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 085410人工智能专硕317求调剂（0854都可以） +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
信息提示
关闭
请填处理意见
关闭