| 5 | 1/1 | 返回列表 |
| 查看: 399 | 回復: 2 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
legend85銀蟲 (小有名氣)
|
[交流]
【求助】隱函數(shù)求解問題 已有2人參與
|
||
|
各位大蝦,我現(xiàn)在有三組實驗數(shù)據(jù) x=[3.2034 2.5439 2.2442 2.1043 1.9245 1.7846];y=[0.5989 0.4912 0.4593 0.4137 0.3844 0.3573];z=[0.0201 0.0192 0.0191 0.0188 0.0184 0.0180]; 現(xiàn)在我想根據(jù)x和z的值用公式x-(y/(y+z))*((1.6504*y+1.4592*z)/1.6504/0.2969)^1.6504=0對y值進行求解并與實驗值進行比較,應該如何操作并畫出計算值y的趨勢圖???圖以x軸為橫坐標,y為縱坐標!用matlab或origin能行嗎? |

至尊木蟲 (正式寫手)
|
我用二分法寫了個程序,不知是否滿足你的需要: % 隱函數(shù)求解問題 clc;clear;close all; xt=[3.2034 2.5439 2.2442 2.1043 1.9245 1.7846]; zt=[0.0201 0.0192 0.0191 0.0188 0.0184 0.0180]; y=1; for k=1:length(xt); x=xt(k);z=zt(k); y1=0; y2=1; f1=x-(y1/(y1+z))*((1.6504*y1+1.4592*z)/1.6504/0.2969)^1.6504; f2=x-(y2/(y2+z))*((1.6504*y2+1.4592*z)/1.6504/0.2969)^1.6504; while max(abs(f1),abs(f2))>1e-6; f1=x-(y1/(y1+z))*((1.6504*y1+1.4592*z)/1.6504/0.2969)^1.6504; f2=x-(y2/(y2+z))*((1.6504*y2+1.4592*z)/1.6504/0.2969)^1.6504; y=[y1+y2]/2; f=x-(y/(y+z))*((1.6504*y+1.4592*z)/1.6504/0.2969)^1.6504; if f>0; y1=y; else y2=y; end end yt(k)=y; k end yt0=[0.5989 0.4912 0.4593 0.4137 0.3844 0.3573] plot(yt0,'k+');hold on; plot(yt,'rx');hold on; %%%%%%%%%%% 圖解確定根的個數(shù) % figure; for k=1:6; x=xt(1);y=yt(k); ys=linspace(0,1,101); f=x-(ys./(ys+z)).*((1.6504*ys+1.4592*z)/1.6504/0.2969).^1.6504; subplot(2,3,k);plot(ys,f); end %%% over %% |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調劑 +14 | 一樣YWY 2026-04-01 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 284求調劑 +12 | 小熊~~ 2026-03-31 | 12/600 |
|
|
[考研] 一志愿北交大材料工程,總分358 +3 | cs0106 2026-04-01 | 3/150 |
|
|
[考研] 324分 085600材料與化工 +17 | 呆鵝oor 2026-03-27 | 17/850 |
|
|
[考研] 324求調劑 +4 | 想上學求調 2026-04-01 | 5/250 |
|
|
[考研] 考研調劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
[考研] 08工科,295,接受跨專業(yè)調劑 +6 | lmnlzy 2026-03-31 | 6/300 |
|
|
[考研] 318求調劑 +3 | 篤行致遠. 2026-03-31 | 3/150 |
|
|
[考研] 材料工程085601數(shù)二英一335求調劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
[考研] 322求調劑:一志愿湖南大學 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 334求調劑 +7 | Trying] 2026-03-31 | 7/350 |
|
|
[考研] 276求調劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
[考研] 262求調劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 085601一志愿西北工業(yè)大學初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[考研] 303求調劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 一志愿中南大學化學0703總分337求調劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 305求調劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|