±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1116  |  »Ø¸´: 4

zhaoxyah

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

[ÇóÖú] matlab½â·½³Ì½á¹ûÊä³öÎÊÌâ ÒÑÓÐ1È˲ÎÓë

matlab½â·½³ÌÈçÏ£¬Èý¸öδ֪Êý£¬ÏÖÔÚÔÚmatlabÃüÁî´°¿ÚÔËÐнá¹ûÏëÖ»Êä³ö½á¹ûx(3)£¬¸ÃÔõô°ì£¿ÇóÖú£¬Ð»Ð»
function solve_noneqs2
clear all;clc
format long
global einc
x0=[0.8 -2 0.9];
for i=1:1:50
  einc=i;  
x(i,1:3)=fsolve(@objfun,x0)

end

x
function f=objfun(x)



global einc
deltam=1.4;
em=700.0;
thita=10*3.14/180.0;  
gam=deltam*einc*exp(-2.0*(sqrt(einc/em)-1.0))/em;

gama=gam*cos(thita)/(1-gam+gam*cos(thita));  

mie=1836.0*1.0;
miu=1.0/mie;
tie=1.0;
nn=5;
uiy0=-0.01;
e0=-0.01;
beta=0.05;
b=1.0;
rou=0.02302*b;
alfa=1.0;
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+x(1))*(1-sqrt(1+x(1))*sin(thita)*sin(thita)+sqrt(1+x(1))*beta*rou*x(3)/e0)/((1+...
+erf(sqrt(-x(2)))+exp(x(2))/sqrt(-pi*x(2)))/(1+erf(sqrt(-x(2))))...
+(-x(1)/(miu*nn*nn-2*x(2))))+alfa*tie-x(3)*x(3);
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

CelestialCYJ

ľ³æ (СÓÐÃûÆø)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
zhaoxyah: ½ð±Ò+30 2015-02-10 22:25:22
Êä³öxµÄµÚÈýÁÐÂð£¿x(:,3)
2Â¥2015-02-10 18:39:10
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhaoxyah

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

ÒýÓûØÌû:
2Â¥: Originally posted by CelestialCYJ at 2015-02-10 18:39:10
Êä³öxµÄµÚÈýÁÐÂð£¿x(:,3)

Êǽ«x(i,1:3)=fsolve(@objfun,x0)¸Ä³Éx(:£¬3)=fsolve(@objfun,x0)ô£¬ÕâÑùºÃÏñ²»Ðа¡
3Â¥2015-02-10 21:21:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

CelestialCYJ

ľ³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
3Â¥: Originally posted by zhaoxyah at 2015-02-10 21:21:51
Êǽ«x(i,1:3)=fsolve(@objfun,x0)¸Ä³Éx(:£¬3)=fsolve(@objfun,x0)ô£¬ÕâÑùºÃÏñ²»Ðа¡...

Êǽ«x¸Ä³É(:,3)
4Â¥2015-02-10 22:03:50
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhaoxyah

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

ÒýÓûØÌû:
4Â¥: Originally posted by CelestialCYJ at 2015-02-10 22:03:50
Êǽ«x¸Ä³É(:,3)...

лл
5Â¥2015-02-10 22:25:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zhaoxyah µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 262Çóµ÷¼Á +7 ÌìϵÚÒ»ÎÄ 2026-04-04 8/400 2026-04-05 21:31 by ¼¤Á÷Ó¶É
[¿¼ÑÐ] 283Çóµ÷¼Á +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by à£à£à£0119
[¿¼ÑÐ] Çóµ÷¼ÁÇóµ÷¼Á +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[¿¼ÑÐ] ±¾¿Æ211£¬293·ÖÇëÇóµ÷¼Á +8 Á«²Ë¾ÍÊÇź°É 2026-04-03 9/450 2026-04-05 19:12 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 0817»¯Ñ§¹¤³ÌÓë¼¼ÊõÇóµ÷¼Á£¬Ò»Ö¾Ô¸Öк£Ñó319 +13 lv945 2026-04-04 13/650 2026-04-05 18:14 by Öí»á·É
[¿¼ÑÐ] 318Çóµ÷¼Á +11 ykyhsa 2026-04-05 13/650 2026-04-05 12:44 by aidndnjck
[¿¼ÑÐ] Çóµ÷¼Á +4 wos666 2026-04-03 4/200 2026-04-05 11:48 by arrow8852
[¿¼ÑÐ] ²ÄÁϹ¤³Ì085601Êý¶þÓ¢Ò»335Çóµ÷¼Á +6 Ë«ÂíβƦÀϰå2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[¿¼ÑÐ] µ÷¼Á +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[ÂÛÎÄͶ¸å] ÇóÎÄÏ× 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[¿¼ÑÐ] ÍÁľ304Çóµ÷¼Á +4 ÍÃͻͻͻ£¬ 2026-03-31 4/200 2026-04-04 13:34 by 1753564080
[¿¼ÑÐ] Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-04 3/150 2026-04-04 12:19 by Éá¶øºóµÃ
[¿¼ÑÐ] 311Çóµ÷¼Á +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[¿¼ÑÐ] ²ÄÁÏר˶322·Ö +13 ¹þ¹þ¹þºðºðºð¹þ 2026-04-01 13/650 2026-04-03 16:08 by ŶŶ123
[¿¼ÑÐ] 085501Ò»Ö¾Ô¸Ì칤´ó£¬»úеר˶Çóµ÷¼Á£¬¿ç²ÄÁÏ +3 33ÉÏ 2026-04-03 3/150 2026-04-03 14:08 by 1753564080
[¿¼ÑÐ] Ò»Ö¾Ô¸Öйú¿ÆÑ§Ôº´óѧ265Çóµ÷¼Á +9 Ìñµ­ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +6 ÏàÐűػá¹ââÍòÕ 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸ÄϾ©Ê¦·¶´óѧ¼ÆËã»úר˶£¬³õÊÔ373£¬Áù¼¶Í¨¹ý£¬ +3 ¼ÆËã»ú×·ÃÎÈË 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[¿¼ÑÐ] 085600 Ò»Ö¾Ô¸9 ×Ü·Ö351 Çóµ÷¼ÁѧУ +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Ò»Ö¾Ô¸ ÄϾ©º½¿Õº½Ìì´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶ +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û