matlab矩阵二阶微分方程求解,MATLAB解含参数方程、矩阵方程、二阶微分方程组
一、以下公式: 其中x、y、m、n为参数,a、b为未知数,利用MATLAB求解方程;spa
syms x y m n a b;code
[a,b]=solve('x=m*cos(a)+n*cos(a+b)','y=m*sin(a)+n*sin(a+b)','a','b');blog
会获得四组解;对a(3)、b(3)、a(4)、b(4)进行化简,it
simplify(a(3));class
simplify(b(3));im
simplify(a(4));db
simplify(a(4));img
获得:di
a(3)= 2*atan((2*m*y - (- m^4 + 2*m^2*n^2 + 2*m^2*x^2 + 2*m^2*y^2 - n^4 + 2*n^2*x^2 + 2*n^2*y^2 - x^4 - 2*x^2*y^2 - y^4)^(1/2))/(m^2 + 2*m*x - n^2 + x^2 + y^2));dba
b(3)= 2*atan(((- m^2 + 2*m*n - n^2 + x^2 + y^2)*(m^2 + 2*m*n + n^2 - x^2 - y^2))^(1/2)/(- m^2 + 2*m*n - n^2 + x^2 + y^2));
a(4)= 2*atan((2*m*y + (- m^4 + 2*m^2*n^2 + 2*m^2*x^2 + 2*m^2*y^2 - n^4 + 2*n^2*x^2 + 2*n^2*y^2 - x^4 - 2*x^2*y^2 - y^4)^(1/2))/(m^2 + 2*m*x - n^2 + x^2 + y^2));
b(4)= -2*atan(((- m^2 + 2*m*n - n^2 + x^2 + y^2)*(m^2 + 2*m*n + n^2 - x^2 - y^2))^(1/2)/(- m^2 + 2*m*n - n^2 + x^2 + y^2));
以a(3)为例,假设x=20,其它参数未知,则:
a(3)=subs(a(3),x,20) ; 便可将x取值替换为20;
假设 x=20 , y=30 ,其它参数未知,则:
a(3)=subs(a(3),{x,y},{20,30});
二、求解矩阵方程:
a=[1 0;0 2];
b=[1;3];
syms x1 x2;
y=a*[x1;x2]+b;
s=solve(y(1),y(2),'x1','x2');
s.x1 %w2值
s.x2 %w2值
三、
Matlab求解二阶微分方程组:
m=1;
g=9.8;
k=100;
l0=1.1;
x0=[0.1 0 0 0]; %初始值;
% 定义 x(1)=l, x(2)=l', x(3)=a, x(4)=a';
dx=@(t,x)[x(2); (m*x(1)*x(4)^2-m*g*cos(x(3))-k*(x(1)-l0))/m; x(4); (-2*m*x(1)*x(2)*x(4) + m*g*x(1)*sin(x(3)))/(m*x(1)^2)];
[t,x]=ode15s(dx,[0 10],x0);
lstr = { '\itl', '{\itl}''', '\alpha', '\alpha''' };
for i=1:length(lstr)
subplot(2,2,i)
plot(t, x(:,i));
xlabel('Time')
ylabel( lstr{i} )
end
matlab矩阵二阶微分方程求解,MATLAB解含参数方程、矩阵方程、二阶微分方程组相关推荐
- 微分方程求解 matlab,4MATLAB常微分方程求解.ppt
4MATLAB常微分方程求解 MATLAB微分方程 1 求简单微分方程的解析解 2 求微分方程的数值解 3 建模实例 1 求简单微分方程的解析解 求微分方程(组)的解析解命令: dsolve('方程1 ...
- 【AFSA TSP】基于matlab人工鱼群算法求解旅行商问题【含Matlab源码 422期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[TSP]基于matlab人工鱼群算法求解旅行商问题[含Matlab源码 422期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2 ...
- matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组
1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...
- 微分方程求解matlab冲激信号,用Matlab求冲激响应的几种方法
·实验技术与方法· 用 Matlab 求冲激响应的几种方法 贺富堂,应柏青,张 锋,孙 敏 ( 西安交通大学 电工电子教学实验中心,陕西 西安 710049) 摘 要: 用 Matlab 软件求冲激响 ...
- 二阶微分方程的matlab解法,以动力学方程为例
过冷水最近有接触一点点动力学的知识.作为动力学入门,当然的会解动力学方程了.于是本期过冷就教大家解动力学微分方程. 上图是两个小车通过弹簧链接起来的做来回摆动运动.应用拉克朗日方程建立系统的运动微分方 ...
- matlab 矩阵与数比较,MATLAB 对矩阵中的数据进行大小比较
用matlab对矩阵进行高斯消元 %求解例题2.1%高斯法求解线性方程组Ax=b%A为输入矩阵系数,b为方程组右端系数%方程组的解保存在x变量中%先输入方程系数A=[1 2 3;2 7 5;1& ...
- matlab矩阵处理实验报告,matlab实验报告一二三
<matlab实验报告一二三>由会员分享,可在线阅读,更多相关<matlab实验报告一二三(37页珍藏版)>请在金锄头文库上搜索. 1.2015 秋秋 2013 级级MATLA ...
- matlab矩阵的第一列,matlab提取矩阵第一列
第2章 MATLAB矩阵及其运算 2.1 变量和数据操作 2.2 MATLAB矩阵 2.3 MATLAB运算 2.4 矩阵分析 2.5 矩阵的超越函数 2.6 字符串 2.7 结构数据和单元数据 2. ...
- MATLAB矩阵乘法、MATLAB矩阵的转置、MATLAB串联矩阵
MATLAB矩阵乘法 MATLAB中如果有两个矩阵 A 和 B ,其中 A 是 m*n 矩阵,B 是 n*p 矩阵,那么他们相乘能够产生一个 m*p 的矩阵 C. MATLAB矩阵乘法只发生在矩阵 A ...
- matlab 矩阵与向量乘,matlab问题,关于矩阵向量的乘问题
问题描述: matlab问题,关于矩阵向量的乘问题 mx=wtbx*T; my=wtby*T; mz=wtbz*T;%%%通过四元素解微分方程 p=mx.^2+my.^2+mz.^2;//这个地方有问 ...
最新文章
- Linux终端C语言实现图片拷贝
- 加载BeanFactory
- Mysql学习之order by的工作原理
- Jmeter接口测试进阶
- Python网络数据采集2-wikipedia
- Express框架Restful API Ajax 跨域 开启Cookie支持
- mac docker安装linux,Mac上使用docker安装centos
- tf.name_scope()详解【命名空间其实就是给几个变量包一层名字,方便变量管理】
- android compress函数,Linux Kernel(Android) 加密算法小结(cipher、compress、digest)
- OculusQuest2 直接连电脑调试
- 3D图像生成和编辑研究成果大放送!朱俊彦团队放出两篇论文实现代码 | 资源...
- nessus8.9.0百度网盘_免费百度网盘SVIP共享20.1.19
- 如何优雅的调用第三方接口
- LoadRunner测试139邮箱的脚本
- QT界面主题风格设置
- cmf+ThinkPHP6.0 配置使用多应用路由
- HTML、js实现图片绕中心旋转
- matlab设置固定的窗宽窗位,如何设定窗宽窗位,附正常人体组织CT值
- 编译原理——自下而上的语法分析方法(LR分析法)
- Android的虚拟设备的缩写,Android虚拟设备的英语缩写是
热门文章
- python 获取像素颜色_python如何读取像素值
- 工作中ibatis中的连表查询及in()的使用案例
- ADMM随堂笔记(3):一致和协同
- 计算机网络图论,网络图论
- 网易云音乐Android一面面经
- 自动控制原理学习笔记(二)线性定常连续控制系统的数学模型
- A星寻路算法详解(完整代码+图片演示)
- igs无法分配驱动器映射表_明白步进驱动器原理,就知道为什么步进电机要加驱动器了...
- 【干货】区块链技术生态的设计|《白话区块链》作者蒋勇分享实录
- APP原型设计利器-墨刀MockingBot