版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

wave0403

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1033
帖子: 174
在線: 260小時(shí)
蟲號(hào): 310161
注冊(cè): 2006-12-23
性別: GG
專業(yè): 化學(xué)動(dòng)力學(xué)

[求助] 積分并擬合實(shí)驗(yàn)數(shù)據(jù)

最近用matlab處理一個(gè)實(shí)驗(yàn)數(shù)據(jù)，是極坐標(biāo)下的.dat文件（以附件為例），文件格式：（r, theta, I）。這里半徑0.5<r<237.5, 9.5度<theta<89度，I代表（r, theta）的強(qiáng)度。我想對(duì)任意兩個(gè)半徑0.5<r1<r2<237.5強(qiáng)度求和得intensity（theta），然后對(duì)這個(gè)intensity隨theta變化的曲線進(jìn)行擬合并最后做誤差分析，具體擬合公式見代碼，編了下面代碼，但是運(yùn)行不了。煩請(qǐng)高手指正錯(cuò)誤，謝謝了。

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);

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 1_polar.dat

2014-10-24 15:58:55, 1.53 M

» 猜你喜歡

299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(fù)

1樓 2014-10-24 15:45:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

用含有定積分的隱函數(shù)擬合數(shù)據(jù)，matlab或者origin 已經(jīng)有18人回復(fù)
數(shù)據(jù)關(guān)系的擬合已經(jīng)有9人回復(fù)
跪求高手幫忙看看這組數(shù)據(jù)如何進(jìn)行擬合已經(jīng)有7人回復(fù)
怎樣擬合數(shù)據(jù)求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D參數(shù) 已經(jīng)有18人回復(fù)
origin 非線性自定義擬合問題~~方程組編寫及運(yùn)行~~ 已經(jīng)有8人回復(fù)
速率方程的問題求教！已經(jīng)有4人回復(fù)
通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復(fù)
數(shù)據(jù)擬合多項(xiàng)式已經(jīng)有6人回復(fù)
對(duì)數(shù)據(jù)進(jìn)行曲線擬合后如何得到擬合誤差已經(jīng)有30人回復(fù)
用matlab擬合雙因變量實(shí)驗(yàn)結(jié)果求助已經(jīng)有7人回復(fù)
怎么由測(cè)試數(shù)據(jù)得到擬合方程已經(jīng)有11人回復(fù)
實(shí)驗(yàn)數(shù)據(jù)的數(shù)學(xué)模型構(gòu)建、曲線擬合方面問題求助已經(jīng)有6人回復(fù)
【一個(gè)數(shù)據(jù)擬合的問題，比較麻煩】已經(jīng)有17人回復(fù)
【求助】大量數(shù)據(jù)程序化擬合？已經(jīng)有9人回復(fù)
【求助】用ZsimpWin擬合數(shù)據(jù)后如何將擬合的數(shù)據(jù)導(dǎo)出并在origin中做到一張圖上呢？已經(jīng)有10人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wave0403 的主題更新