一区二区成人免费视频网站观看 ,亚洲淫妇被狗狗狂操,后入爆操黑丝美女白丝av

ybnav

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 681.3
帖子: 65
在線(xiàn): 47.6小時(shí)
蟲(chóng)號(hào): 508614
注冊(cè): 2008-02-21
性別: GG
專(zhuān)業(yè): 信息安全

引用回帖:

Originally posted by change0618 at 2011-03-11 20:03:41:
z(t) 是個(gè)啥玩意

不好意思是y(t)

回復(fù)此樓

11樓2011-03-11 22:29:17

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

ybnav

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 681.3
帖子: 65
在線(xiàn): 47.6小時(shí)
蟲(chóng)號(hào): 508614
注冊(cè): 2008-02-21
性別: GG
專(zhuān)業(yè): 信息安全

引用回帖:

Originally posted by zzg_gm at 2011-03-11 20:33:49:
這要用If語(yǔ)句判斷，確定一些量的值然后用dsolve函數(shù)求解！要是想畫(huà)圖，可將求解出的函數(shù)表達(dá)式用plot命令畫(huà)出（因?yàn)槭且粋€(gè)個(gè)函數(shù)畫(huà)圖，只是二維的，所以用plot,而不是plot3),！你再試試，命令的格式一般的書(shū)上都有 ...

您的意思是不是要分段來(lái)做？求出各函數(shù)等于0的那個(gè)點(diǎn)？

贊一下

回復(fù)此樓

12樓2011-03-11 22:31:03

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

ybnav

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 681.3
帖子: 65
在線(xiàn): 47.6小時(shí)
蟲(chóng)號(hào): 508614
注冊(cè): 2008-02-21
性別: GG
專(zhuān)業(yè): 信息安全

引用回帖:

Originally posted by xiegangmai at 2011-03-11 15:42:54:
把你的函數(shù)貼出來(lái)，我稍改改就可以求解了。

或者參考：http://www.gaoyang168.com/bbs/viewthread.php?tid=2764636###

把微分方程寫(xiě)成分段函數(shù)的形式即可。

我在matlab里是這樣做的：
如果不考慮f1(t)、f2(t)、f3(t)，令它們的值都等于1
函數(shù)
function dydt=xyz(t,y)

dydt=[  -y(5)
      -0.1*y(2)+y(5)
      0.1*y(2)
      -y(2)+y(5)+y(5)+0.5*y(5)
      y(2)-y(5)-y(5)-0.5*y(5)
   ];

腳本
clear all
close all

clc

tspan = [0 60];
y0 = [1000;1;0;200;0];

[t,Y]=ode45(@xyz,tspan,y0);
plot(t,Y(:,1),'r-','LineWidth',1.5);
hold on;
plot(t,Y(:,2),'g-','LineWidth',1.5);
hold on;
plot(t,Y(:,3),'b.','LineWidth',1.5);
hold on;
plot(t,Y(:,4),'g.','LineWidth',1.5);
hold on;
plot(t,Y(:,5),'r-.','LineWidth',1.5);

axis([0 60 0 1000])
legend('Y(1)','Y(2)','Y(3)','Y(4)','Y(5)')
grid on

這是沒(méi)有加上f1(t)、f2(t)、f3(t)這三個(gè)的情況，如果加上這三個(gè)函數(shù)，應(yīng)該怎么改呢？

贊一下

回復(fù)此樓

13樓2011-03-11 22:37:59

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

change0618

鐵桿木蟲(chóng) (著名寫(xiě)手)

方丈大師

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線(xiàn): 546.7小時(shí)
蟲(chóng)號(hào): 496517
注冊(cè): 2008-01-19
專(zhuān)業(yè): 化學(xué)反應(yīng)工程

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-03-11 23:43:01

我覺(jué)得你這個(gè)方程的主要問(wèn)題有：
在t=24.52之后x就開(kāi)始小于0了，也就是說(shuō)x=0時(shí)t在24.52-24.53之間，但是又出不來(lái)真實(shí)值。
假設(shè)在某t時(shí)x=0了，則f2(t)=0了，dxdt=-a*f2*v+a*f1*y+a*f3*y+c*y 卻轉(zhuǎn)向變成了正的，x又開(kāi)始變正了，那f2(t)就有等于1了,此時(shí)dx/dt確實(shí)負(fù)的，x又開(kāi)始降低變負(fù).于是病態(tài)了

贊一下(2人)

回復(fù)此樓

14樓2011-03-11 23:39:52

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

xiegangmai

版主 (職業(yè)作家)

我沒(méi)頭銜

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
散金: 2394
紅花: 63
沙發(fā): 4
帖子: 4046
在線(xiàn): 1330.2小時(shí)
蟲(chóng)號(hào): 653607
注冊(cè): 2008-11-13
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 仿真模擬

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by ybnav at 2011-03-11 22:37:59:
我在matlab里是這樣做的：
如果不考慮f1(t)、f2(t)、f3(t)，令它們的值都等于1
函數(shù)
function dydt=xyz(t,y)

dydt=[  -y(5)
      -0.1*y(2)+y(5)
      0.1*y(2)
      -y(2)+y(5)+y(5)+0.5 ...

如樓上分析的，到24后就算不下去了。

CODE:

function ybnav

clc

clear

tspan = [ 0 60 ];

y0 = [ 1000; 1; 0; 200; 0 ];

figure

[ t, Y ] = ode45( @xyz, tspan, y0 );

plot( t, Y( :,1 ), 'r-', 'LineWidth', 1.5 );

hold on;

plot( t, Y( :, 2 ), 'g-', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 3 ), 'b.', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 4 ), 'g.', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 5 ), 'r-.', 'LineWidth', 1.5 );

