解决Eqn不同变量(Ms)——MATLAB的答案——MATLAB中央 - 金宝app,下载188bet金宝搏,金宝搏官方网站

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

1
链接

直接链接到这个答案

https://it.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;%使用以前的聚合值下初始猜测

米(j,我)= fzero (@ (M) fn (M, a1a4 (j), p4p1(我)),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 - (gamma1-1)) /(γ+ 1);

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

t3 = 2 * gamma4 / (gamma4 - 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,我)= fzero (@ (M) fn (M, a1a4 (j), p4p1(我)),M0);
                             
                              % m = m (j,我);
                             
                              结束
                             
                              结束

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

i循环运行通过的值p4p1一次(1,2,3,…,10000)。

对于这些,为了得到fzero找到M的值,使得函数fn等于零,我们需要给它一个初始猜测。有时你可以逃脱一个初始猜测循环中的每一步。然而,对于你的方程,最好使用一个猜测,来自以前聚合值。即初始猜测一步我将从张一步M的聚合值。1.01的第一个猜测是任意的,但是,幸运的是,它成功了!

fzero调用函数fn除了M0和两个参数。它通过a1a4相等价值的,第i p4p1的价值。fzero然后它的计算来确定,在最佳状态,M, fn的价值的价值尽可能接近于零。

M的返回值被存储在一个4 x10000数组(a1a4代表不同的值的行和列p4p1)的不同的值。

我希望这可以帮助。

Benneth佩雷斯 2021年9月27日

谢谢你我有点明白我需要更多的练习,我看过了你的代码好几次了。

如果我想找到p4p1价值假设M = 3。换句话说向后的工作。我试着调用函数fn得到一个答案,但我遇到了一些麻烦,因为我们最初p4p1时解决所有的M值从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但我想要一个确切的解决方案p4p1当M = 3。

你能斯蒂尔我正确的方向吗?

艾伦·史蒂文斯 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;%使用以前的聚合值下初始猜测

米(j,我)= fzero (@ (M) fn (M, a1a4 (j), p4p1(我)),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)

Z = fn2 (M, a1a4)——p4p1;

结束

函数a1a4 pratio = fn2(米)

γ= 1.667;

gamma4 = 1.4;

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

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

t3 = 2 * gamma4 / (gamma4 - 1);

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

结束

登录置评。

解决Eqn不同变量(女士)

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

接受的答案

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

更多的答案(0)

另请参阅

类别

标签

社区寻宝

解决Eqn不同变量(女士)

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

接受的答案

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

更多的答案(0)

另请参阅

类别

标签

社区寻宝

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

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