| 5 | 1/1 | 返回列表 |
| 查看: 944 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
zhaoqian59捐助貴賓 (小有名氣)
|
[求助]
這個(gè)程序是不是哪里有問題?初始值跟實(shí)驗(yàn)值是一樣的,但是模擬結(jié)果不對(duì) 已有1人參與
|
||
|
f(x) pl=1000; phyd=110000; P0=29.39E9; a0=0.0341; v=1.5; x1=x(1); x2=x(2); pg=P0*(a0/x1)^(3*v); y=[x2 (1/pl*(pg-phyd)-3/2*x2^2)/x1]'; Ts = 0.001; % 仿真步長(zhǎng) Tn = 1; % 仿真終止時(shí)間 t = 0:Ts:Tn; % 仿真時(shí)間范圍 N = length(t); x=[0.034 0]'; % 系統(tǒng)初值:x=[a a的導(dǎo)數(shù)] %% Variable Declaration and Initialization a=0.034*ones(1,N); a_dt=zeros(1,N); %x = [a a_d]'; % intial states %% Simulation Iteration for n = 1:N a(n)=x(1); a_dt(n)=x(2); k1 = fx(x); k2 = fx(x+Ts/2*k1); k3 = fx(x+Ts/2*k2); k4 = fx( x+Ts *k3); x = x + Ts/6 * ( k1 + 2*k2 + 2*k3 + k4 ); end %% Plot Response plot(t,a) xlabel('t/(s)') ylabel('a') 是一個(gè)二階非線性常微分方程 |
matlab編程繪圖 |

捐助貴賓 (小有名氣)

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

捐助貴賓 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|