MATLAB的答案

一个方程的根当两个参数改变

24日视图(30天)
Yakoob说道
Yakoob说道 在2023年6月19日15:28
评论道: Yakoob说道2023年6月22日11:13
我有这个方程
r0 = 0.05;
k1 = 0.5;
k2 = 0.5;
μ= 0.5;
ρ= 0.5;
ε= 0.25;
K = 1;
α= 0.1;
q = 0.1;
b = 0.8;
ζ= 0.075;
ω= 0.001
信谊Mσ埃塔Msol积极
螺母(M) =μ+(ρ* M / (1 + M));
gro (M) = r0 * (1 + k1 *螺母(M) * (1-k2 *螺母(M)));
林(M) = 1 /(1 +螺母(M));
P (M) =ε* M / (1 + M);
eqn = (q / b)。* gro (M)。*(η+ P (M)) *(1 +螺母(M)) *(1 -(((ασ*。*。*(η+ P (M)) *(1 +螺母(M)) + b。* (q +α)。*(ζ。* M-omega0)))。/σα(b。*。*。* K)))——(q /σ)。*(ζ。* M-omega0) = = 0;
我可以计算这个方程的根当σ多变和η= 0.05使用命令
(num,窝)= numden (lh (eqn));
sigma_num = linspace (0.01、1、20);
为u = 1:元素个数(sigma_num)
Msol (u) = vpa(解决(潜艇(σ,num sigma_num (u))));
结束
会产生20根,每一个与每个值σ。现在,我想计算方程的根当两个参数不同,σ和eta这样eta_num = linspace (0.01、1、20)。我认为我应该有400根,但我能怎么做呢?任何帮助都是赞赏!很多谢谢!

在回答这个问题。

接受的答案

Torsten
Torsten 在2023年6月19日17:02
编辑:Torsten 在2023年6月19日十七13
跑:
“num”是一个多项式的学位7 m .因此,有七根。但是你的代码只给其中一个七。所以我希望你得到一根七你之后。
r0 = 0.05;
k1 = 0.5;
k2 = 0.5;
μ= 0.5;
ρ= 0.5;
ε= 0.25;
K = 1;
α= 0.1;
q = 0.1;
b = 0.8;
ζ= 0.075;
ω= 0.001;
信谊Mσ埃塔Msol积极
螺母(M) =μ+(ρ* M / (1 + M));
gro (M) = r0 * (1 + k1 *螺母(M) * (1-k2 *螺母(M)));
林(M) = 1 /(1 +螺母(M));
P (M) =ε* M / (1 + M);
eqn = (q / b)。* gro (M)。*(η+ P (M)) *(1 +螺母(M)) *(1 -(((ασ*。*。*(η+ P (M)) *(1 +螺母(M)) + b。* (q +α)。*(ζ。* M-omega0)))。/σα(b。*。*。* K)))——(q /σ)。*(ζ。* M-omega0) = = 0;
(num,窝)= numden (lh (eqn));
sigma_num = linspace (0.01、1、20);
eta_num = linspace (0.01、1、20);
n =元素个数(eta_num);
m =元素个数(sigma_num);
M_num = 0 (n, m);
% M_num =细胞(n, m);
为i = 1: n
为j = 1: m
M_num (i, j) = vpa(解决(潜艇(num(ησ),[eta_num(我),sigma_num (j)]) = = 0, M));
% M_num {i, j} =根(sym2poly(潜艇(num(ησ),[eta_num(我),sigma_num (j)))));
结束
结束
M_num
M_num = 20×20
0.0135 0.0146 0.0157 0.0168 0.0180 0.0193 0.0206 0.0219 0.0233 0.0248 0.0263 0.0279 0.0296 0.0314 0.0333 0.0352 0.0373 0.0395 0.0418 0.0442 0.0143 0.0194 0.0247 0.0302 0.0360 0.0420 0.0482 0.0547 0.0615 0.0686 0.0760 0.0837 0.0918 0.1002 0.1090 0.1181 0.1276 0.1376 0.1479 0.1587 0.0150 0.0240 0.0333 0.0431 0.0532 0.0638 0.0748 0.0862 0.0982 0.1106 0.1235 0.1369 0.1509 0.1654 0.1804 0.1959 0.2120 0.2286 0.2458 0.2635 0.0157 0.0284 0.0416 0.0554 0.0697 0.0847 0.1003 0.1164 0.1332 0.1507 0.1687 0.1874 0.2068 0.2268 0.2474 0.2687 0.2905 0.3130 0.3360 0.3597 0.0164 0.0326 0.0495 0.0672 0.0856 0.1047 0.1246 0.1452 0.1666 0.1887 0.2116 0.2352 0.2595 0.2845 0.3102 0.3366 0.3636 0.3912 0.4194 0.4481 0.0170 0.0366 0.0571 0.0785 0.1007 0.1238 0.1478 0.1727 0.1983 0.2249 0.2522 0.2803 0.3091 0.3387 0.3690 0.4000 0.4315 0.4637 0.4965 0.5297 0.0176 0.0405 0.0644 0.0893 0.1152 0.1421 0.1699 0.1987 0.2284 0.2590 0.2905 0.3227 0.3558 0.3895 0.4240 0.4591 0.4948 0.5311 0.5679 0.6053 0.0182 0.0442 0.0713 0.0996 0.1289 0.1594 0.1909 0.2234 0.2569 0.2913 0.3266 0.3627 0.3995 0.4371 0.4754 0.5142 0.5537 0.5937 0.6342 0.6752 0.0187 0.0477 0.0779 0.1094 0.1420 0.1759 0.2108 0.2468 0.2838 0.3217 0.3605 0.4001 0.4405 0.4816 0.5233 0.5656 0.6085 0.6519 0.6958 0.7401 0.0193 0.0510 0.0842 0.1187 0.1545 0.1915 0.2296 0.2689 0.3091 0.3503 0.3924 0.4353 0.4789 0.5232 0.5681 0.6135 0.6595 0.7060 0.7529 0.8002

更多的答案(0)

类别

找到更多的在创建和连接矩阵在帮助中心和文件交换

标签

社区寻宝

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

开始狩猎!