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]

回复此楼

» 猜你喜欢

体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有18人回复
面上可以超过30页吧？已经有7人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有5人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有6人回复
版面费该交吗已经有13人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有10人回复
什么是人一生最重要的？已经有4人回复

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

查看全部散金贴

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人)