| 查看: 417 | 回复: 0 | ||
[求助]
积分并拟合实验数据
|
|
最近用matlab处理一个实验数据，是极坐标下的.dat文件（以附件为例），文件格式：（r, theta, I）。这里半径0.5<r<237.5, 9.5度<theta<89度，I代表（r, theta）的强度。我想对任意两个半径0.5<r1<r2<237.5强度求和得intensity（theta），然后对这个intensity随theta变化的曲线进行拟合并最后做误差分析，具体拟合公式见代码，编了下面代码，但是运行不了。烦请高手指正错误，谢谢了。 function b = intensity(r1,r2,'1_polar.dat') start_r = ((r1/0.5)-1)*160+1; end_r = ((r2/0.5)-1)*160 +1; global intensity theta a = load('1_polar.dat '); p = 1; % Loop for every value of theta for i = 9.5:0.5:89 j = 1; for k = start_r:160:end_r r(j) = a(k,1); I(j) = a(k,3); j = j+1; end %Method of Trapezoidal integration is used for integrating the intenstiy values intensity(p) = trapz(r',I); theta(p) = i; clear r I %increment radius bounds for getting next values start_r = start_r +1; end_r = end_r + 1; p = p+1; end %Guess for two parameters C, beta guess = [1 1]; param = fminsearch(@fun, guess); c = param(1); beta = param(2); %Generate the fitted values yfit = c.*(1+b.* (5*(cosd(theta)).^5 - 1)); plot(theta,intensity, 'o',theta,yfit,'-'); function sse = fun(param) global intensity theta c = param(1); b = param(2); %error and sum of squared error error = intensity - c.*(1+b.* (5*(cosd(theta)).^5 - 1)); sse = sum(error.^2); |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 1_polar.dat
2014-10-24 15:58:55, 1.53 M
» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有12人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
国社科系统bug了，是不是要放榜了？
已经有6人回复
-
申博发邮件
已经有6人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
用含有定积分的隐函数拟合数据，matlab或者origin
已经有18人回复
- 数据关系的拟合
已经有9人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- origin 非线性自定义拟合问题~~方程组编写及运行~~
已经有8人回复
- 速率方程的问题求教！
已经有4人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- matlab数据拟合问题求助
已经有7人回复
- 数据拟合多项式
已经有6人回复
- 对数据进行曲线拟合后如何得到拟合误差
已经有30人回复
- 用matlab拟合双因变量实验结果求助
已经有7人回复
- 怎么由测试数据得到拟合方程
已经有11人回复
- 实验数据的数学模型构建、曲线拟合方面问题求助
已经有6人回复
- 【一个数据拟合的问题，比较麻烦】
已经有17人回复
- 【求助】大量数据程序化拟合？
已经有9人回复
- 【求助】用ZsimpWin拟合数据后 如何将拟合的数据导出 并在origin中做到一张图上呢？
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
