Matlab的多項式运算和数据分析.doc

课程名称： Matlab语言程序设计 实验题目： Matlab的多项式运算和数据分析 一、实验目的 1. 掌握Matlab中多项式函数的使用； 2. 掌握利用Matlab实现数据拟合的方法； 3. 掌握Matlab中的数据分析函数的使用。 二、实验内容和结果 1.求多项式函数的根，求以得到的根向量为根的多项式；求的导函数系数向量和以5为常数项的原函数系数向量；在同一绘图窗口中用plot绘制,,的图像，并加图例。 clc p=[2 3 0 4]; x=0:0.01:2; p1=roots(p); yy=poly2sym(p1) y=polyval(p,x); plot(x,y,g) hold on q1=polyder(p) y1=polyval(q1,x); plot(x,y1,b) hold on q2=[2 3 0 5]; y2=polyval(q2,x); plot(x,y2,r) hold on legend(p(x),q1(x),q2(x)) xlabel(x) ylabel(y) 2.用电压伏的电池给电容器充电，电容器上时刻的电压为，其中是电容器的初始电压，是充电常数。试由下面一组，数据确定和。 （秒） 0.5 1 2 3 4 5 7 9 （伏） 6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63 t=[.5 1 2 3 4 5 7 9]; v=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63]; u1=10-v; u3=log(u1); u4=exp(u3); plot(t,u3,*r); hold on q=polyfit(t,u3,1) y=polyval(q,t); plot(t,y); tt=1/-.2835 v0=10-exp(1.4766) v0 = 5.6220 =-3.5273 =5.6220 拟合图 3.利用多项式函数和eig求3阶Pascal矩阵的特征值，并比较二者的结果。 clc A=pascal(3); z1=eig(A); p=poly(A); z2=roots(p); isequal(sort(z1),sort(z2)) ans = 0 4.输入向量，按从大到小的顺序排序。 x=[1 2 5 6 81 4]; fliplr(sort(x)) 5.求数组中所有绝对值大于3的元素的个数。 A=[-4 -2 0 2 4;-3 -1 1 3 5]; B=(A3)|(A-3); sum(sum(B)) ans = 3 三、思考及总结 上课时要认真听讲，上机实验时要多操作，不会的要多问多练习，不懂的多问，多查书和课件。多思考，多上机实践。  - 4 -