炼油厂将abc三种原油 matlab求解,matlab 上机作业1,中国石油大学

第 一 次 上 机 实 习

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

三条曲线