| ²é¿´: 1167 | »Ø¸´: 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 Çë´ó¼Ò°ï°ïæ·ÖÎöһϴíÎóÔÒò |
» ²ÂÄãϲ»¶
µ÷¼Á
ÒѾÓÐ17È˻ظ´
347²ÄÁÏר˶Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
Çó¿¼ÑвÄÁϵ÷¼Á
ÒѾÓÐ3È˻ظ´
ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£©×ö¹ý·Ö×ÓʵÑé
ÒѾÓÐ6È˻ظ´
¼ÆËã»ú408£üÔÚУ¶à´Î¹ú¼Ò¼¶¾ºÈü»ñ½±£üÉêÇëµ÷¼Á
ÒѾÓÐ4È˻ظ´
388Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
288Çóµ÷¼Á
ÒѾÓÐ12È˻ظ´
Çóµ÷¼Á
ÒѾÓÐ29È˻ظ´
266µ÷¼Á
ÒѾÓÐ6È˻ظ´
304Çóµ÷¼Á
ÒѾÓÐ12È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
¸ßÊÖ½øÀ´ÓÃ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














»Ø¸´´ËÂ¥