版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6303)
>考研 (1412)
>導師招生 (1277)
>文獻求助 (480)
>蟲友互識 (382)
>基金申請 (234)
>休閑灌水 (220)
>考博 (214)
>論文投稿 (158)
>碩博家園 (153)
>招聘信息布告欄 (110)
>博后之家 (93)
>教師之家 (68)
>綠色求助(高懸賞) (44)
>公派出國 (42)
>留學生活 (35)

返回列表

有累人1

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 744
帖子: 5
在線: 16.9小時
蟲號: 3385453
注冊: 2014-08-27
性別: GG
專業(yè): 機構學與機器人

[求助] 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人回復
294求調劑材料與化工專碩已經有5人回復
一志愿華中科技大學，080502，354分求調劑已經有4人回復
一志愿吉林大學材料學碩321求調劑已經有6人回復
085410人工智能專碩317求調劑（0854都可以）已經有3人回復
330求調劑已經有3人回復
一志愿中海洋材料工程專碩330分求調劑已經有5人回復
304求調劑已經有5人回復
311求調劑已經有4人回復
354求調劑已經有7人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

改變

1樓 2015-07-06 09:26:09

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

xzq220120

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1424.5
紅花: 1
帖子: 25
在線: 26.4小時
蟲號: 2391611
注冊: 2013-03-31
專業(yè): 質譜分析

好東西，謝謝分享

回復此樓

2樓2015-07-07 18:37:40

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主有累人1 的主題更新