exn554.m

来自「基天子matlab的双缝干涉从单缝(线光源)发出的单色光射到与之平行的双缝上」· M 代码 · 共 24 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-5-4程序
% 曲线的插值和拟合
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月

clear,clf,format compact
C=[1,0,0,0;1,1,1,1;1,2,4,8;1,3,9,27], 
b=[3;0;-1;6],U=rref([C,b])
pause
t1=1.5,f1=3-2*t1-2*t1^2+t1^3
pause
ezplot('t^3-2*t^2-2*t+3',[-1,4]),
grid on, hold on
plot(0:3,[3,0,-1,6],'*')
plot(1.5,-1.125,'or')
pause
A1=[1,0,0,0;1,1,1,1;1,2,4,8;1,3,9,27;1,-1,1,-1], 
b1=[3;0;-1;6;5], U01=rref([A1,b1])
pause
a1=pinv(A1)*b1,
E=A1*a1-b1,EE=norm(E)
t=-1:0.1:4;
plot(t,a1(1)+a1(2)*t+a1(3)*t.^2+a1(4)*t.^3,':r')
plot(-1,5,'x')