我用我的代码解出来变量一直保持在初始值求大侠指教!

来源：学生作业帮编辑：作业帮分类：综合作业时间：2024/05/24 17:29:11

我用我的代码解出来变量一直保持在初始值求大侠指教!
我想用matlab解微分方程组,为什么我无法解出来?（解出来变量一直保持在初始值）求大侠指教!下面是我的代码.
clear all
clc
global k1 k2 C0 kdEG kaEG kdCBH kaCBH Pf EGt CBHt
k1 = 40;
k2 = 1*10^14;
C0 = 5.56*10^-6;
kdEG = 1*10^-6;
kaEG = 5*10^6;
kdCBH = 1*10^-7;
kaCBH = 5*10^6;
Pf = 0.01;
EGt = input('Enter EGt:');
CBHt = input('Enter CBHt:');
S10 = Pf;
S20 = S10*Pf;
[t,y]=ode45('functionall',[0,15],[0 0 S10 S20 0]);
plot(t,y(3),t,y(4),t,y(5));
%%%%%%%%%%%%%%%%%
function z = functionall(t,y)
global k1 k2 C0 kdEG kaEG kdCBH kaCBH Pf EGt CBHt
z = zeros(5,1);
dydt(1) =((y(3)-y(2)-y(1))*(EGt-y(1)*C0)-kdEG*y(1))*kaEG;
dydt(2) =((y(3)-y(2)-y(1))*(CBHt-y(2)*C0)-kdCBH*y(2))*kaCBH;
dydt(3) = -Pf*(k1*y(1)+k2*y(2)*y(4)*C0);
dydt(4) = k1*y(1)-k2*y(4)*y(4)*y(2)*C0/y(3);
dydt(5) = k2*y(2)*y(4)*C0;
end
有没有大侠可以帮我看看哪儿错了.

function dydt= functionall(t,y)
global k1 k2 C0 kdEG kaEG kdCBH kaCBH Pf EGt CBHt
dydt= zeros(5,1);
dydt(1) =((y(3)-y(2)-y(1))*(EGt-y(1)*C0)-kdEG*y(1))*kaEG;
dydt(2) =((y(3)-y(2)-y(1))*(CBHt-y(2)*C0)-kdCBH*y(2))*kaCBH;
dydt(3) = -Pf*(k1*y(1)+k2*y(2)*y(4)*C0);
dydt(4) = k1*y(1)-k2*y(4)*y(4)*y(2)*C0/y(3);
dydt(5) = k2*y(2)*y(4)*C0;
end
把z换成dydt
再问：这样修改之后反而运行不了了。。是有其他什么错误吗？
再答： [t,y]=ode45(@functionall,[0,15],[0 0 S10 S20 0]);

第二张是答案第三张是我的错解,确实是算不出来,但是不知道我具体错在哪里,求指教PS：是第二问望远镜能录像吗我用望远镜观察鸟类,但能不能把看到的东西录下来啊,求各位大侠指教... 谢谢你一直陪着我,用英文翻译出来,要标准的! 槃字在贴吧打不出来?求代码帮我看看这段lingo代码,我不明白为什么一直出现error code11的错误我正在考该死的大学物理试验这题怎么搞的各位大侠指教一下杭电acm 1005 我觉得很简单,不知道怎么出现stack_overflow,求大神指教.下面是代码. 我英语在初中一直很好,可是在高中我的英语成绩一直保持在90分左右,我该怎么办我的代码的错误在哪里,请以改正一下我的代码, 谁能帮我编段优化设计程序啊,一直目标函数、约束、变量,求满足约束且目标最小时的变量取值,万分感谢车辆识别代码错误?我一直用在泉州交警网里面查询违章记录,但是现在突然不知道为什么查不了,一直显示车辆识别代码错误,车牌号手心有痣有什么含义?我左手手心有两颗痣不很清晰,右手有一颗痣很清晰,还请大侠们赐教……在此先谢过指教的朋友!

我用我的代码解出来变量一直保持在初始值 求大侠指教!

我用我的代码解出来变量一直保持在初始值求大侠指教!