用变变量求解方程(Ms) - MATLAB答案- MATLAB中心 - 金宝app,下载188bet金宝搏,金宝搏官方网站

艾伦·史蒂文斯 2021年9月27日

1
链接

直接链接到这个答案

https://fr.mathworks.com/matlabcentral/answers/1461289-solving-eqn-with-varying-variable-ms#answer_795959

是这样的:

a1a4 = 1。/[1, 2, 4, 10];

n = 10000;

p4p1 = 1: n;

M = 0(元素个数(a1a4), n);

为j = 1:元素个数(a1a4)

m = 1.01;

为i = p4p1

M0 = m;%使用以前的收敛值作为下一次初始估计

M(j,i) = f0 (@(M) fn(M,a1a4(j)，p4p1(i))，M0);

m = m (j,我);

结束

情节(p4p1, M),网格

包含(“p4p1”), ylabel (“米”）

传奇(“a4 / a1 = 1”，“a4 / a1 = 2”，“a4 / a1 = 4”，“a4 / a1 = 10”）；

函数Z = fn (M, a1a4, p4p1)

γ= 1.667;

gamma4 = 1.4;

t1 = (2* M^2-(1-1))/(1 + 1);

t2 = a1a4 * (gamma4-1) / (gamma1-1) * (M - 1 / M);

T3 = -2* (- 1) /(- 1);

Z = t1*(1 - t2)^t3 - p4p1;

结束

4评论
显示隐藏 3年长的评论

Benneth佩雷斯 2021年9月27日

你好艾伦,

这是伟大的!这正是我要找的。就我所知，如果这要求不过分的话。你能给我解释一下下面的部分吗?我真的想在循环中做得更好，我通常会重复复制和粘贴我想做的事情，以得到相同的响应。

我理解J只是一个计数器，基本上会运行4次。

我猜m=1.01是第一个假设，所以方程从m开始。我最难理解的部分是注释部分，以供参考:

谢谢你！

                              为j = 1:元素个数(a1a4)
                             
                              m = 1.01;
                             
                              为i = p4p1
                             
                              % M0 = m;%使用以前的收敛值作为下一次初始估计
                             
                              % M(j,i) = f0 (@(M) fn(M,a1a4(j)，p4p1(i))，M0);
                             
                              % m = m (j,i);
                             
                              结束
                             
                              结束

艾伦·史蒂文斯 2021年9月27日

i循环每次一个遍历p4p1的值(即1,2,3，…, 10000)。

对于每一个，为了让f0找到一个M的值使函数fn等于0，我们需要给它一个初始的猜测。有时，您可以对循环中的每一步进行一个初始猜测。然而，对于你的方程，最好使用来自先前的收敛值的猜测。即第i步的初始猜想设为第i-1步M的收敛值。最初对1.01的猜测是随意的，但幸运的是，它成功了!

除了M0, fzero还使用两个参数调用fn函数。它传递a1a4的第j个值和p4p1的第i个值。然后fzero进行计算，尽可能地确定M的值，从而使fn的值尽可能接近于零。

M的返回值然后存储在4x10000数组中(行代表a1a4的不同值，列代表p4p1的不同值)。

我希望这对你有所帮助。

Benneth佩雷斯 2021年9月27日

谢谢你，我知道我需要更多的练习，我已经看过你的代码好几次了。

如果我想求出p4p1的值，假设M=3。换句话说，就是逆向工作。我尝试调用函数fn来得到答案，但我遇到了一些麻烦，因为我们最初是在解所有的M值，当p4p1从1-10000。

我回到我最初的简单的解决代码，但它给了我一个疯狂的结果

sympref ( “FloatingPointOutput” ,真正的)

γ= 1.667;

gamma4 = 1.667;

M = 3;

信谊 p4p1

a1a4 = 1。/[1, 2, 4, 10];

eqn =(((γ2。* . * M ^ 2-gamma1 + 1) /(γ+ 1))。* (1 - (a1a4)。* ((gamma4-1)。/(γ+ 1))。* ((M . ^ 21)。/ M))。^(2。* gamma4. / (gamma4-1))) -p4p1 = = 0;

p4p1 =解决(eqn p4p1)

从图中我可以看到p4p1应该在2600左右，但我想知道M=3时p4p1的精确解。

你能告诉我正确的方向吗?

艾伦·史蒂文斯 2021年9月28日

注意，p4p1将取决于你为a4a1和m选择的值。你可以按如下方式重新安排程序:

a1a4 = 1。/[1, 2, 4, 10];

n = 10000;

p4p1 = 1:10: n;

M = 0(元素个数(a1a4),元素个数(p4p1));

为j = 1:元素个数(a1a4)

m = 1.01;

为i = 1:元素个数(p4p1)

M0 = m;%使用以前的收敛值作为下一次初始估计

M(j,i) = f0 (@(M) fn(M,a1a4(j)，p4p1(i))，M0);

m = m (j,我);

结束

情节(p4p1, M),网格

包含(“p4p1”), ylabel (“米”）

传奇(“a4 / a1 = 1”，“a4 / a1 = 2”，“a4 / a1 = 4”，“a4 / a1 = 10”）；

%调用函数fn2与所需的值M和a1a4

M_desired = 3;a4a1_desired = 3;a1a4_desired = 1 / a4a1_desired;

disp (['M = 3, a4a1 = 'num2str (a4a1_desired)'， p4p1 = 'a1a4_desired num2str (fn2 (3))))

M = 3, a4a1 = 3, p4p1 = 2273.2071

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

函数Z = fn (M, a1a4, p4p1)

stickline (c = ref (c, 1) - ref (c, 1);

结束

函数a1a4 pratio = fn2(米)

γ= 1.667;

gamma4 = 1.4;

t1 = (2* M^2-(1-1))/(1 + 1);

t2 = a1a4 * (gamma4-1) / (gamma1-1) * (M - 1 / M);

T3 = -2* (- 1) /(- 1);

ratio = t1*(1 - t2)^t3;

结束

登录评论。

变变量(Ms)求解Eqn

0评论
显示隐藏1以上的评论

接受的答案

4评论
显示隐藏 3年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

变变量(Ms)求解Eqn

0评论 显示隐藏1以上的评论

接受的答案

4评论 显示隐藏 3年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

0评论
显示隐藏1以上的评论

4评论
显示隐藏 3年长的评论