Znn3bq.jpeg
²é¿´: 993  |  »Ø¸´: 3

zhaoxyah

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

[ÇóÖú] ÇëÎÊmatlab½â·½³Ì×é½á¹ûÈçºÎ±£´æ ÒÑÓÐ1È˲ÎÓë

matlab³ÌÐòÈçÏ£¬Ïë°ÑÔËÐнá¹û±£´æÖÁÒ».txtÎļþ£¬±ãÓÚÓÃORIGIN»­Í¼·ÖÎö£¬ÇëÎʸÃÈçºÎÉèÖÃÈçÏÂÔËËã½á¹û±£´æ£¬ÆÚ´ý´óÏÀÖ¸µ¼
function solve_noneqs2
clear all;clc
format long
global thita
x0=[0.8 -5 0.1];
for i=1:1:88
    thita=i*3.14/180.0;
[x]=fsolve(@objfun,x0)
end
function f=objfun(x)
global thita
gam=0.2;
gama=gam/(cos(thita)+gam*(1-cos(thita)));  
nn=0;

mie=1836.0*4.0;
miu=1.0/mie;
tie=0.01;
rou=0.2;
uiy0=-0.01;
e0=-0.01;
f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));
f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));
f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...
+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
½ö¹©²Î¿¼£¬ÔËÐÐÍê±Ïºó£¬Çëµ½DÅÌÕÒÃûΪ shuju.xls  µÄexcelÎĵµ
CODE:
function solve_noneqs2
clear all;clc
format long
global thita
x0=[0.8 -5 0.1];
for i=1:1:88
    thita=i*3.14/180.0;
x(1:3,i)=fsolve(@objfun,x0);

end
xlswrite('D:\shuju.xls',x)


function f=objfun(x)
global thita
gam=0.2;
gama=gam/(cos(thita)+gam*(1-cos(thita)));  
nn=0;

mie=1836.0*4.0;
miu=1.0/mie;
tie=0.01;
rou=0.2;
uiy0=-0.01;
e0=-0.01;
f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));
f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));
f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...
+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);

ÇëÎÊmatlab½â·½³Ì×é½á¹ûÈçºÎ±£´æ
¸½Í¼1.png

MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2015-01-28 14:23:11
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
zhaoxyah: ½ð±Ò+30, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸ 2015-01-28 14:34:46
»òÕ߸ü¼òµ¥µã£¬°Ñx½á¹ûÏÔʾÔÚÖ÷³ÌÐò´°¿ÚÖУ¬¸´ÖÆÕ³Ìùµ½Origin¼´¿É£º
CODE:
function solve_noneqs2
clear all;clc
format long
global thita
x0=[0.8 -5 0.1];
for i=1:1:88
    thita=i*3.14/180.0;
x(i,1:3)=fsolve(@objfun,x0);

end
x



function f=objfun(x)
global thita
gam=0.2;
gama=gam/(cos(thita)+gam*(1-cos(thita)));  
nn=0;

mie=1836.0*4.0;
miu=1.0/mie;
tie=0.01;
rou=0.2;
uiy0=-0.01;
e0=-0.01;
f(1)=(1+x(1))*x(3)-(1-gama)*sqrt(2*mie*pi)*exp(x(2))/(1+erf(sqrt(-x(2))));
f(2)=gama*x(3)/(1-gama)-x(1)*sqrt(nn*nn-2*mie*x(2))/(1+x(1));
f(3)=(1+erf(sqrt(-x(2)))+2*exp(x(2)))/(1+erf(sqrt(-x(2))))...
+(-x(1)/(miu*nn*nn-2*x(2)))-(1-rou*uiy0*sin(thita)/e0)/(x(3)*x(3)-tie);

