| 2 | 1/1 | 返回列表 |
| 查看: 457 | 回復(fù): 1 | |||
[交流]
【求助】用matlab給出參數(shù)k的具體數(shù)值
|
|
請高手幫忙根據(jù)下面(1)、(2)兩組數(shù)據(jù),運(yùn)用matlab編程求解兩個k值.我對matlab完全不懂,麻煩高手給出參數(shù)k的具體數(shù)值。非常感謝!!! 方程:y=y0+(ylimt-y0)/(2*c0)*(x+1/k+c0-((x+c0+1/k)^2-4*c0*x)^0.5) 其中 (1) y0=1.77235; ylimt=0.93066; c0=0.00002 xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125 0.00145] ydata=[1.72427 1.58145 1.37302 1.23911 1.18416 1.13857 1.10483 1.04227 1.02577 0.98729 0.97012 0.96478 0.96057 0.9512 0.94087 0.93066] (2) y0=1.17094; ylimt=0.30983; c0=0.000015 xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125 0.00145 0.00155] ydata=[1.1414 1.08394 1.02934 1.00169 0.98691 0.92352 0.91823 0.90804 0.83408 0.80894 0.75992 0.71518 0.64931 0.64825 0.56648 0.34087 0.30983] |
» 搶金幣啦!回帖就可以得到:
+5/305
+1/85
+1/84
+1/78
+2/76
+1/66
+1/53
+1/32
+1/30
+1/17
+1/16
+1/10
+1/7
+1/6
+1/3
+1/3
+1/2
+1/1
+1/1
+1/1
|
先換元,令X=(x+1/k+c0-((x+c0+1/k)^2-4*c0*x)^0.5)有X=(2*c0*(y0data-ydata))./(ylimtdata-ydata) k=1./(((2*c0*xdata)+X.*X)/(2*X)-xdata-c0*ones(1,16)) 下面是代碼: 第一組數(shù)據(jù): y0=1.77235; ylimt=0.93066; c0=0.00002 ; xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125] ydata=[1.72427 1.58145 1.37302 1.23911 1.18416 1.13857 1.10483 1.04227 1.02577 0.98729 0.97012 0.96478 0.96057 0.9512 0.94087] y0data=ones(1,15)*y0; ylimtdata=ones(1,15)*ylimt; X=(2*c0*(y0data-ydata))./(ylimtdata-ydata) X.*X k=1./(((2*c0*xdata)+X.*X)/(2*X)-xdata-c0*ones(1,15)) 第二組數(shù)據(jù) y0=1.17094; ylimt=0.30983; c0=0.000015; xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125 0.00145] ydata=[1.1414 1.08394 1.02934 1.00169 0.98691 0.92352 0.91823 0.90804 0.83408 0.80894 0.75992 0.71518 0.64931 0.64825 0.56648 0.34087] y0data=ones(1,16)*y0; ylimtdata=ones(1,16)*ylimt; X=(2*c0*(y0data-ydata))./(ylimtdata-ydata) k=1./(((2*c0*xdata)+X.*X)/(2*X)-xdata-c0*ones(1,16)) 自己在MATLAB上運(yùn)行一下吧 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +11 | 煙汐憶海 2026-03-16 | 21/1050 |
|
|---|---|---|---|---|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +3 | sbdksD 2026-03-19 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|