[求助]菜鸟级问题：线性拟合怎样做？help里的英文我不太会！

mckof

原始数据：<BR>UH=[-6.22,-8.21,-10.23,-12.21,-14.21,-16.16];<BR>d=0.5;<BR>Is=3.00;<BR>C=0.449;<BR>Im=[0.300,0.400,0.500,0.600,0.700,0.800];<BR>B=Im.*C;<BR>需要拟合的方程RH=(UH.*d)/(Is*B.)<BR>我真的不会啊！希望大家帮帮忙啦！谢谢！

happy

以什么做自变量阿

mckof

不太知道哦！这是一条物理实验题，我用我理解的内容转化出来的。只知道d=0.5;<BR>Is=3.00;<BR>C=0.449;<BR>是固定不变的，自变量是否能有两个啊？如果是的话那就是Ｂ和UH了！

happy

如果是两个自变量的话你给的RH不够<BR>另外看你的数据UH和Im是相关的

mckof

<P>那请问按你认为还差什么呢？</P>

mckof

<P>如果是UH=RH*（（Is*B.）./d），Ｂ是自变量，求系数ＲＨ呢？还有相关系数怎样求啊？谢谢！</P>

[此贴子已经被作者于2005-12-19 19:58:29编辑过]

sorry

polyfit（x,y,n），x,y为要拟合的数据，n为希望最佳拟合数据的多项式的阶次

mckof

能否举例啊？好抽象哦！而且好像和那个表达式没有什么关系似的！

sorry

UH=[-6.22,-8.21,-10.23,-12.21,-14.21,-16.16];<BR>d=0.5;<BR>Is=3.00;<BR>C=0.449;<BR>Im=[0.300,0.400,0.500,0.600,0.700,0.800];<BR>B=Im.*C;<BR>RH=(UH.*d)./(Is.*B);<BR>Pm=polyfit(B,RH,5);<BR><BR>求出来的Pm就你你要拟合的系数<BR>RH1=Pm(6)+Pm(5).*B+Pm(4).*B.^2+Pm(3).*B.^3+Pm(2).*B.^4+Pm(1).*B.^5<BR><BR>

mckof

<P>请问“求出来的Pm就你你要拟合的系数”这句话是什么意思啊？是不是漏了几个字？我看不懂！<BR>另外，介不介意我问多几个弱智的问题？我按你所说出的那样做，得到的结果为：<BR>Pm =</P>
<P>  1.0e+004 *</P>
<P>    0.8617   -1.0788    0.5287   -0.1267    0.0149   -0.0015</P>
<P><BR>RH1 =</P>
<P>   -7.6961   -7.6188   -7.5947   -7.5538   -7.5353   -7.4981<BR><BR>但我只需要一个RH值和一个相关系数r就可以了！请问应该选择哪一个啊？而且线性度的相关系数绝对值是少于1的。希望能帮上忙啦！<BR>谢谢！</P>

happy

<DIV class=quote><B>以下是引用<I>mckof</I>在2005-12-19 21:49:46的发言：</B><BR>
<P>请问“求出来的Pm就你你要拟合的系数”这句话是什么意思啊？是不是漏了几个字？我看不懂！<BR>另外，介不介意我问多几个弱智的问题？我按你所说出的那样做，得到的结果为：<BR>Pm =</P>
<P>  1.0e+004 *</P>
<P>    0.8617   -1.0788    0.5287   -0.1267    0.0149   -0.0015</P>
<P><BR>RH1 =</P>
<P>   -7.6961   -7.6188   -7.5947   -7.5538   -7.5353   -7.4981<BR><BR>但我只需要一个RH值和一个相关系数r就可以了！请问应该选择哪一个啊？而且线性度的相关系数绝对值是少于1的。希望能帮上忙啦！<BR>谢谢！</P></DIV>
<P>线性拟合的话<BR>用Pm=polyfit(B,RH,2);就行了<BR>Pm的到两个数<BR>RH=Pm(2)+Pm(1)*B;<BR></P>

mckof

<P>哦哦！有一点点明白啦！那么那个线性相关系数怎样求呢？</P>

[此贴子已经被作者于2005-12-20 14:34:39编辑过]

cora

corrcoef

mckof

<P>举个例子可以吗？就用上面的材料，我英语超差啊！我对着那些英文硬套了一下，<br>clc,clear all<br>UH=[-6.22,-8.21,-10.23,-12.21,-14.21,-16.16];<br>d=0.5;<br>Is=3.00;<br>C=0.449;<br>Im=[0.300,0.400,0.500,0.600,0.700,0.800];<br>B=Im.*C;<br>RH=(UH.*d)./(Is.*B);<br>Pm=polyfit(B,RH,2)<br>RH=Pm(2)+Pm(1)*B;<br>RH<br>r=corrcoef(B,UH)<br>得到的结果是<br>r =</P>
<P>    1.0000   -1.0000<br>   -1.0000    1.0000<br>请问怎么回事？实际上不可能这么好的吧？应该略少于1的吧？而且不可能同时出现正负值的！</P>

[此贴子已经被作者于2005-12-20 15:09:54编辑过]

tammy

这就不是matlab问题了