为什么不M2函数被调用?有人能解释这个错误是什么以及如何修复它吗?

32视图(30天)

显示旧的评论

cindyawati cindyawati 约20小时前

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1984614-why-doesn-t-m2-function-get-called-could-somebody-please-explain-what-this-error-is-and-how-to-fix

编辑: cindyawati cindyawati约1小时前

我有函数平方米

                          函数M2 = fRK4M2 (M1、M2)
                         
                          M1 = 10;
                         
                          M2 = 0;
                         
                          mu_2 = 10 ^ 3;
                         
                          K1 = 10 ^ 4;
                         
                          K2 = 5 * 10 ^ 4;
                         
                          M2 = (K1 * M1 * M1) - M1 (K2 * *平方米)(mu_2 *平方米);

这是我的M1函数

                          函数M1 = fRK4M1 (M1, M2, M3, O, P)
                         
                          δ= 50;
                         
                          K1 = 10 ^ 4;
                         
                          Ko = 0.1;
                         
                          n = 3;
                         
                          Oa = 10;
                         
                          Pa = 100;
                         
                          mu_1 = 10 ^ 3;
                         
                          K2 = 5 * 10 ^ 4;
                         
                          K3 = 10 ^ 3;
                         
                          γ= 75;
                         
                          M1 =(δ* M1 * (1 - (M1 /γ))2 * K1 * M1 * M1-M1 * (K2。*平方米)- ((Oa-n) * K3 * M1 * M3)——(Ko * M1 (Pa-Oa) * * O) - (mu_1 * M1));

这是我的主要功能

                          clc;明确;
                         
                          %输入时间
                         
                          t (1) = 0;
                         
                          dt = 0.1;%的时间间隔
                         
                          t = 0: dt: 100;%时间跨度
                         
                          %输入空数组
                         
                          T = 0(长度(T), 1);% t空数组
                         
                          M1 = 0(长度(t), 1);%为M1空数组
                         
                          M2 = 0(长度(t), 1);%对M2空数组
                         
                          为j = 1:长度((t))
                         
                          T (j + 1) = T (j) + dt;
                         
                          M1 (j + 1) = M1 (j) + 1. / (1 + exp (- t (j)));
                         
                          k1M2 = dt * fRK4M2 (M2 (j), M1 (j));
                         
                          k2M2 = dt * fRK4M2 (M2 (j) + k1M2/2, M1 (j) + k1M2/2);
                         
                          k3M2 = dt * fRK4M2 (M2 (j) + k2M2/2, M1 (j) + k2M2/2);
                         
                          k4M2 = dt * fRK4M2 (M2 (j) + k3M2, M1 (j) + k3M2/2);
                         
                          平方米(j + 1) =平方米(j) + 1/6 * (k1M2 + 2 * k2M2 + 2 * k3M2 + k4M2);
                         
                          k1M1 = dt * fRK4M1 (M1 (j), M2 (j));
                         
                          k2M1 = dt * fRK4M1 (M2 (j) + k1M2/2, M1 (j) + k1M1/2);
                         
                          k3M1 = dt * fRK4M1 (M2 (j) + k2M2/2, M1 (j) + k2M1/2);
                         
                          k4M1 = dt * fRK4M1 (M2 (j) + k3M2/2, M1 (j) + k3M1/2);
                         
                          M1 (j + 1) = M1 (j) + (1/6 * (k1M1 + (2 * k2M1) + (2 * k3M1) + k4M1));
                         
                          结束

我得到的结果并不如预期的,有人请告诉我我的错和如何解决它。谢谢

2的评论
显示1年长的评论藏1年长的评论

cindyawati cindyawati 大约8小时前

编辑:cindyawati cindyawati 大约8小时前

谢谢你的回应 @Torsten 是的我的代码来解决系统的微分方程。如果我用龙格库塔奥德4日是可能的吗?

登录置评。

在回答这个问题。

答案(1)

Lakshya 约19小时前

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1984614-why-doesn-t-m2-function-get-called-could-somebody-please-explain-what-this-error-is-and-how-to-fix answer_1257864

你好,

