求助编个matlab程序

lhq213

哪位高手帮忙编个程序，matlab的Ｍ文件
有两组数据V,D，利用关系式V=a*D^b,求参数a,b的估计量.
很急啊！！！在线等！！！

branch

不知道你的数据维数？ab是矩阵？？
你可以参照数值计算的办法：

lhq213

数值计算的办法?
V,D都是100维的向量.是已知的数据。
拜托帮忙写个matlab的M文件，我要a,b的估计值，是不是还需要收敛性条件？我不会用matlab，但是现在有任务，必须要找到a,b的估计值。
现在真的很急！！！拜托！！！

branch

load m.mat    %将保存数据的文件调入；也可直接输入数据V=[]; D=[];
v=log(V);                %log(V)=log(a)+b*log(D);
n=size(D);
d=[log(D);ones(n,1)];
[loga,b]=v/d;
a=10^(loga);

[ 本帖最后由 branch 于 2007-4-30 15:58 编辑 ]

xjzuo

这个问题作个回归分析就可完成.
(拟合即得参数)

lhq213

y=[log(D);ones(n,1)]；
这句好像有问题，
Warning: Input arguments must be scalar.
??? Error using ==> vertcat
All rows in the bracketed expression must have the same
number of columns.
是什么原因啊？
我真是太菜了！！！惭愧！！！

zionward

我想在程序执行时在commond window里显示“please enter the value ***"     然后将这个值送到程序中。
而且要输入多个值    哪位高手能指教一下啊。在此谢了

branch

需要根据你输入的D矩阵的行列进行变换，我这里默认为行向量。

branch

a=input('请输入××的值：')

zjzza

提供数据，1分钟后给出参数估计结果