版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (3815)
>基金申请 (1364)
>论文道贺祈福 (517)
>导师招生 (292)
>招聘信息布告栏 (228)
>论文投稿 (169)
>虫友互识 (139)
>学术会议 (132)
>文献求助 (126)
>硕博家园 (73)
>考研 (53)
>休闲灌水 (51)
>教师之家 (43)
>考博 (41)
>公派出国 (26)
>博后之家 (25)

返回列表

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

回复此楼

» 猜你喜欢

filecode变化情况已经有16人回复
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑已经有7人回复
关于豆爷回答的JTJC与%2F数量已经有6人回复
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！已经有31人回复
filecode 已经有19人回复
基金中了已经有11人回复
听说今天filecode变了已经有44人回复
固定端突然变了，今天已经有8人回复
filecode 已经有11人回复
求各位大神看下已经有33人回复

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

查看全部散金贴

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