版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

有累人1

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 744
帖子: 5
在線: 16.9小時(shí)
蟲號(hào): 3385453
注冊(cè): 2014-08-27
性別: GG
專業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人

[求助] MATLA B求解微分方程

這個(gè)程序得不到想要的結(jié)果，運(yùn)行可能導(dǎo)致電腦死機(jī)，求大神幫助，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]

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助二階偏微分方程求解已經(jīng)有13人回復(fù)
這種偏微分方程如何求解�。壳笾更c(diǎn) 已經(jīng)有20人回復(fù)
非線性二階微分方程組求解，matlab 已經(jīng)有12人回復(fù)
偏微分方程Matlab求解已經(jīng)有7人回復(fù)
變系數(shù)偏微分方程matlab求解已經(jīng)有6人回復(fù)
matlab 非線性微分方程求解已經(jīng)有3人回復(fù)
微分方程求解（通解中含hypergeom函數(shù)）已經(jīng)有17人回復(fù)
常微分方程求解公式已經(jīng)有7人回復(fù)
偏微分方程求解！已經(jīng)有4人回復(fù)
一個(gè)二階常微分方程的求解方法？已經(jīng)有5人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
復(fù)系數(shù)微分方程求解已經(jīng)有3人回復(fù)
Matlab求解二階偏微分方程組，希望能給出相應(yīng)的求解代碼已經(jīng)有10人回復(fù)
分?jǐn)?shù)階微分方程的求解（用差分法，matlab編程求解）已經(jīng)有9人回復(fù)
求助matlab求解復(fù)雜微分方程已經(jīng)有7人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)

改變

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

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

xzq220120

金蟲 (初入文壇)

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

好東西，謝謝分享

回復(fù)此樓

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主有累人1 的主題更新