| 查看: 1627 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
qq295118688银虫 (小有名气)
|
[求助]
求助高人：使用Matlab或Origin进行非线性拟合的方法已有1人参与
|
|
|
求助，现在已知公式如下： y=m*k*c*x/((1-k*x)*(1-k*x+c*k*x)) 其中m,k,c为待求参数， 请问怎么使用Matlab或者Origin求的以上三个参数？ 感谢！ 数据如下 x. y 0 0 5 2.45198 10 3.41962 15 4.25807 20 4.98196 25 5.69367 30 6.31764 40 7.78736 50 9.24003 60 10.73413 70 12.10637 75 13.16662 80 14.33655 85 16.40343 90 19.0894 95 22.68451 98 25.90914 感恩！ |
回复此楼» 猜你喜欢
北伊利诺伊大学(Northern Illinois University)机械系招收全奖博士生
已经有0人回复
-
第五届粉末衍射晶体结构解析及Rietveld精修高级研讨会
已经有45人回复
-
物理学I论文润色/翻译怎么收费?
已经有179人回复
-
FDTD软件里的时域图像画图
已经有2人回复
-
复杂量子系统的新物态和新机理讲习班
已经有0人回复
-
导电原子力显微镜cAFM如何测试mcAFM(即磁导AFM)性能？
已经有0人回复
-
感谢小木虫，已完成目标
已经有0人回复
-
零点能
已经有2人回复
-
科研求助
已经有0人回复
-
【25申博】博导推荐求助
已经有5人回复
-
青岛大学李山东教授招生
已经有2人回复
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
qq295118688: 金币+30, ★★★★★最佳答案 2019-07-23 04:08:20
独孤神宇: 金币+5, 鼓励交流 2019-07-23 06:59:42
感谢参与，应助指数 +1
qq295118688: 金币+30, ★★★★★最佳答案 2019-07-23 04:08:20
独孤神宇: 金币+5, 鼓励交流 2019-07-23 06:59:42
|
a=[0 0; 5 2.45198; 10 3.41962; 15 4.25807; 20 4.98196; 25 5.69367; 30 6.31764; 40 7.78736; 50 9.24003; 60 10.73413; 70 12.10637; 75 13.16662; 80 14.33655; 85 16.40343; 90 19.0894; 95 22.68451; 98 25.90914] f=@(k,x) k(1)*k(2)*k(3)*x./((1-k(2)*x).*(1-k(2)*x+k(3)*k(2)*x)) % k(1)=m,k(2)=k,k(3)=c k0=[1 1 1] for ii=1:10 k=lsqcurvefit(f,k0,a(:,1),a(:,2)) k0=k end plot(a(:,1),a(:,2),'o',a(:,1),f(k,a(:,1))) |
» 本帖已获得的红花（最新10朵）
qq295118688
5楼2019-07-23 03:45:31
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
|
后面再加一行： R_sq=1-sum((a(:,2)-f(k,a(:,1))).^2)/sum((a(:,2)-mean(a(:,1))).^2) 算完R^2=0.9994 再给朵花花呗：） |
» 本帖已获得的红花（最新10朵）
7楼2019-07-23 04:03:17
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
9楼2019-07-23 04:16:51
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
11楼2019-07-23 16:36:19
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43440.6
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.8小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
|
a=[0 0; 5 2.45198; 10 3.41962; 15 4.25807; 20 4.98196; 25 5.69367; 30 6.31764; 40 7.78736; 50 9.24003; 60 10.73413; 70 12.10637; 75 13.16662; 80 14.33655; 85 16.40343; 90 19.0894; 95 22.68451; 98 25.90914] f=@(k,x) k(1)*k(2)*x.*((1-k(4)*(1-k(3))*x))./((1+k(2)*x).*(1-k(4)*x)) % k(1)=m,k(2)=K,k(3)=w,k(4)=k k0=[1 1 1 1] for ii=1:10 k=lsqcurvefit(f,k0,a(:,1),a(:,2)) k0=k end plot(a(:,1),a(:,2),'o',a(:,1),f(k,a(:,1))) R_sq=1-sum((a(:,2)-f(k,a(:,1))).^2)/sum((a(:,2)-mean(a(:,1))).^2) |
14楼2019-07-25 18:32:39
50