24小时热门版块排行榜

>论坛更新日志 (3464)
>文献求助 (361)
>虫友互识 (355)
>导师招生 (196)
>考博 (128)
>硕博家园 (101)
>休闲灌水 (95)
>博后之家 (91)
>论文投稿 (83)
>公派出国 (68)
>基金申请 (61)
>教师之家 (40)
>找工作 (40)
>考研 (32)
>有机交流 (31)
>绿色求助(高悬赏) (22)

返回列表

天边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]

回复此楼

» 猜你喜欢

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

查看全部散金贴

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