MATLAB求三元回归方程代码

用Matlab的solve()函数求解.例如eq1='x+(y+z)/2=w';eq2='y+x+w-10=z';eq3='x+w=z+y/4';eq4='w+z=y+8*x-1';[X1,X2,X3

A=[13362;1111]';b=[22526]';c=A\b%y=c(1)*x+c(2)

你把太多无关的部分混进来了.从代码上看,S=f(A1,a,A2,b,A3,t,c)而A1=f(k),k已知；a=f(t),已知；……其他的也都已知,根本不需要方程啊,直接求出来了.

试一下这个：symsxyz[x,y,z]=meshgrid(-2:0.2:2,-2:0.25:2,-2:0.16:2);v=((50.*(1+0.0001.*x).^0.5)./(0.0484.*pi

figure;t=[1999200020012002200320042005200620072008];c=[399.72506.97754.98989.41202.481473.2915251717

wavread函数读取Microsoft的wav格式音频文件,返回值x是音频的数据向量samples{i}{j}=wavread(fname);读取\\train\\dgt（i-1）\\dgt（j-1

t=sym('t','real');q=1500;m=mod(t,30);n=(t-m)/30;u=0.00.07095t^3-6.443t^2+92.19t+734.2;y=mod(n+1,2)*(

程序是没错的,但显然拟合效果很不理想.y=[354355474357262728292229191114232022138327265];x1=[222222555555777777101010101

用regress(y,x)；其中y是nx1阶x是n乘p阶向量能得到系数.用cftool(x,y)在提供的界面直接拟合!.还有多项式拟合polyfit（x,y,n);n是多项式的阶数.

A = [1,2,3;2,3,4;4,5,5];b = [1;3;2];X = A\b;这样得到的X是一个列向量,正好对应于[a0;a1;a

2B再问：你个2B

可以不用拟合工具箱,直接用矩阵除法即可!因为为线性求a1,a2即把a1,a2当成未知数,x1,x2,Y-a0当成已知量则x1*a1+x2*a2=Y-a0,即[x1,x2]*[a1;a2]=Y-a0令矩

实话说,没心思看这样的代码,全是计算.估计mx和my是某点的坐标（我才是圆心）R应该就是半径mse就不知道了.再问：。。。就是想知道mse是什么。。。。。。。。。帮忙看看行吗。。给你追加。。。今天早上

x=[2.82,2.02,1.57,1.28,1.08,0.94];z=[40,55.9,62.8,66.2,68.7,70.6];Z=z';X1=[x;x]';%c=ax-bx+zz=c-ax+bx

如果用xz-x+>0与xz-y-1>0约束一下,第一个约束不清晰xz-x+>0====》xz-x+1>0functionaay=@(x)-(x(1)-1).*(1+x(2)-x(1).*x(3))./

A=[1,2,3,4,5;3,4,5,6,7]A=1234534567>>B=(A-min(A(:)))./(max(A(:))-min(A(:)))B=00.16670.33330.50000.66

请用regress函数求解x1=[0.61340.61220.66500.65070.6434];x2=[0.25070.02040.02870.01430.0286];x3=[0.05670.051

参考代码：>> syms w kc real>> Tjw = 7.8*kc/(7.8*kc-0.11*w^

可运行代码见附件.再问：大神，为什么你使用clearall;closeall，为什么大多数此类代码都是x=[ones(34,1),x1,x2];，而你给我发的代码确是x=[ones(1,34);x1;

而广告费与销售量成二次多项式（抛物线关系）关系LinearmodelPoly2:f(x)=p1*x^2+p2*x+p3Coefficients(with95%confidencebounds):p1=