| 查看: 509 | 回复: 1 | |||
[交流]
【求助】用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] |
» 猜你喜欢
关于豆爷回答的JTJC与%2F数量
已经有4人回复
基金中了
已经有11人回复
filecode
已经有18人回复
听说今天filecode变了
已经有44人回复
固定端突然变了,今天
已经有8人回复
filecode
已经有11人回复
大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑
已经有6人回复
一个有机合成实验室都需要哪些设备?
已经有12人回复
求各位大神看下
已经有33人回复
咨询面上基金
已经有7人回复
» 抢金币啦!回帖就可以得到:
1-甲基-1,2-双(三甲基硅氧基)乙烷(CAS:17887-27-3)
+1/87
1,2-双(三甲基硅氧基)乙烷(CAS:7381-30-8 )
+1/87
乙二醇二甲基环硅醚(CAS:57915-66-9)
+1/87
三(三甲基硅烷)磷酸酯(CAS:10497-05-9)
+1/84
四(三甲基硅基)赤藓醇醚(CAS:18547-29-0 )
+1/83
中国科学技术大学-苏州国家实验室联培招生-本科直博/硕博连读/申请考核三种方式
+1/82
实验室富余
+1/80
坐标重庆 诚征男友
+1/63
哈工大(深圳)何自开团队诚招化学工程与技术专业博士生
+1/45
哈工大(深圳)何自开团队诚招2027级化学工程与技术专业博士生
+1/45
2026年量子算法、密码学与数据分析国际会议(QACDA 2026)
+1/37
上海交通大学诚聘有机化学方向博士生/科研助理
+1/29
澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域
+1/12
北理工集成电路杰青团队 | 诚招科助理
+1/10
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/9
有合作举办scopus会议的吗?
+1/6
北理工集成电路杰青团队 | 诚招科助理
+1/5
高弯曲强度玻璃制品的特殊加工工艺
+1/4
北京大学临床医学高等研究院袁琳洁课题组诚招科研助理(推荐读博)
+1/3
致每一位等待国自然结果的科研人:我们想听听你的心声
+1/1
★ ★
robert2020(金币+2):多谢应助!欢迎常来仿真模拟版交流! 2010-11-09 22:40:19
天边5328(金币+100):谢谢! 2010-11-15 13:05:07
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上运行一下吧 |
2楼2010-11-09 08:49:32










回复此楼