| 3 | 1/1 | 返回列表 |
| 查看: 396 | 回復: 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 %% |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 招08考數(shù)學 +3 | laoshidan 2026-03-20 | 7/350 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 材料學碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 070300化學319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|