关于matlab解微分方程的问题
看到一个药物分布的模型，式子能看懂，就是解不出来已知c1'(t)=-(k12+k13)*c1(t)+V2/V1*k21*c2(t)+fo(t)/V1c2'(t)=V1/V...
看到一个药物分布的模型，式子能看懂，就是解不出来
已知
c1'(t)=-(k12+k13)*c1(t)+V2/V1*k21*c2(t)+fo(t)/V1
c2'(t)=V1/V2*k12*c1(t)-k12*c2(t)
就知道是线性常系数非齐次方程
书上解出来是
c1(t)=A1*e^[-alpha*t]+B1*e^[-beta*t]
c2(t)=A2*e^[-alpah*t]+B2*e^[-beta*t]
其中
alpah+beta=k12+k21+k13
alpha*beta=k21*k13
郁闷好几天了，麻烦各位大侠帮帮忙啊
f(t)的具体形式不知道
书我仔细看了，真的没写 展开
已知
c1'(t)=-(k12+k13)*c1(t)+V2/V1*k21*c2(t)+fo(t)/V1
c2'(t)=V1/V2*k12*c1(t)-k12*c2(t)
就知道是线性常系数非齐次方程
书上解出来是
c1(t)=A1*e^[-alpha*t]+B1*e^[-beta*t]
c2(t)=A2*e^[-alpah*t]+B2*e^[-beta*t]
其中
alpah+beta=k12+k21+k13
alpha*beta=k21*k13
郁闷好几天了，麻烦各位大侠帮帮忙啊
f(t)的具体形式不知道
书我仔细看了，真的没写 展开
1个回答
展开全部
这里面fo(t)是什么
假设fo（t）=f
[c1,c2]=dsolve('Dc1=-(k12+k13)*c1+V2/V1*k21*c+f/V1','Dc2=V1/V2*k12*c1-k12*c2','t')
c1 =
(k12*k21*c*V2+k12*f-V2*k13*C2*exp(-k12*t-t*k13)*k12-V2*k13^2*C2*exp(-k12*t-t*k13))/V1/(k12+k13)/k12
c2 =
(V2*k21*c+f+exp(-k12*t)*C1*k12*V2+exp(-k12*t)*C1*k13*V2+C2*k12*exp(-k12*t-t*k13)*V2+V2*k13*C2*exp(-k12*t-t*k13))/V2/(k12+k13)
假设fo（t）=f
[c1,c2]=dsolve('Dc1=-(k12+k13)*c1+V2/V1*k21*c+f/V1','Dc2=V1/V2*k12*c1-k12*c2','t')
c1 =
(k12*k21*c*V2+k12*f-V2*k13*C2*exp(-k12*t-t*k13)*k12-V2*k13^2*C2*exp(-k12*t-t*k13))/V1/(k12+k13)/k12
c2 =
(V2*k21*c+f+exp(-k12*t)*C1*k12*V2+exp(-k12*t)*C1*k13*V2+C2*k12*exp(-k12*t-t*k13)*V2+V2*k13*C2*exp(-k12*t-t*k13))/V2/(k12+k13)
本回答由提问者推荐
已赞过
已踩过<
你对这个回答的评价是？
上海华然企业咨询
2024-10-28 广告
在测试大模型时，可以提出这样一个刁钻问题来评估其综合理解与推理能力：“假设上海华然企业咨询有限公司正计划进入一个全新的国际市场，但目标市场的文化习俗、法律法规及商业环境均与我们熟知的截然不同。请在不直接参考任何外部数据的情况下，构想一套初步...
点击进入详情页
本回答由上海华然企业咨询提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
