| ²é¿´: 1148 | »Ø¸´: 3 | |||
wpynwpuľ³æ (ÕýʽдÊÖ)
|
[ÇóÖú]
²ÉÓÃmatlabÇó½â΢·Ö·½³ÌÓöµ½ÎÊÌ⣬Çë°ïæ·ÖÎöһϴíÎó-ÖØÐ»
|
|
ÎÊÌâÊÇ£ºÒºµÎÔ˶¯µÄ΢·Ö·½³Ì£»ÒºµÎ³õʼËÙ¶È10m/s£¬Çó½âʱ¼äʱ0.15s£¬²ÉÓÃÁú¸ñ-¿âËþ·¨Çó½â£¬³ÌÐòÈçÏ£º function y=fly1(t,U) a=1.1616; %µªÆøÃÜ¶È b=1.873e-5; %µªÆøð¤¶È d=3e-4; %ÈÛµÎÖ±¾¶ c=2385; %ÈÛµÎÃÜ¶È g=9.8; %ÖØÁ¦¼ÓËÙ¶È f=d*a*U/b %ReÊý h=0.28+6/f.^0.5+21/f %ÍÏҷϵÊý y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c) options=odeset('outputfcn','odeplot') [t,U]=ode45('fly1',[0:0.001:0.15],[0;10],options) U1=U(:,1) plot(t,U1) title('Droplet Velocity VS Time') xlabel('Time/s') ylabel('Droplet Velocity/m/s') ³öÏֵĴíÎóÊÇ£º ??? Maximum recursion limit of 2000 reached. Use set(0,'RecursionLimit',N) to change the limit. Be aware that exceeding your available stack space can crash MATLAB and/or your computer. Error in ==> odeset Çë´ó¼Ò°ï°ïæ·ÖÎöһϴíÎóÔÒò |
» ²ÂÄãϲ»¶
281Çóµ÷¼Á£¨0805£©
ÒѾÓÐ15È˻ظ´
085600²ÄÁÏÓ뻯¹¤µ÷¼Á 324·Ö
ÒѾÓÐ8È˻ظ´
ÄÜÔ´²ÄÁÏ»¯Ñ§¿ÎÌâ×éÕÐÊÕ˶ʿÑо¿Éú8-10Ãû
ÒѾÓÐ10È˻ظ´
267Ò»Ö¾Ô¸ÄϾ©¹¤Òµ´óѧ0817»¯¹¤Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
274Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑоÀú
ÒѾÓÐ10È˻ظ´
328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑоÀú
ÒѾÓÐ7È˻ظ´
²ÄÁÏרҵÇóµ÷¼Á
ÒѾÓÐ5È˻ظ´
Ò»Ö¾Ô¸Ìì½ò´óѧ»¯Ñ§¹¤ÒÕרҵ£¨081702£©315·ÖÇóµ÷¼Á
ÒѾÓÐ10È˻ظ´
344Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
¸ßÊÖ½øÀ´ÓÃmatlabÇó½â³£Î¢·Ö·½³Ì£¬Ê®·Ö¸Ðл£¡
ÒѾÓÐ17È˻ظ´
matlab΢·Ö·½³Ì×é²ÎÊýÄâºÏ£¬ÒÔÖÜΪµ¥Î»Çó½â£¬»ã×ܺóÒÔÄêΪµ¥Î»½øÐÐÊýÖµ±È½Ï
ÒѾÓÐ4È˻ظ´
matlab ÏßÐÔ³£Î¢·Ö·½³ÌÇó½â
ÒѾÓÐ5È˻ظ´
matlab½â΢·Ö·½³Ì
ÒѾÓÐ10È˻ظ´
MatlabÇó½âƫ΢·Ö·½³Ì×é
ÒѾÓÐ22È˻ظ´
matlabÊýÖµÇó½â±ß½çÌõ¼þ΢·Ö·½³Ì×é
ÒѾÓÐ7È˻ظ´
¸ß½ðÇóÖúmatlab½â΢·Ö·½³Ì×é
ÒѾÓÐ12È˻ظ´
matlab½â΢·Ö·½³Ì×é
ÒѾÓÐ15È˻ظ´
¡¾ÇóÖú¡¿matlab Çó½â΢·Ö·½³ÌÖеÄδ֪²ÎÊý
ÒѾÓÐ20È˻ظ´
¡¾ÇóÖú¡¿Ïò¸÷λ´óÏÀÇóÖúmatlabÇó½â΢·Ö·½³Ì×éÓöµ½µÄÒ»¸öÎÊÌâ
ÒѾÓÐ21È˻ظ´
¡¾ÇóÖú¡¿MATLAB ÓÐÏÞ²î·Ö·¨£¨FDM£©Çó½âƫ΢·Ö·½³Ì
ÒѾÓÐ22È˻ظ´
¡¾ÇóÖú¡¿»ý·Ö΢·Ö·½³ÌmatlabÇó½â
ÒѾÓÐ6È˻ظ´
¡¾ÇóÖú¡¿matlabÔõôÇó½âƫ΢·Ö·½³Ì×é°¡£¬ÏÈллÁË
ÒѾÓÐ13È˻ظ´
change0618
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
·½ÕÉ´óʦ
- Ó¦Öú: 44 (СѧÉú)
- ½ð±Ò: 17724.5
- ºì»¨: 17
- Ìû×Ó: 2413
- ÔÚÏß: 546.7Сʱ
- ³æºÅ: 496517
- ×¢²á: 2008-01-19
- רҵ: »¯Ñ§·´Ó¦¹¤³Ì
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ...
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
wpynwpu: ½ð±Ò+50, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ÕæÊǸßÊÖ£¬Ð»Ð» 2012-05-26 11:04:43
³ôË®¹µ: ½ð±Ò+3, лл½»Á÷~~ 2012-05-26 11:06:04
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
wpynwpu: ½ð±Ò+50, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ÕæÊǸßÊÖ£¬Ð»Ð» 2012-05-26 11:04:43
³ôË®¹µ: ½ð±Ò+3, лл½»Á÷~~ 2012-05-26 11:06:04
|
function exercise11 options=odeset('outputfcn','odeplot'); [t,U]=ode45(@fly1,[0:0.001:0.15],10,options); U1=U(:,1); plot(t,U1) title('Droplet Velocity VS Time') xlabel('Time/s') ylabel('Droplet Velocity/m/s') function y=fly1(t,U) a=1.1616; %µªÆøÃÜ¶È b=1.873e-5; %µªÆøð¤¶È d=3e-4; %ÈÛµÎÖ±¾¶ c=2385; %ÈÛµÎÃÜ¶È g=9.8; %ÖØÁ¦¼ÓËÙ¶È f=d*a*U/b; %ReÊý h=0.28+6/f.^0.5+21/f; %ÍÏҷϵÊý y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c); |
2Â¥2012-05-26 10:51:19
wpynwpu
ľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 4 (Ó×¶ùÔ°)
- ½ð±Ò: 1469.1
- ºì»¨: 1
- Ìû×Ó: 550
- ÔÚÏß: 262.5Сʱ
- ³æºÅ: 438036
- ×¢²á: 2007-09-13
- ÐÔ±ð: GG
- רҵ: ½ðÊô²ÄÁϵÄÄý¹ÌÓë½á¾§Ñ§
3Â¥2012-05-26 11:04:02
change0618
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
·½ÕÉ´óʦ
- Ó¦Öú: 44 (СѧÉú)
- ½ð±Ò: 17724.5
- ºì»¨: 17
- Ìû×Ó: 2413
- ÔÚÏß: 546.7Сʱ
- ³æºÅ: 496517
- ×¢²á: 2008-01-19
- רҵ: »¯Ñ§·´Ó¦¹¤³Ì
4Â¥2012-05-26 12:31:24













»Ø¸´´ËÂ¥