第 一 次 上 机 实 习
1. 编写一个程序ex1_1.m,对输入的方阵计算其伴随矩阵。
a) 要求第一行帮助行注明这个函数的用途,能在matlab命令窗口可以用help或lookfor查到该帮助;
b) 对输入的变量进行检测,如果输入的不是方阵或者矩阵元素不是数,则报错,退出程序。
function A=ex1_1(X)
% 建立函数A=ex1_1(X),求矩阵X的伴随矩阵
[a,b]=size(X);
if a~=b||~isnumeric(X) || ~ismatrix(X)
error('输入的不是方阵或者矩阵元素不是数!');
end
for i=1:a
for j=1:b
A(i,j)=(-1)^(i+j)*X(i,j)*det(X([1:i-1,i+1:end],[1:j-1,j+1:end]));
end
end
A=A';
2. 已知表达式y=6*x^5+4*x^3+2*x^2-7x+10,x的范围是[0,100],使用三阶拟合和五阶的方法得出多项式的表达式,并编程在图中绘制出原曲线、三阶拟合和五段拟合的曲线。 clc
clear all
x=linspace(0,100,2000);
p1=[6 0 4 2 -7 0];
y=polyval(p1,x);
p3=polyfit(x,y,3);
p5=polyfit(x,y,5);
y1=polyval(p3,x);
y2=polyval(p5,x);
plot(x,y,'r-o')
hold on
grid on
plot(x,y1,'b-.',x,y2,'y--*');
grid on
三条曲线