您提供的代码初始化M1和M2到特定值(M1 = 10平方米= 0)内 fRK4M2 函数,它可能覆盖的值作为参数传递。如果你打算使用传递的参数,您可以删除这些作业的 fRK4M2 函数。

希望这个有帮助。

6个评论
显示5年长的评论隐藏5年长的评论

cindyawati cindyawati 约1小时前

编辑:cindyawati cindyawati 约1小时前

你好 @Lakshya 我已经输入参数对M3、O和P但我对M2的阴谋,O和P不像预期的(我期望得到像乙状结肠图)。为什么我的图不显示乙状结肠图虽然我已经包括在主函数吗? 这是我的主要功能

                                   clc;明确;
                                  
                                   %输入时间
                                  
                                   t (1) = 0;
                                  
                                   dt = 0.1;%的时间间隔
                                  
                                   t = 0: dt: 100;%时间跨度
                                  
                                   %输入空数组
                                  
                                   T = 0(长度(T), 1);% t空数组
                                  
                                   M1 = 0(长度(t), 1);%为M1空数组
                                  
                                   M2 = 0(长度(t), 1);%对M2空数组
                                  
                                   M3 = 0(长度(t), 1);%对M3空数组
                                  
                                   O = 0(长度(t), 1);O %空数组
                                  
                                   P = 0(长度(t), 1);% P空数组
                                  
                                   %初始条件
                                  
                                   M1 = 10;
                                  
                                   M2 = 0;
                                  
                                   M3 = 0;
                                  
                                   O = 0;
                                  
                                   P = 0;
                                  
                                   为j = 1:长度((t))
                                  
                                   T (j + 1) = T (j) + dt;
                                  
                                   M1 (j + 1) = M1 (j) + 1. / (1 + exp (- t (j)));
                                  
                                   k1M2 = dt * fRK4M2 (M2 (j), M1 (j));
                                  
                                   k2M2 = dt * fRK4M2 (M2 (j) + k1M2/2, M1 (j) + k1M2/2);
                                  
                                   k3M2 = dt * fRK4M2 (M2 (j) + k2M2/2, M1 (j) + k2M2/2);
                                  
                                   k4M2 = dt * fRK4M2 (M2 (j) + k3M2, M1 (j) + k3M2/2);
                                  
                                   平方米(j + 1) =平方米(j) + 1/6 * (k1M2 + 2 * k2M2 + 2 * k3M2 + k4M2);
                                  
                                   k1M3 = dt * fRK4M3 (M1 (j), M3 (j));
                                  
                                   k2M3 = dt * fRK4M3 (M3 (j) + k1M3/2, M1 (j) + k1M3/2);
                                  
                                   k3M3 = dt * fRK4M3 (M3 (j) + k2M3/2, M1 (j) + k2M3/2);
                                  
                                   k4M3 = dt * fRK4M3 (M3 (j) + k3M3, M1 (j) + k3M3);
                                  
                                   M3 (j + 1) = M3 (j) + 1/6 * (k1M3 + 2 * k2M3 + 2 * k3M3 + k4M3);
                                  
                                   k1O = dt * fRK4O (M1 (j), O (j));
                                  
                                   k2O = dt * fRK4O (O (j) + k1O / 2, M1 (j) + k1O / 2);
                                  
                                   k3O = dt * fRK4O (O (j) + k2O / 2, M1 (j) + k2O / 2);
                                  
                                   k4O = dt * fRK4O (O (j) + k3O M1 (j) + k3O);
                                  
                                   O (j + 1) = O (j) + 1/6 * (k1O + 2 * k2O + 2 * k3O + k4O);
                                  
                                   k1P = dt * fRK4P (P (j), M1 (j));
                                  
                                   k2P = dt * fRK4P (P (j) + k1P / 2, M1 (j) + k1P / 2);
                                  
                                   k3P = dt * fRK4P (P (j) + k2P / 2, M1 (j) + k2P / 2);
                                  
                                   k4P = dt * fRK4P (P (j) + k3P, M1 (j) + k3P / 2);
                                  
                                   P (j + 1) = P (j) + 1/6 * (k1P + 2 * k2P + 2 * k3P + k4P);
                                  
                                   k1M1 = dt * fRK4M1 (M1 (j), M2 (j), M3 (j), O (j)、P (j));
                                  
                                   k2M1 = dt * fRK4M1 (M2 (j) + k1M2/2, M3 (j) + k1M3/2, O (j) + k1O / 2, P (j) + k1P / 2, M1 (j) + k1M1/2);
                                  
                                   k3M1 = dt * fRK4M1 (M2 (j) + k2M2/2, M3 (j) + k2M3/2, O (j) + k2O / 2, P (j) + k2P / 2, M1 (j) + k2M1/2);
                                  
                                   k4M1 = dt * fRK4M1 (M2 (j) + k3M2/2, M3 (j) + k3M3/2, O (j) + k3O / 2, P (j) + k3P / 2, M1 (j) + k3M1/2);
                                  
                                   M1 (j + 1) = M1 (j) + (1/6 * (k1M1 + (2 * k2M1) + (2 * k3M1) + k4M1));
                                  
                                   结束
                                  
                                      未识别的函数或变量“fRK4M2”。

                                   图;
                                  
                                   情节(T, M1,“r”,“线宽”3)
                                  
                                   包含(“时间”)
                                  
                                   ylabel (“M1”)
                                  
                                   图
                                  
                                   情节(T, M2,“b”,“线宽”3)
                                  
                                   包含(“时间”)
                                  
                                   ylabel (“平方米”)
                                  
                                   图
                                  
                                   情节(T, M3,‘g’,“线宽”3)
                                  
                                   包含(“时间”)
                                  
                                   ylabel (“M3”)
                                  
                                   图
                                  
                                   情节(T,啊,“b”,“线宽”5)
                                  
                                   包含(“时间”)
                                  
                                   ylabel (“O”)
                                  
                                   图
                                  
                                   情节(T, P,“r”,“线宽”5)
                                  
                                   包含(“时间”)
                                  
                                   ylabel (“P”)