axis( [ 0 60 0 1000 ] )

legend( 'Y(1)', 'Y(2)', 'Y(3)', 'Y(4)', 'Y(5)' )

grid on 

figure

tspan = [ 0 24 ];

[ t, Y ] = ode45( @xyz1, tspan, y0 );

plot( t, Y( :,1 ), 'r-', 'LineWidth', 1.5 );

hold on;

plot( t, Y( :, 2 ), 'g-', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 3 ), 'b.', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 4 ), 'g.', 'LineWidth', 1.5 );

% hold on;

plot( t, Y( :, 5 ), 'r-.', 'LineWidth', 1.5 );

axis( [ 0 60 0 1000 ] )

legend( 'Y(1)', 'Y(2)', 'Y(3)', 'Y(4)', 'Y(5)' )

grid on 

function dydt = xyz( t, y )

dydt = zeros( 5, 1 );

dydt = [  -y( 5 )

        -0.1 * y( 2 ) + y( 5 )

         0.1 * y( 2 )

        - y( 2 ) + y( 5 ) + y( 5 ) + 0.5 * y( 5 )

         y( 2 ) - y( 5 ) - y( 5 ) - 0.5 * y( 5 )

     ];

 function dydt = xyz1( t, y )

dydt = zeros( 5, 1 );

dydt = [  -y( 5 ) * 0 * ( y( 1 ) == 0 ) + -y( 5 ) * 1 * ( y( 1 ) > 0 ) 

        -0.1 * y( 2 ) + y( 5 ) * 0 * ( y( 1 ) == 0 ) + y( 5 ) * 1 * ( y( 1 ) > 0 )

         0.1 * y( 2 )

        - y( 2 ) * 0 * ( y( 4 ) == 0 ) - y( 2 ) * 1 * ( y( 4 ) > 0 ) + y( 5 ) * 0 * ( y( 1 ) == 0 ) + y( 5 ) * 1 * ( y( 1 ) > 0 )  + y( 5 )  * 0 * ( y( 2 ) == 0 ) + y( 5 ) * 1 * ( y( 2 ) > 0 )+ 0.5 * y( 5 )

         y( 2 ) - y( 5 ) - y( 5 ) - 0.5 * y( 5 )

     ];

贊一下(1人)

回復(fù)此樓

明德厚學(xué)、求是創(chuàng)新

15樓2011-03-11 23:43:47

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

ybnav

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 681.3
帖子: 65
在線(xiàn): 47.6小時(shí)
蟲(chóng)號(hào): 508614
注冊(cè): 2008-02-21
性別: GG
專(zhuān)業(yè): 信息安全

感謝兩位，這是一篇文獻(xiàn)上的，我想還原一下，因?yàn)槲乙黄撐睦镎靡部梢杂蒙�，我再好好檢查檢查方程，看看是哪里的問(wèn)題，非常感謝！

贊一下

回復(fù)此樓

16樓2011-03-12 00:08:14

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

信彼南山

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線(xiàn): 233.5小時(shí)
蟲(chóng)號(hào): 1133529
注冊(cè): 2010-10-27
專(zhuān)業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

★
xiegangmai(金幣+1): 鼓勵(lì)討論交流！用Simulink是可以解的，也簡(jiǎn)便，但也需要方程不奇異。 2011-03-12 13:27:23

simulink一下試試

贊一下(1人)

回復(fù)此樓

17樓2011-03-12 11:01:03

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

信彼南山

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線(xiàn): 233.5小時(shí)
蟲(chóng)號(hào): 1133529
注冊(cè): 2010-10-27
專(zhuān)業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

★ ★ ★ ★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiegangmai(金幣+5): 謝謝應(yīng)助！辛苦了！ 2011-03-12 16:56:10

simulink做了一下

贊一下(2人)

回復(fù)此樓

18樓2011-03-12 15:17:23

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

信彼南山

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線(xiàn): 233.5小時(shí)
蟲(chóng)號(hào): 1133529
注冊(cè): 2010-10-27
專(zhuān)業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

f1和f2一直都是1
f3有變化

贊一下

回復(fù)此樓

19樓2011-03-12 15:18:14

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

ybnav

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 681.3
帖子: 65
在線(xiàn): 47.6小時(shí)
蟲(chóng)號(hào): 508614
注冊(cè): 2008-02-21
性別: GG
專(zhuān)業(yè): 信息安全

只會(huì)簡(jiǎn)單的matlab操作，simulink還沒(méi)接觸過(guò)，要好好學(xué)習(xí)學(xué)習(xí)，感謝信彼南山大力相助

贊一下

回復(fù)此樓

20樓2011-03-12 23:09:43

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 290求調(diào)劑 +5	孔志浩 2026-03-12	10/500	2026-03-16 09:01 by 余暉&
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專(zhuān)家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5	瑪卡巴卡啊哈 2026-03-11	5/250	2026-03-13 21:40 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 321求調(diào)劑（食品/專(zhuān)碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321
[考博] 2026年博士申請(qǐng) +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43
[考研] 求調(diào)劑材料專(zhuān)碩293 +6	段_(:з」∠)_ 2026-03-10	6/300	2026-03-10 18:22 by ms629