作者neilabc127 (康康)
看板MATLAB
标题[问题]如何使用matlab解laplace 联立?
时间Tue May 7 17:50:57 2013
抱歉 我工数不好 想请问要如何用laplace解未知时间函数的联立解
程式如下 但我不知该怎麽写 才能解出答案
c1=sym('c1(t)')
c2=sym('c2(t)')
c3=sym('c3(t)')
comfd=sym('comfd(t)')
cfd=sym('cfd(t)')
syms k1 k2 k3 k4 k5 k6 k7 k11 k0 t
ans1=laplace(k1*cfd-(k2+k3)*c1-diff(c1,t));
ans2=laplace(k3*c1-k4*c2-diff(c2,t));
ans3=laplace(k5*comfd-k6*c3-k7*c3-diff(c3,t));
ans4=laplace(k0*cfd-k1*cfd+k2*c1-diff(cfd,t));
ans5=laplace(k0*cfd+k6*c3-k11*comfd-k5*comfd-diff(comfd,t));
整理後会得到
laplace(c1(t), t, s)= k1*laplace(cfd(t), t, s)/(s+k2+k3)
laplace(c2(t), t, s)= k3*laplace(c1(t), t, s)/(k4+s)
laplace(c3(t), t, s)= k5*laplace(comfd(t), t, s)/(k6+k7+s)
laplace(cfd(t), t, s)= k2*laplace(c1(t), t, s)/(k1+s-k0)
laplace(comfd(t), t, s)= (k6*laplace(c3(t), t, s) + k0*laplace(cfd(t), t,
s))/(k11+k5+s)
slove(‘laplace(c1(t), t, s)= k1*laplace(cfd(t), t, s)/(s+k2+k3)’,’
laplace(c2(t), t, s)= k3*laplace(c1(t), t, s)/(k4+s)’,’ laplace(c3(t), t,
s)= k5*laplace(comfd(t), t, s)/(k6+k7+s)’,’ laplace(cfd(t), t, s)=
k2*laplace(c1(t), t, s)/(k1+s-k0)’,’ laplace(comfd(t), t, s)=
(k6*laplace(c3(t), t, s) + k0*laplace(cfd(t), t, s))/(k11+k5+s)’)
我想得到解联立後的laplace(c1(t),t,s)和其他几个s domain变数的关系
想问这要怎写?因为他只给我Error: The expression to the left of the equals sign is not a valid target
for an assignment.
感谢 感谢(磕头 Q0Q
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.251.182
※ 编辑: neilabc127 来自: 140.112.251.182 (05/07 17:52)