附1
31
附2 直流测功机特性方程
一、方程式
T2?CT?IF?T0?a1IF?a2其中，Te?CT?IF 二、最小二乘原理
Jmin???(k)?ξTξ?(z?Φθ)T(z?Φθ)?z?Φθ
k?1N22daTbdaTdbT?b?a dθdθdθ??[ΦTΦ]?1ΦTz θ三、编程
If=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0]; T2=[0.182 0.290 0.395 0.505 0.612 0.720 0.828 0.936 1.045 1.155 1.265 1.375 1.482 1.588 1.690 1.789 1.887 1.986 2.08 2.170 2.251]; Y=T2';
one=ones(length(If),1); X=[If' one]; Theta=inv(X'*X)*X'*Y plot(x,y) hold on plot(If,T2,'r*')
结果 Theta = 1.0513 0.1978
四、估计曲线与实验曲线比较
2.521.510.5000.20.40.60.811.21.41.61.82
32