matlab期末课程设计.doc

学习资料收集于网络，仅供参考情（亲情）（友情）（爱情）（热情）（情意）很清很清的河水 很大很大的西瓜 很长很长的路7、字的结构分析笔画笔顺填空。A B C D E F G H I J K L M N O P Q5、根据表达的需要，学会使用逗号、句号、问号和感叹号。（以后）（以前） （冷淡）（热情） 黑暗（光明、明亮）（ ）月（ ）日是元旦节。 （ ）月（ ）日 是中秋节。暖和的衣裳 暖和的被子 暖和的天气课件|教案|试卷|无需注册a b c d e f g h i j k l m n o p q中山大学南方学院电气与计算机工程学院课程名称：MATLAB语言及应用项目名称：期末课程设计姓名：叶健行学号：1820180102018 年 12 月 25 日学习资料目录1 课程设计目的12 课程设计题目及要求13 MATLAB程序及结果34 总结191 课程设计目的（1）掌握MATLAB基本知识及技能。（2）掌握数值分析方法及其常见的数值算法MATLAB编程实现的能力。（3）具备选择合适的开发平台进行MATLAB程序设计的能力。（4）具备解决问题，进行初步MATLAB程序设计的能力。（5）初步了解学习如何使用MATLAB进行计算、求根、绘图、拟合。（6）掌握数学建模问题的一般解决步骤。（7）掌握数据插值等问题解决方法。2 课程设计题目及要求选择第1、4、5、9、10、11、12共7题进行作答。1 给出某地一个月内温度的测定值（内含平均温度、最低、最高温度），并求其平均温度、最低、最高温度的平均值、标准差。（参考下表) （15分） 四月份每天的温度报表温度日期 平均温度（） 最高温度（）最低温度（）110.2198214.31910313.0199415.8191259.715768.0164715.2208814.8207910.31771012.11971114.72291215.327121322.629171420.727161522.331191613.015121716.925121813.51610198.7116207.31142111.81582210.91482315.621122417.825132518.624142622.729172723.528182824.529182918.627143018.927154 试依序回答下面的问题 （20分）(a) 用plot绘制函数f(x)=xcos(x)与g(x)=sin(x)的图形于同一个视窗内（Figure 1）；(5分)(b) 用subplot将两个函数图形绘制在同一视窗的左右两幅子图内(Figure 2)。绘图范围取。(5分)(c) 在Figure 1和Figure 2中用属性编辑区将f(x)图形改成红色，g(x)的图形改为紫色，线条粗细为3.0。(4分)(d) 在Figure 1加入图例说明，其中f(x)的图例为xcos(x)，g(x)的图例为sin(x)。(2分)(e) 在Figure 1加上图形的标题，标题名称为Function Plots，字体为Helvetica，大小为16。(2分) (f) 在Figure 1和Figure 2设定x轴的文字解为x，y轴的文字解说为f(x) & g(x)。(2分)5 试画出，的图形：（10分）(a) 分别用mesh函数和surf函数画出三维立体图，描述其区别。(5分)(b) 画出等高线图。(5分)9 弦割法求方程在区间-2，-1内的一个实根近似值，使.（15分）10 求的根。（10分）11用高斯算法编程求解下列线性方程组的解。（15分）12 在实验中测得如下10组数据: （20分）X 1 2 3 4 5 6 10 12 15 16 Y 7 15 19 30 38 37 23 69 39 28（1）求最多能拟合多项式的次数是多少？并求出各项系数。（2） 将数据点和拟合曲线在同一图中绘出。3 MATLAB程序及结果1 给出某地一个月内温度的测定值（内含平均温度、最低、最高温度），并求其平均温度、最低、最高温度的平均值、标准差。（参考下表) （15分）四月份每天的温度报表温度日期 平均温度（） 最高温度（）最低温度（）110.2198214.31910313.0199415.8191259.715768.0164715.2208814.8207910.31771012.11971114.72291215.327121322.629171420.727161522.331191613.015121716.925121813.51610198.7116207.31142111.81582210.91482315.621122417.825132518.624142622.729172723.528182824.529182918.627143018.92715代码：A=10.2198;14.31910;13.0199;15.81912;9.7157;8.0164;15.2208;14.8207;10.3177;12.1197;14.7229;15.32712;22.62917;20.72716;22.33119;13.01517.3114;11.8158;10.9148;15.62112;17.82513;18.62414;22.72917;23.52818;24.52918;18.62714;18.92715mean(A)std(A)2;16.92512;13.51610;8.7116;代码截图：结果：4 试依序回答下面的问题 （20分）(g) 用plot绘制函数f(x)=xcos(x)与g(x)=sin(x)的图形于同一个视窗内（Figure 1）；(5分)(h) 用subplot将两个函数图形绘制在同一视窗的左右两幅子图内(Figure 2)。绘图范围取。(5分)(i) 在Figure 1和Figure 2中用属性编辑区将f(x)图形改成红色，g(x)的图形改为紫色，线条粗细为3.0。(4分)(j) 在Figure 1加入图例说明，其中f(x)的图例为xcos(x)，g(x)的图例为sin(x)。(2分)(k) 在Figure 1加上图形的标题，标题名称为Function Plots，字体为Helvetica，大小为16。(2分) (l) 在Figure 1和Figure 2设定x轴的文字解为x，y轴的文字解说为f(x) & g(x)。(2分)(a)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);plot(x,y1,x,y2);代码截图：结果：Figure 1(b)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);subplot(121); plot(x,y1);subplot(122); plot(x,y2);代码截图：结果：Figure 2(c)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);plot(x,y1,r,x,y2,m);代码截图：结果：Figure 1x=0:0.1:18;y1=x.*cos(x);y2=sin(x);subplot(121); plot(x,y1,r);subplot(122); plot(x,y2,m);代码截图：结果：Figure 2(d)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);plot(x,y1,x,y2);legend(xcos(x),sinx);代码截图：结果：Figure 1(e)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);plot(x,y1,x,y2);legend(xcos(x),sinx);title(Function Plots,fontname,Helvetica,FontSize,16);代码截图：结果：Figure 1(f)代码：x=0:0.1:18;y1=x.*cos(x);y2=sin(x);subplot(121); plot(x,y1);xlabel(x)ylabel(f(x)subplot(122); plot(x,y2);xlabel(x)ylabel(g(x)代码截图：结果：5 试画出，的图形：（10分）(c) 分别用mesh函数和surf函数画出三维立体图，描述其区别。(5分)(d) 画出等高线图。(5分)(a)Mesh:代码：x=-10:0.5:10;y=-10:0.5:10;X,Y=meshgrid(-10:0.5:10);r=sqrt(X.2+Y.2);Z=cos(r);mesh(X,Y,Z)代码截图：结果：surf:代码:x=-10:0.5:10;y=-10:0.5:10;X,Y=meshgrid(-10:0.5:10);r=sqrt(X.2+Y.2);Z=cos(r);surf(X,Y,Z)代码截图：结果：区别：mesh命令绘制的图形是一个一排排的彩色曲线组成的网格图；而surf命令绘制得到的是着色的三维曲面。(b)代码：x=-10:0.5:10;y=-10:0.5:10;X,Y=meshgrid(-10:0.5:10);r=sqrt(X.2+Y.2);Z=cos(r);contour(X,Y,Z)代码截图：结果：9 弦割法求方程在区间-2，-1内的一个实根近似值，使.（15分）代码：function i,A=n(f,x0,x1,eps,NMAX)f=inline(x3-3*x2-x-9);x0=-2;x1=-1;NMAX=100;eps=0.00001;Tn=(f(x1)-f(x0)/(x1-x0);x2=x1-f(x1)/Tn;i=1;A=;A=A,x0,x1,x2;while (abs(x2-x1)eps & iNMAX) x0=x1; x1=x2; Tn=(f(x1)-f(x0)/(x1-x0); x2=x1-f(x1)/Tn; i=i+1; A=A,x2;endi,A代码截图：结果：10 求的根。（10分）代码：t=-5:1:5;x1=(sin(t).2;x2=-0.1.*t;y=x1.*exp(x2)-0.5.*abs(t);roots(y)代码截图：结果：11用高斯算法编程求解下列线性方程组的解。（15分）代码：function X=A1(A,b)%A 是一个 N 阶矩阵，b 是 n 维向量，X 是解A=1 5 9;2 6 10;3 7 11;4 8 12;b=13;14;15;16;N N=size(A);X=zeros(1,N+1);B=A b; %增广矩阵for p=1:N-1Y,j=max(abs(B(p:N,p); C=B(p,:); B(j+p-1,:)=C;if B(p,p)=0 A 是奇异阵，方程无唯一解 end for k=p+1:N m=B(k,p)/B(p,p); B(k,p:N+1)=B(k,p:N+1)-m*B(p,p:N+1) endend %程序函数 hui_dai 进行回代 X=A1_hui_dai(B(1:N,1:N),B(1+N,N+1); function X=A1_hui_dai(A,b) n=length(b); X=zeros(n,1); X(n)=b(n)/A(n,n); for k=n-1:-1:1X(k)=(b(k)-A(k,k+1:n)*X(k+1:n)/A(k,k); end代码截图：结果：12 在实验中测得如下10组数据: （20分）X 1 2 3 4 5 6 10 12 15 16 Y 7 15 19 30 38 37 23 69 39 28（1）求最多能拟合多项式的次数是多少？并求出各项系数。（3） 将数据点和拟合曲线在同一图中绘出。解析：因题目给出的实验数据是10组，可以看做是10个多项式方程组，仅当多项式的次数小于方程组的个数时，方程组为超定方程，此时才可以进行最小二乘多