这是我的结果图。为什么我的O和P (P图相同的像阿图)图没有变化吗

cindyawati cindyawati 约1小时前

这是我的函数M2, M3, O和P @Lakshya

                                   函数M2 = fRK4M2 (M1, M2)
                                  
                                   M1 = 10;
                                  
                                   % M2 = 0;
                                  
                                   mu_2 = 10 ^ 3;
                                  
                                   K1 = 10 ^ 4;
                                  
                                   K2 = 5 * 10 ^ 4;
                                  
                                   M2 = (K1 * M1 * M1) - M1 (K2 * *平方米)(mu_2 *平方米);

M3函数

                                   函数M3 = fRK4M3 (M1, M3)
                                  
                                   % M1 = 10;
                                  
                                   M2 = 0;
                                  
                                   % M3 = 0;
                                  
                                   mu_3 = 10 ^ 3;
                                  
                                   K2 = 5 * 10 ^ 4;
                                  
                                   K3 = 10 ^ 3;
                                  
                                   M3 = (K2 * M1 * M2) - (K3 * M1 * M3) - (mu_3 * M3);

O函数

                                   函数O = fRK4O (O, M1)
                                  
                                   % M1 = 10;
                                  
                                   M3 = 0;
                                  
                                   O = 0;
                                  
                                   Ko = 0.1;
                                  
                                   mu_o = 10 ^ 4;
                                  
                                   K3 = 10 ^ 3;
                                  
                                   O = (K3 * M1 * M3) - (Ko * M1 * O) - (mu_o * O);

登录置评。

在回答这个问题。

类别

MATLAB 语言基础知识矩阵和数组创建和连接矩阵

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

社区寻宝

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

开始狩猎!

为什么不M2函数被调用?有人能解释这个错误是什么以及如何修复它吗?

2的评论
显示1年长的评论藏1年长的评论

答案(1)

6个评论
显示5年长的评论隐藏5年长的评论

另请参阅

类别

标签

社区寻宝

为什么不M2函数被调用?有人能解释这个错误是什么以及如何修复它吗?

2的评论 显示1年长的评论藏1年长的评论

答案(1)

6个评论 显示5年长的评论隐藏5年长的评论

另请参阅

类别

标签

社区寻宝

2的评论
显示1年长的评论藏1年长的评论

6个评论
显示5年长的评论隐藏5年长的评论