作者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)