Matlab高等代数一..docVIP

《Matlab→高等代数》 实验一 学号： 姓名： 年级专业班级： 实验时间： 年 月 日 时— 时 实验地点： 授课教师：谢乐平 实验一 Matlab与高等代数，多项式及行列式 一、实验目的 1、了解Matlab在高等代数中应用的基本操作； 2、掌握Matlab应用于多项式的运算，带余除法，求最大公因式，求根及因式分解。 3、掌握Matlab应用于计算行列式及cramer法则求解线性方程组。 二、实验内容 一、基本操作 注意: 后面是要输入电脑的内容。%后面是表示注释内容。 ★基本数学运算，只需将运算式直接打入提示号（）之後，并按Enter键即可。例如： (5*2+1.3-0.8)*10/25 %显示：ans =4.2000 x = (5*2+1.3-0.8)*10^2/25 %显示：x = 42 若不想让MATLAB显示运算结果，只需在运算式後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2); %若要显示变数y的值，直接键入y即可 ★要进行符号运算须先要对符号变量进行定义，一般有下面的方法进行定义： (1)直接用引号‘ ’创建（包括代数式）。 (2)用函数sym(‘ ’)创建（包括代数）。 (3)用命令syms 创建多个符号变量，变量之间用空格分开。 注意：‘ ’中包括的空格等都视为符号。并且Matlab严格区分大小写。 符号函数赋值：subs( )；解 dsolve（）见下例： clear %表示清除内存中的变量 syms x a b %定义三个变量x a b f=a*x^3-x+b f1=subs(f,2) %f中x赋值2 f2=subs(f,a,2) %f中a赋值2 a=2; b=1; g=subs(f) t=1:4; g1=subs(g,t) plot(t,g1) %画图函数plot ★输入语句 x=input(please input a number:) %显示：please input a number:要求输入数给x x=input(please input a string:,s) %显示：please input a string: ★输出语句 自由格式 (disp) disp(23+454-29*4) disp([11 22 33; 44 55 66; 77 88 99]) disp(this is a string) 格式化输出 (fprintf)。 fprintf(The area is %8.5f\n, 12.345678) %输出格式前有%，跳行符有\符号。 ★如果对某个函数不明白，输入‘help 函数名’获得帮助，如 help fprintf 二、多项式 一、多项式的创建 clear %方法一(由多项式的系数确定多项式) p=[1 -2 3] %直接给出多项式的系数，按次数从高到低排列 poly2sym(p) %给出以p中数为系数的多项式的表达式 %方法二(先定义变量，再定义多项式) syms x %定义变量x f=3*x^4-4*x^2+8 %给出多项式 poly2sym(f)%给出f多项式的表达式 %方法三（由多项式的根确定的多项式） x=[-1 2] %给出多项式的根 px=poly(x)%以x的元素为根的多项式。 poly2sym(px)%给出ps多项式的表达式 ★引用多项式函数求多项式的函数值 Y=polyval(px,x) %引用函数为polyval。px为多项式，x自变量取值，Y对应的函数值。 clear d=[-1 2] px=poly(d) y=polyval(px,4)%求多项式px在自变量等于4时的函数值 x=-4:0.5:8 yx=polyval(px,x)%求多项式px在自变量等于x序列时的函数值序列 plot(x,yx)%作出两个变量的函数图 a=roots(px)%求多项式px的根 二、多项式的相关函数 乘:conv；除：deconv；微分：polyder；求根：roots或solve；取得系数：sym2poly； 带余除法：deconv；求最大公因式：gcd；符号表达式的展开：expand；符号因式分解：factor。 三、多项式运算及带余除法 ◆讨论. clear a1=[1 0] %表示多项式x a2=[2 3] %表示多项式2x+3 a3=[1 18] %表示多项式x+18 p1=conv(a1,a2) %表示a1与a2相乘 p2=conv(p1