24小时热门版块排行榜

返回列表

天边5328

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 890.1
帖子: 104
在线: 43.8小时
虫号: 698719

[交流] 【求助】用matlab给出参数k的具体数值

请高手帮忙根据下面(1)、(2)两组数据,运用matlab编程求解两个k值.我对matlab完全不懂，麻烦高手给出参数k的具体数值。非常感谢!!!
方程：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]

回复此楼

» 猜你喜欢

请教限项目规定已经有4人回复
拟解决的关键科学问题还要不要写已经有8人回复
存款400万可以在学校里躺平吗已经有26人回复
最失望的一年已经有7人回复
推荐一本书已经有16人回复
国自然申请面上模板最新2026版出了吗？已经有20人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
疑惑？已经有5人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-09 04:26:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snowman09

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1619.1
帖子: 148
在线: 125.1小时
虫号: 1060657

★ ★
robert2020(金币+2):多谢应助！欢迎常来仿真模拟版交流！ 2010-11-09 22:40:19
天边5328(金币+100):谢谢！ 2010-11-15 13:05:07

先换元，令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))
下面是代码：
第一组数据：
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))

第二组数据

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上运行一下吧

赞一下(1人)