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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】关于matlab循环
71/1返回列表
查看: 947  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
信息提示
关闭
请填处理意见
关闭