如何用matlab解出如下微分方程组并作图?

dφ/dZ=1/R[n+(n+1)(2cosφ+cosφ*cosφ);dR/dZ=-1/tanφ;n为1到2之间的常数得到Z和R的关系并作图,最好附上过程,谢谢。... dφ/dZ=1/R[n+(n+1)(2cosφ+cosφ*cosφ);
dR/dZ=-1/tanφ;
n为1到2之间的常数
得到Z和R的关系并作图,最好附上过程,谢谢。
展开
匿名用户
2014-03-23
展开全部

代码供参考:

n = 1.5;
% x1 = φ, x2 = R
dx = @(z,x)[1/(x(2)*(n+(n+1)*(2*cos(x(1))+cos(x(1))*cos(x(1))))); -1/tan(x(1))];
[t,x]=ode45(dx,[0 1],[0.5 0.5]);
plot(t,x(:,2))

 

你现在的问题有很多不清楚的地方:

1、第一个表达式很乱,我不确定是不是正确理解了你的意图;

2、没给初始条件;

3、常数n是多大?这个微分方程似乎属于刚性的,不同的n可能会有影响,而使用不同的ode求解器结果可能也不一样。

富港检测技术(东莞)有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共振点会同时激发... 点击进入详情页
本回答由富港检测技术(东莞)有限公司_提供
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式