Matlab相图

最新推荐文章于 2023-10-12 09:01:33 发布
最新推荐文章于 2023-10-12 09:01:33 发布
阅读量1k 收藏 5
点赞数
文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_46129398/article/details/125419417
版权
此博客展示了通过数值模拟研究水平动力学系统的相图，使用ode45函数对不同初始角度（-4.5到4.5）和常数K=5、R=0.5下的角度变化进行分析。图形展示了角度与角度变化量随时间的演化，有助于理解系统的动态行为。
摘要由CSDN通过智能技术生成

%% 水平动力学运动的相图

clc;clear;figure('numbertitle','off','name','水平动力学运动的相图');

global K R

K=5;R=0.5;t1=[-10:0.1:10];

for i=[-0.5 -1.5 -3 -4.5 0.5 1.5 3 4.5]

    H1=[0;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[0.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    title('K=5,R=0.5');

    hold on

    H1=[-0.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[1.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[-1.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[3;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[-3;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[4.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

    H1=[-4.5;i];

    [t,H] = ode45(@dH,t1,H1);

    plot(H(1:end,1),H(1:end,2));

    axis([-5 5 -4 4]);

    xlabel('角度');ylabel('角度变化量');

    hold on

end

function dHdt=dH(t,H)

global K R 

dHdt=zeros(2,1);

dHdt(1)=H(2);

dHdt(2)=-sin(H(1))*(1-K*cos(H(1)))-R*H(2);

end

d757a5e0ecbe41d49968a267e80ad876.jpg

 

计算物理
关注
  • 0
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
MATLAB常见非线性可视化绘制方法-分岔图与庞加莱截面（混沌可视化、Poincare截面、Logistic、Henon、Lorenz、Rossler、Duffing系统）
m0_67394360的博客
06-06 3237
本文首发于 matlab爱好者 微信公众号，欢迎关注。惯例声明：本人没有相关的工程应用经验，只是纯粹对相关算法感兴趣才写此博客。所以如果有错误，欢迎在评论区指正，不胜感激。本文主要关注于算法的实现，对于实际应用等问题本人没有任何经验，所以也不再涉及。本文大概介绍一下常用的非线性可视化的方法，由于篇幅问题，分了两章。第一篇介绍一下相图与相空间的概念。第二篇介绍一下庞加莱截面和分岔图的绘制方法。后面第三篇就不打算写可视化了，看看要不要写一个频谱分析、李雅普诺夫指数分析、维数分析之类的时间序列分析，但这一部分有趣
MATLAB画相图 compass函数
echo_xiao9的博客
06-09 1979
电路实验需要画三相电源相图 MATLAB 代码 Uun=220;//相量 Uwn=220*exp(1j*-120/180*pi); Uvn=220*exp(1j*120/180*pi); V=[Uun,Uwn,Uvn]; figure(1); compass(V); hold on; Iu=Uun*(j*2*pi*50*4*10^(-6))*300; Iv=80/Uvn*300; Iw=80/Uwn*300; In=Iu+Iv+Iw; W=[Iu,Iv,Iw,In]; compass(W,'r');
MATLAB绘制微分方程的相图/方向场/向量场
热门推荐
陌雨’的博客
02-09 2万+
MATLAB 笔记
MATLAB图形操作3--相图
m0_47017197的博客
11-28 2324
【代码】MATLAB图形操作3--相图。
如何画matlab的相图
最新发布
shddjdjn的博客
10-12 723
h=streamslicetitle('相图')`
matlab画相图
JISANSAN的博客
06-06 1万+
用matlab自带函数streamslice画极限环 [x,y]=meshgrid(linspace(-2,2)); streamslice(x,y,y-x.*(x.^2+y.^2-1),-x-y.*(x.^2+y.^2-1)); title('Limit Circle') xlabel('x');ylabel('y');
matlab画相图离散,【原创】给大家一个画相图的MATLAB程序
weixin_30125497的博客
03-16 1211
CODE:function phaseplot(varargin)type=varargin{end};switch typecase'2't=varargin{1};%温度y_1=varargin{2};%液相线y_g=varargin{3};%气相线y2i=linspace(0,1);ti=spline(y_1,t,y2i);plot(y_1,t,'o',y2i,ti);hold ony1i=...
lorenz混沌系统matlab程序,Lorenz混沌系统相图
weixin_33318222的博客
04-07 1110
来源： http://blog.csdn.net/superdont建立函数1：function xdot=lorrenz(t,x)for t=0:0.1:100;if(cos(5.3*t)>=0)p=1;elsep=-1;endendxdot=[-8/3*x(1)+p*x(2)*x(3);(-x(2)+x(3))*(25-10*cos(5.3*t));-p*x(1)*x(2)-(17.5*...
matlab蔡氏电路相图
07-08
基于matlab的简单蔡氏电路二维相图。
MATLAB | 如何用MATLAB如何绘制各式各样精致的三元相图(ternary plot)
slandarer的博客
04-18 4841
~
系统相图matlab代码-geneus:通过Usadel扩散方程模拟一般非平衡超导体/铁磁体系统
06-08
系统相图matlab代码GENEUS：一般非均衡 Usadel 求解器 GENEUS 是一组数值程序，用于求解一维超导纳米结构中的 Usadel 扩散方程，无论是处于平衡状态还是非平衡状态。 这些程序非常灵活，能够处理具有例如磁性元件、自旋轨道耦合、自旋相关散射、强极化磁界面、电压梯度和温度梯度的系统。 该套件还包含用于计算所有这些混合结构的临界温度和相图的专门程序。 这些程序被构建为用户友好：它们使用简单的配置文件进行配置，其中可以包括数学表达式来初始化物理系统，并且输出很容易导入到 Gnuplot 或 Matlab 等中。 最后，代码是用现代面向对象的 Fortran 编写的，使其既简单又高效。 该软件是由他在博士研究期间开发的。 该研究由挪威 NTNU 监督。 该软件本身是在 ; 这基本上意味着你可以自由地将它用于任何目的，只要你在适当的时候给予信任。 该项目还依赖于外部库 、 和 。 如源文件本身所示，这些可以在混合开源许可证下使用。 为方便起见，这些库已与此软件页面捆绑在一起，它们的来源位于 . 有关更多信息，请参阅 .
达芬振子相图分叉图.rar_matlab相图_ode设置相图_分叉图matlab_分叉图程序_达芬 分叉图
07-14
达芬振子matlab编程相图分叉图庞加莱映射主程序ODE45
matlab代码图-ternplot:在Matlab中绘制三元图
05-26
matlab代码图Matlab的三元图 此代码可帮助您在三元相图上绘制三元相数据。 此包装也可提供。 它的作用类似于绘图（响应保持等），还包括ternlabel.m以标记所有三个轴。 还具有绘制三维图和轮廓的功能。
用Matlab绘制相平面图
Ltter1的博客
03-28 1万+
MATLAB相图程序 Eg1: clc;clear;close; [x,y]=meshgrid(linspace(-6,6)); %定义[x,y]网络采样点的范围 h=streamslice(x,y,y,-2.*x); %描述系统坐标的微分方程 title('极限环的相平面图') %定义标题范围 xlabel('x','Color',[0.6 0.5 0.4]);ylabel('y','C
matlab的给平面图,matlab画相平面图
weixin_35651578的博客
03-17 2178
输出函数的二维相平面图 odephas3 ODE 输出函数的三维相空间图 odeplot ODE 输出函数的时间轨迹图 odeprint 在 Matlab 指令窗显示结果 odeset 创建或改写 ODE ......的相平面轨迹、 本例演示:(A)SIMULINK模型和MATLAB指令的配合使用...(0)=2.8,0i ( 0 ) = i0 , s ( 0 ) = s 0 在相平面 s ~ i...
matlab 绘 洛伦兹系统 3D相图
颹蕭蕭
07-17 1万+
洛伦兹系统 相图
振动系统庞加莱截面、分岔、相平面图matlab程序
qq_41621140的博客
02-19 5478
function ydot = zjzdfun(t,y,flag,u,x0,w0,v,w) ydot = [ y(2); u*(x02-y(1)2)y(2)-y(1)w0^2-vcos(wt) ]; end %求解庞加莱截面 u = [0.85 1.02 0.66 1.08]; x0 =1; w0 = 1; v = 1; w = 0.44; T = 2pi/w; str{1}=‘庞加莱截面-周期1吸引子’; str{2}=‘庞加莱截面-周期2吸引子’; str{3}=‘庞加莱截面-不变环面吸引子’; str
【MATLAB基础绘图第15棒】绘制三元相图
WW、forever的博客
04-18 2488
MATLAB绘制三元相图
用MATLAB编程正弦稳态相量图,matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析...
weixin_42514627的博客
03-18 405
matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析 课程设计任务书 学生姓名： 专业班级： 指导教师： 刘 新 华 工作单位：信息工程学院 题 目: 利用MATLAB对线性电路正弦稳态特性分析 要求完成的主要任务: 1.利用MATLAB强大的图形处理功能，符号运算功能和数值计算功能，实现线性电路正弦稳态特性的仿真波形。 2.绘出表现其正弦稳态特性的时域波形图和向量图。 课程设计的目...
2020-10-04 Matlab中相图的画法
mingtao379的博客
10-04 5553
Matlab中相图的画法Matlab中相图的画法演示结果 Matlab中相图的画法 最近在Youtube中看到一段很好的相图代码，把它拿出来分享给大家。 [t,y]=meshgrid(0:.2:4,0:.2:4); m=t+sin(y); L=sqrt(1+m.^2); quiver(t,y,1./L,m./L) axis tight hold on f=@(u,v) u+sin(v) [u,v]=ode45(f,[0 5],1:3); plot(u,v) hold off axis([0 4 0 4])
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值