关于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-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大（共振点）；正弦振动在任一瞬间只包含一种频率的振动，而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率，所以样品上的共振点会同时激发...
点击进入详情页
本回答由富港检测技术（东莞）有限公司_提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
