| 查看: 396 | 回复: 1 | ||
[求助]
MATLA B求解微分方程
|
|
这个程序得不到想要的结果，运行可能导致电脑死机，求大神帮助，TXT在附件中。 clc;clear all; V=textread('D:\Program Files (x86)\MATLAB71\work\1.txt'); i=1; c=1; re=1; U=[V(i,1) V(i,2) V(i,3);V(i,4) V(i,5) V(i,6);V(i,7) V(i,8) V(i,9)]; r=(U+U')/2; w=(U-U')/2; Dr=(0.0184*exp(-0.7148*c*re))*(sqrt(2*sum(diag(r*r)))); A=[1/3 0 0;0 1/3 0;0 0 1/3]; syms a11 a12 a13 a21 a22 a23 a31 a32 a33; f1=eval((w(1,1)*a11+w(1,2)*a21+w(1,3)*a31-w(1,1)*a11-w(2,1)*a12-w(3,1)*a13)/(-2)+(r(1,1)*a11+r(1,2)*a21+r(1,3)*a31+r(1,1)*a11+r(2,1)*a12+r(3,1)*a13-2*r(1,1)*(6*a11/7-3/35)-2*r(1,2)*3*a12/7-2*r(1,3)*3*a13/7-2*r(2,1)*(2*a12+a21)/7-2*r(2,2)*((a11+a22)/7-1/35)-2*r(2,3)*a23/7-2*r(3,1)*(2*a13+a31)/7-2*r(3,2)*a32/7-2*r(3,3)*((a11+a33)/7-1/35))/2+Dr*(1-3*a11)); f2=eval((w(1,1)*a12+w(1,2)*a22+w(1,3)*a32-w(1,2)*a11-w(2,2)*a12-w(3,2)*a13)/(-2)+(r(1,1)*a12+r(1,2)*a22+r(1,3)*a32+r(1,2)*a11+r(2,2)*a12+r(3,2)*a13-2*r(1,1)*(a12+2*a21)/7-2*r(1,2)*((a11+a22)/7-1/35)-2*r(1,3)*a23/7-2*r(2,1)*((a11+a22)/7-1/35)-2*r(2,2)*3*a12/7-2*r(2,3)*a13/7-2*r(3,1)*a23/7-2*r(3,2)*a13/7-2*r(3,3)*a12/7)/2-3*Dr*a12); f3=eval((w(1,1)*a13+w(1,2)*a23+w(1,3)*a33-w(1,3)*a11-w(2,3)*a12-w(3,3)*a13)/(-2)+(r(1,1)*a13+r(1,2)*a23+r(1,3)*a33+r(1,3)*a11+r(2,3)*a12+r(3,3)*a13-2*r(1,1)*(a13+2*a31)/7-2*r(1,2)*a32/7-2*r(1,3)*((a11+a33)/7-1/35)-2*r(2,1)*a32/7-2*r(2,2)*a13/7-2*r(2,3)*a12/7-2*r(3,1)*((a11+ a33)/7-1/35)-2*r(3,2)*a12/7-2*r(3,3)*3*a13/7)/2-3*Dr*a13); f4=eval((w(2,1)*a11+w(2,2)*a21+w(2,3)*a31-w(1,1)*a21-w(2,1)*a22-w(3,1)*a23)/(-2)+(r(2,1)*a11+r(2,2)*a21+r(2,3)*a31+r(1,1)*a21+r(2,1)*a22+r(3,1)*a23-2*r(1,1)*3*a21/7-2*r(1,2)*((a11+a22)/7-1/35)-2*r(1,3)*a23/7-2*r(2,1)*((a11+a22)/7-1/35)-2*r(2,2)*(a21+2*a12)/7-2*r(2,3)*a13/7-2*r(3,1)*a23/7-2*r(3,2)*a13/7-2*r(3,3)*a21/7)/2-3*Dr*a21); f5=eval((w(2,1)*a12+w(2,2)*a22+w(2,3)*a32-w(1,2)*a21-w(2,2)*a22-w(3,2)*a23)/(-2)+(r(2,1)*a12+r(2,2)*a22+r(2,3)*a32+r(1,2)*a21+r(2,2)*a22+r(3,2)*a23-2*r(1,1)*((a11+a22)/7-1/35)-2*r(1,2)*(2*a21+a12)/7-2*r(1,3)*a13/7-2*r(2,1)*3*a21/7-2*r(2,2)*(6*a22/7-3/35)-2*r(2,3)*3*a23/7-2*r(3,1)*a31/7-2*r(3,2)*(2*a23+a32)/7-2*r(3,3)*((a22+a33)/7-1/35))/2+Dr*(1-3*a22)); f6=eval((w(2,1)*a13+w(2,2)*a23+w(2,3)*a33-w(1,3)*a21-w(2,3)*a22-w(3,3)*a23)/(-2)+(r(2,1)*a13+r(2,2)*a23+r(2,3)*a33+r(1,3)*a21+r(2,3)*a22+r(3,3)*a23-2*r(1,1)*a23/7-2*r(1,2)*a31/7-2*r(1,3)*a21/7-2*r(2,1)*a31/7-2*r(2,2)*(a23+2*a32)/7-2*r(2,3)*((a22+a33)/7-1/35)-2*r(3,1)*a21/7-2*r(3,2)*((a22+a33)/7-1/35)-2*r(3,3)*3*a23/7)/2-3*Dr*a22); f7=eval((w(3,1)*a11+w(3,2)*a21+w(3,3)*a31-w(1,1)*a31-w(2,1)*a32-w(3,1)*a33)/(-2)+(r(3,1)*a11+r(3,2)*a21+r(3,3)*a31+r(1,1)*a31+r(2,1)*a32+r(3,1)*a33-2*r(1,1)*3*a31/7-2*r(1,2)*a32/7-2*r(1,3)*((a11+a33)/7-1/35)-2*r(2,1)*a32/7-2*r(2,2)*a31/7-2*r(2,3)*a12/7-2*r(3,1)*((a11+a33)/7-1/35)-2*r(3,2)*a12/7-2*r(3,3)*(a31+2*a13)/7)/2-3*Dr*a31); f8=eval((w(3,1)*a12+w(3,2)*a22+w(3,3)*a32-w(1,2)*a31-w(2,2)*a32-w(3,2)*a33)/(-2)+(r(3,1)*a12+r(3,2)*a22+r(3,3)*a32+r(1,2)*a31+r(2,2)*a32+r(3,2)*a33-2*r(1,1)*a32/7-2*r(1,2)*a31/7-2*r(1,3)*a21/7-2*r(2,1)*a31/7-2*r(2,2)*3*a32/7-2*r(2,3)*((a22+a33)/7-1/35)-2*r(3,1)*a21/7-2*r(3,2)*((a22+a33)/7-1/35)-2*r(3,3)*(a32+2*a23)/7)/2-3*Dr*a32); f9=eval((w(3,1)*a13+w(3,2)*a23+w(3,3)*a33-w(1,3)*a31-w(2,3)*a32-w(3,3)*a33)/(-2)+(r(3,1)*a13+r(3,2)*a23+r(3,3)*a33+r(1,3)*a31+r(2,3)*a32+r(3,3)*a33-2*r(1,1)*((a11+a33)/7-1/35)-2*r(1,2)*a12/7-2*r(1,3)*(2*a31+a13)/7-2*r(2,1)*a21/7-2*r(2,2)*((a22+a33)/7-1/35)-2*r(2,3)*(2*a32+a23)/7-2*r(3,1)*3*a31/7-2*r(3,2)*3*a32/7-2*r(3,3)*(6*a33/7-3/35))/2+Dr*(1-3*a33)); [a11,a12,a13] = dsolve(strcat('Da11 =',char(f1)),strcat('Da12 =',char(f2)),strcat('Da13 =',char(f3)),'a11(0) = num2str(A(1,1))','a12(0) = num2str(A(1,2))','a13(0) = num2str(A(1,3))','t') [a21,a22,a23] = dsolve(strcat('Da21 =',char(f4)),strcat('Da22 =',char(f5)),strcat('Da23 =',char(f6)),'a21(0) = num2str(A(2,1))','a22(0)= num2str(A(2,2))','a23(0) = num2str(A(2,3))','t') [a31,a32,a33] = dsolve(strcat('Da31 =',char(f7)),strcat('Da32 =',char(f8)),strcat('Da33 =',char(f9)),'a31(0)=num2str(A(3,1))','a32(0)=num2str(A(3,2))','a33(0)=num2str(A(3,3))','t') [a11,a12,a13,a21,a22,a23,a31,a32,a33] = dsolve(strcat('Da11 =',char(f1)),strcat('Da12 =',char(f2)),strcat('Da13 =',char(f3)),strcat('Da21 =',char(f4)),strcat('Da22 =',char(f5)),strcat('Da23 =',char(f6)),strcat('Da31 =',char(f7)),strcat('Da32 =',char(f8)),strcat('Da33 =',char(f9)),'a11(0)=num2str(A(1,1))','a12(0)=num2str(A(1,2))','a13(0)=num2str(A(1,3))','a21(0)=num2str(A(2,1))','a22(0)=num2str(A(2,2))','a23(0)=num2str(A(2,3))','a31(0)=num2str(A(3,1))','a32(0)=num2str(A(3,2))','a33(0)=num2str(A(3,3))','t') t=0.01325; A=[a11 a12 a13;a21 a22 a23;a31 a32 a33] |
回复此楼» 猜你喜欢
论文撤稿了
已经有8人回复
-
博士申请
已经有7人回复
-
青B发送上会通知了吗
已经有8人回复
-
化学专业申博
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审有感
已经有27人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助二阶偏微分方程求解
已经有13人回复
- 这种偏微分方程如何求解啊？求指点
已经有20人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- 偏微分方程Matlab求解
已经有7人回复
- 变系数偏微分方程matlab求解
已经有6人回复
- matlab 非线性微分方程求解
已经有3人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- 常微分方程求解公式
已经有7人回复
- 偏微分方程求解！
已经有4人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 复系数微分方程求解
已经有3人回复
- Matlab求解二阶偏微分方程组，希望能给出相应的求解代码
已经有10人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 求助matlab求解复杂微分方程
已经有7人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
2楼2015-07-07 18:37:40
10