Ö÷³ÌÐò´°¿ÚÏÔʾÊý¾Ý£º
x =

   0.000975910082976  -4.456299390614189   0.997581989816602
   0.000974491887485  -4.457933083051565   0.995861494999457
   0.000973373594347  -4.459510833495430   0.994138700333996
   0.000972555181984  -4.461031939273196   0.992414123062480
   0.000972037100394  -4.462495608478736   0.990688283093975
   0.000971820034219  -4.463901127867802   0.988961492335664
   0.000971905879388  -4.465247179177315   0.987234706347581
   0.000972296276528  -4.466532872934121   0.985508231684999
   0.000972993613533  -4.467757048737792   0.983782597604362
   0.000974000785205  -4.468918449757659   0.982058335328984
   0.000975321206273  -4.470015720425237   0.980335977889230
   0.000976958826632  -4.471047403729782   0.978616059975027
   0.000978918148859  -4.472011938243347   0.976899117752067
   0.000981204248167  -4.472907654760871   0.975185688683845
   0.000983822795018  -4.473732772516348   0.973476311357398
   0.000986780080491  -4.474485395033108   0.971771525292660
   0.000990083044619  -4.475163505485457   0.970071870734375
   0.000993739308027  -4.475764961634439   0.968377888460546
   0.000997757207001  -4.476287490185926   0.966690119559326
   0.001002145832424  -4.476728680662956   0.965009105219735
   0.001006915072800  -4.477085978604960   0.963335386490364
   0.001012075661864  -4.477356678171161   0.961669504064674
   0.001017639231092  -4.477537914015268   0.960011998006846
   0.001023618367741  -4.477626652358884   0.958363407518746
   0.001030026678853  -4.477619681242085   0.956724270669641
   0.001036878861917  -4.477513599847341   0.955095124115890
   0.001044190782928  -4.477304806759326   0.953476502823156
   0.001051979562526  -4.476989487155612   0.951868939769982
   0.001060263671332  -4.476563598704270   0.950272965626457
   0.001069063035314  -4.476022856124226   0.948689108451285
   0.001078399152550  -4.475362714215790   0.947117893336593
   0.001088295222606  -4.474578349190693   0.945559842071138
   0.001098776290204  -4.473664638164040   0.944015472763928
   0.001109869404846  -4.472616136541590   0.942485299463856
   0.001121603798543  -4.471427053083234   0.940969831750880
   0.001134011083889  -4.470091222386274   0.939469574311788
   0.001147125475312  -4.468602074436102   0.937985026491452
   0.001160984036531  -4.466952600924484   0.936516681815986
   0.001175626957915  -4.465135317885419   0.935065027485895
   0.001191097867948  -4.463142224192278   0.933630543839829
   0.001207444183753  -4.460964755409359   0.932213703782042
   0.001224717506395  -4.458593732330876   0.930814972167297
   0.001242974067723  -4.456019303561788   0.929434805130577
   0.001262275236763  -4.453230881257064   0.928073649382751
   0.001282688094944  -4.450217069102645   0.926731941425302
   0.001304286091341  -4.446965581400276   0.925410106706889
   0.001327149791125  -4.443463151990092   0.924108558694452
   0.001351367732971  -4.439695431438106   0.922827697847479
   0.001377037414339  -4.435646870739886   0.921567910492499
   0.001404266427297  -4.431300589416395   0.920329567552576
   0.001433173772413  -4.426638225469351   0.919113023142749
   0.001463891383994  -4.421639764261887   0.917918612983278
   0.001496565907280  -4.416283342727678   0.916746652588892
   0.001531360777478  -4.410545024707535   0.915597435219628
   0.001568458661971  -4.404398542185730   0.914471229529014
   0.001608064341769  -4.397814996324893   0.913368276841582
   0.001650408126980  -4.390762510675060   0.912288787997128
   0.001695749925212  -4.383205827368949   0.911232939652124
   0.001744384112993  -4.375105834947984   0.910200869936291
   0.001796645400960  -4.366419013804112   0.909192673290152
   0.001852915937186  -4.357096781816254   0.908208394309342
   0.001913633963974  -4.347084718326792   0.907248020339510
   0.001979304438591  -4.336321638956804   0.906311472489330
   0.002050512157033  -4.324738486258682   0.905398594672314
   0.002127938095470  -4.312256991442435   0.904509140081428
   0.002212379926894  -4.298788049302101   0.903642754405397
   0.002304778009748  -4.284229730930023   0.902798954763590
   0.002406248626099  -4.268464834869572   0.901977103030645
   0.002518126937618  -4.251357844521312   0.901176371702644
   0.002642023134824  -4.232751113689719   0.900395699709655
   0.002779896748912  -4.212460037414608   0.899633734516190
   0.002934156350967  -4.190266871974685   0.898888755266820
   0.003107795337292  -4.165912732114099   0.898158569228656
   0.003304579969084  -4.139087090964472   0.897440370064432
   0.003529314656287  -4.109413800794020   0.896730540332491
   0.003788224095930  -4.076432174511963   0.896024370817481
   0.004089516851567  -4.039570905127388   0.895315652616919
   0.004444239122095  -3.998111347125019   0.894596069133892
   0.004867608652992  -3.951134554796789   0.893854263105887
   0.005381174970585  -3.897442714627276   0.893074356036064
   0.006016468959968  -3.835438671220115   0.892233504032521
   0.006821488777965  -3.762933744132360   0.891297668693953
   0.007872960954398  -3.676826032918474   0.890213869991026
   0.009301380227083  -3.572528644488063   0.888894953575232
   0.011347505935047  -3.442872097900130   0.887186803547526
   0.014508161408529  -3.275768993328194   0.884788750793815
   0.019992655528039  -3.048465485954458   0.881024739323846
   0.031669557869812  -2.709785321526165   0.873989668808664
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
3Â¥2015-01-28 14:28:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhaoxyah

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
3Â¥: Originally posted by ÔÂÖ»À¶ at 2015-01-28 14:28:26
»òÕ߸ü¼òµ¥µã£¬°Ñx½á¹ûÏÔʾÔÚÖ÷³ÌÐò´°¿ÚÖУ¬¸´ÖÆÕ³Ìùµ½Origin¼´¿É£º

