| 查看: 534 | 回复: 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] |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
★ ★
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
5