| 查看: 1100 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xiaowei7032银虫 (小有名气)
|
[交流]
【求助】帮忙编写一个用matlab程序进行曲线二次拟合的程序【已解决】
|
||
|
麻烦大家帮忙编写一个用matlab程序进行曲线二次拟合的程序，提供的两组拟合数据如下： 横轴——人均GDP（元）：1771.264、1729.482、1994.071、2343.039、2652.316、3000.924、3352.7、3710.833、3992.457、4319.191、4622.635、4996.076、5438.174、5909.201、6689.266、7348.099 纵轴——污染物排放强度（吨/万元）：5.09、4.91、4.48、3.97、3.79、3.47、3.12、2.94、2.75、2.51、2.4、2.54、2.48、2.3、2.32、2.24 注意：回帖要有运行的程序，运行结果即拟合图（一定要有相关系数），麻烦大家了。 因为很急这用，11月3号前再回帖就无效了 补充：一定要用MATLAB，就想知道matlab怎么进行二次拟合 [ Last edited by nono2009 on 2009-12-3 at 20:25 ] |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有5人回复
-
重要来源：本周末出结果
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2009-10-30 23:25:35
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-10-30 23:16:22
xiaowei7032
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 915.1
- 红花: 1
- 帖子: 169
- 在线: 49.6小时
- 虫号: 601119
- 注册: 2008-09-12
- 性别: MM
- 专业: 污染物行为过程及其环境效
4楼2009-10-30 23:36:47
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xiaowei7032(金币+10,VIP+0): 10-30 23:50
nono2009(金币+2,VIP+0):谢谢专家！ 10-31 07:32
nono2009(金币+9,VIP+0):代发应助奖励！ 12-3 20:25
xiaowei7032(金币+10,VIP+0): 10-30 23:50
nono2009(金币+2,VIP+0):谢谢专家！ 10-31 07:32
nono2009(金币+9,VIP+0):代发应助奖励！ 12-3 20:25
|
>> x=[1771.264 1729.482 1994.071 2343.039 2652.316 3000.924 3352.7 3710.833 3992.457 4319.191 4622.635 4996.076 5438.174 5909.201 6689.266 7348.099]; >> y=[5.09 4.91 4.48 3.97 3.79 3.47 3.12 2.94 2.75 2.51 2.4 2.54 2.48 2.3 2.32 2.24]; >> p = polyfit(x,y,2); >> yy=polyval(p,x); >> Sy=sumsqr(y-mean(y))/16; >> Se=sumsqr(y-yy)/16; >> Rs=1-Se/Sy; >> scatter(x,y);hold on;plot(x,yy);title(['{\itR}^{2}=' num2str(Rs)]) >> |
5楼2009-10-30 23:47:07