function solve_noneqs2
clear all;clc
format long
global thita
x0=;
for i=1:1:88
    thita=i*3.14/180.0;
x(i,1:3)=fsolve(@objf ...

·Ç³£¸ÐлÄúµÄÿһ´ÎÖ¸µ¼£¡×£Ã¿Ì쿪ÐÄ
4Â¥2015-01-28 14:34:19
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zhaoxyah µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸²ÄÁÏ¿ÆÑ§Ó빤³Ì985£¬365·Ö£¬ +7 ²Ä»¯Àî¿É 2026-04-11 9/450 2026-04-12 01:09 by Çï¶¹²ËÑ¿
[¿¼ÑÐ] 280Çóµ÷¼Á +13 wzzzÍõ 2026-04-09 13/650 2026-04-12 00:31 by ÓÂÅʸ߷å0126
[¿¼ÑÐ] ²ÄÁÏ¿¼Ñе÷¼Á +24 ÔÆÄ¾´ï´ï 2026-04-11 26/1300 2026-04-12 00:23 by À¶ÔÆË¼Óê
[˶²©¼ÒÔ°] ÓÐûÓÐѧУ²ÄÁÏרҵÊÕ¿çµ÷(Ò»Ö¾Ô¸085410) +6 momo(Éϰ¶°æ) 2026-04-06 9/450 2026-04-11 22:38 by wj165256
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +3 ban°àСÆß 2026-04-11 3/150 2026-04-11 20:48 by may_ÐÂÓî
[¿¼ÑÐ] 270Çóµ÷¼Á +14 ÑîÀÖ369 2026-04-11 14/700 2026-04-11 20:16 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 085400 328·Ö Çóµ÷¼Á +10 ιÄãÒ»¸ö´ó³È×Ó 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[¿¼ÑÐ] µ÷¼Á +5 ÎĵÀÐÇ̨ 2026-04-11 5/250 2026-04-11 15:01 by ¿­¿­Òª±ä˧
[¿¼ÑÐ] 085410 273·Öµ÷¼Á +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[¿¼ÑÐ] 295·ÖÇóµ÷¼Á +9 ?ÒªÉϰ¶? 2026-04-10 9/450 2026-04-11 10:48 by qingpingzhu
[¿¼ÑÐ] 0854Çóµ÷¼Á +7 assdll 2026-04-05 7/350 2026-04-11 10:34 by Delta2012
[¿¼ÑÐ] 22408 352·ÖÇóµ÷¼Á0854Àà +4 ŬÁ¦µÄÏÄÄ© 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[¿¼ÑÐ] 291 Çóµ÷¼Á +29 »¯¹¤2026½ì±ÏÒµÉ 2026-04-09 29/1450 2026-04-10 22:55 by dick_runner
[¿¼ÑÐ] 285Çóµ÷¼Á +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[¿¼ÑÐ] µ÷¼Á +24 ²»·ê´º 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[¿¼ÑÐ] Çóµ÷¼Á +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[¿¼ÑÐ] 327Çóµ÷¼Á +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[¿¼ÑÐ] 338Çóµ÷¼Á +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 331Çóµ÷¼Á +5 ÕÅÔªÒ» 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[¿¼ÑÐ] 277Çóµ÷¼Á +5 ¿¼Ñе÷¼Álxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û