matlab 程序 求助

小错误不少，调通了。

%Function [varargout]=xxfx1(x,y)

%方差分析表

x=[76.9,82.1,80.1,77.1,82.4,85.1,87.4,85.2,82.6,86.5,77.6,76.1,73.7,87.6,73.1,69.7,70.2,

68.7,64.1,86.0,104.5,90.7,89.7,67.5,66.1,68.8,69.7,72.5,73.6,73.2,71.6,80.7,68.5,64.2]

y=[65.0,74.2,69.8,76.1,75.3,77.4,76.7,77.4,76.7,75.5,67.2,66.2,64.1,76.1,68.4,67.2,64.9,

65.7,60.8,77.9,95.8,80.2,80.2,63.6,61.9,64.7,62.9,66.4,63.1,67.4,68.9,70.3,64.6,59.6]

alpha1=0.05

alpha2=0.01

x=x(:)

y=y(:)

n=length(y)

SST=sum((y.*y))-(sum(y)).*(sum(y))/n

lxx=sum(x.*x)-sum(x).*sum(x)/n

lxy=sum(x.*y)-sum(x).*sum(y)/n

SSR=lxy.*lxy/lxx

SSE=SST-SSR

df1=1

df2=n-2

df3=n-1

VR=SSR/df1

VE=SSE/df2

Fb=VR/VE

F=finv(1-[alpha1alpha2],df1,df2)

F1=min(F)

F2=max(F)

if Fb>F2

tst='**'

else if (Fb>F1)&(Fb<=F2)

tst='*'

else

tst=' '

end

end

format short g

table=cell(4,7)

table{1,1}='方差来源'

table{2,1}='回归'

table{3,1}='剩余'

table{4,1}='总和'

table{1,2}='偏平方和'

table{2,2}=SSR

table{3,2}=SSE

table{4,2}=SST

table{1,3}='自由度'

table{2,3}=df1

table{3,3}=df2

table{4,3}=df3

table{1,4}='方差'

table{2,4}=VR

table{3,4}=VE

table{1,5}='F比'

table{2,5}=Fb

table{1,6}='Fα'

table{2,6}=F1

table{3,6}=F2

table{1,7}='显著性'

table{2,7}=tst

if nargout>=1

varargout{1}=table

else

disp(table)

end

结果：

...

table =

'方差来源''偏平方和''自由度''方差' 'F比' 'Fα' '显著性'

'回归'[ 1699.9][ 1][1699.9][279.99][4.1491] []

'剩余'[ 194.28][32][6.0712] [][7.4993] []

'总和'[ 1894.1][33] [] [] [] []

table =

'方差来源''偏平方和''自由度''方差' 'F比' 'Fα' '显著性'

'回归'[ 1699.9][ 1][1699.9][279.99][4.1491]'**'

'剩余'[ 194.28][32][6.0712] [][7.4993] []

'总和'[ 1894.1][33] [] [] [] []

>>A = magic(4)

B = A.^2

SSR = sum(B(:))

MSR = A^2

A =

16 2 3 13

5 11 10 8

9 7 6 12

4 14 15 1

B =

256 4 9 169

25 121 100 64

81 49 36 144

16 196 225 1

SSR =

1496

MSR =

345 257 281 273

257 313 305 281

281 305 313 257

273 281 257 345

数值分析

主要分支致力于开发矩阵计算的有效算法，这是一个已持续几个世纪以来的课题，是一个不断扩大的研究领域。 矩阵分解方法简化了理论和实际的计算。 针对特定矩阵结构（如稀疏矩阵和近角矩阵）定制的算法在有限元方法和其他计算中加快了计算。 无限矩阵发生在行星理论和原子理论中。 无限矩阵的一个简单例子是代表一个函数的泰勒级数的导数算子的矩阵。

---