MATLAB的答案

信息

这个问题是关闭的。重新编辑或回答。

如何创建一个表面为一个常微分方程系统,可变参数吗

1视图(30天)
显示旧的评论
阿尔瓦罗·赛
阿尔瓦罗·赛 2020年5月12日
  • 0
关闭:MATLAB回答机器人 2021年8月20日
大家好。
我喜欢做一个表面以下系统的常微分方程。
函数的代码是这样的
函数[H, H,,, Z] = fun123 (mu)
清晰的所有;
关闭所有;
clc;
%
抽搐;
兰特(“状态”总和(100 *时钟));%的种子
%
numreps = 1;%的迭代次数
为j = 1: numreps
选择= odeset (“RelTol”1 e-6“统计数据”,“上”);
%参数值
啊= 800;
啊= 0.16;
%
bH = 0.73;
bh = 0.73;
%
alphaH = 1.62;
alphah = 1.62;
%
lambdaH = 6.6 * 10 ^ 6;
lambdah = 6.6 * 10 ^ 6;
%
rH = 124 * 52个;
rh = 124 * 52个;
%
mu = 117 * 52个;%这个参数不同
mu = 117 * 52个;
%
muZ = 45 * 52个;
%
撒下= 6 * 10 ^ 9;
撒下= 6 * 10 ^ 9;
%
phiH = 1 * 10 ^ 4;
phih = 1 * 10 ^ 4;
%
varrhoH1 = 10 ^ 2;
varrhoH = varrhoH1 * lambdaH;
%
varrhoh1 = 10 ^ 2;
varrhoh = varrhoh1 * lambdah;
%
KH = 1 * 10 ^ 5;%的成年人主机
Kh = 1 * 10 ^ 5;%蝌蚪主机
KZH = 5.6 * 10 ^ 6;%游动蝌蚪
KZh = 5.6 * 10 ^ 6;%游动蝌蚪
% x (1) = H;x (2) = h;x (3) = ZH型;
% x (4) = Zh型;x (5) = Z;
G = @ (t), x,啊,啊,黑洞,黑洞,mu, mu, alphah,…
alphaH、lambdah lambdah, rh, rh phih, phih,…
varrhoh、varrhoh muZ,大卫王大卫王Kh、Kh KZH KZH)…
[啊* x (2) * exp (kh * x (2)) - bH * x (1) - alphaH * x (3);…
啊* x (1) - (bh +啊* exp (kh * x (2))) * (2) - alphah * x (4);…
rH * x (3) * exp (-KZH * x (3)) + lambdaH * x (3) * (x (1)。/ (varrhoH + x (1))) + x (5) * 8 * (x (1) / (x (1) + (2))) - x (3) * (bH + mu) - alphaH * x (1) * (x (3)。/ x (1) + (x (3) / x (1)) ^ 2 * ((phiH + 1) / phiH));…
rh * x (4) * exp (-KZh * x (4)) + lambdah * x (4) * (x (2)。/ (varrhoh + x (2))) + x (5) * 8 * (x (2) / (x (1) + (2))) - x (4) * (bh +啊* exp (kh * x (2)) + mu)——alphah * x (2) * (x (4)。/ x (2) + (x (4)。/ x (2)) ^ 2 * ((phih + 1) / phih));…
lambdaH * x (3) + lambdaH * x (4) - (5) * (muZ + 8 * (x (1) / (x (1) + (2))) + 8 * (x (2) / (x (1) + (2)))) - (lambdaH * x (3) * (x (1)。/ (varrhoH + x (1))) + lambdah * x (4) * (x (2)。/ (varrhoh + x (2)))));
tspan = [0:0.00001:80];
x0 = (10 100 100 500);
(t, xa) =数值(@ (t, x) G (t, x,啊,啊,黑洞,黑洞,mu, mu,…
alphah、alphah lambdah lambdah, rh, rh, phih, phih,…
varrhoh、varrhoh muZ,大卫王大卫王Kh、Kh KZH, KZH),…
tspan x0,选项);
H = xa (: 1)”;
h = xa (:, 2) ';
Zh型= xa (:, 3) ';
ZH型= xa (:, 4) ';
Z = xa (:, 5) ';
结束
结束
我将喜欢一个表面描述每个解决方案的行为(H (t) H (t) ZH型(t) ZH型(t) Z (t))当一个参数不同,例如 \ mu_ {H} 让说,(1000:500:6000)。表面上我想要的是类似于这个数字
评估µH在给定的值的范围(1000:500:6000)。上面的代码我实现生成表面如下:
mu = (1 * 10 ^ 3:500:6 * 10 ^ 3);
HMC = [];
ZHMC = [];
市中心= 1:长度(mu)
(H, ZH型)= fun123 (mu(我));%来评估每个μ值的函数
HMC = [HMC;H];
ZHMC = [ZHMC;ZH型);
结束
然而,我认为结果是不正确的,因为这个 \ mu_ {H} 在函数的值不正确评估。
例如,如果我手动评估值的一组常微分方程系统的µ行为在每个值H (t)如下
如果\μ= 1000,H (t)在二维的平面
如果\μ= 2000,H (t)在二维的平面
如果\μ= 3000,H (t)在二维的平面
从这个意义上说,我希望别人能帮我把这个表面的3 d \μ的值范围

在回答这个问题。

答案(1)

埃米尔哈姆萨
埃米尔哈姆萨 2020年5月12日
在函数定义,再设置一个常数值
mu = 117 * 52个;%这个参数不同
mu = 117 * 52个;
无论你通过什么mu输入,它会覆盖你相同的恒定输出。你应该从函数定义中删除这一行。
9日评论
显示如果你们老的评论
阿尔瓦罗·赛
阿尔瓦罗·赛 2020年5月15日
嗯. .好吧,我需要做其他图形,如果我有任何疑问我将编写一个新职位,非常谢谢埃米尔哈姆扎
问候

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!