Matlab实验报告
基于matlab的实验报告
基于matlab的实验报告实验报告：基于MATLAB 的实验一、实验目的通过使用MATLAB 软件，掌握如何进行数据分析、图像处理、算法实现等一系列实验操作，提高实验者的实践能力和动手能力。
二、实验原理MATLAB 是一种在科学计算和技术开发领域广泛应用的计算机软件。
它能进行矩阵计算、绘制函数和数据图像、实现算法以及进行数据分析等。
通过掌握MATLAB 的使用，能够快速、高效地解决各种科学和工程问题。
三、实验内容1. 数据分析：使用MATLAB 的数据分析工具进行数据的导入、处理和分析。
2. 图像处理：利用MATLAB 的图像处理工具包对图像进行滤波、增强、分割等操作。
3. 算法实现：使用MATLAB 实现常用的算法，如排序、搜索、图像压缩等。
四、实验步骤1. 数据分析：（1）使用MATLAB 的读取数据函数将数据导入MATLAB 环境中。
（2）利用MATLAB 的数据处理函数进行数据清洗和预处理。
（3）使用MATLAB 的统计工具进行数据分析，如求平均值、标准差等。
（4）利用MATLAB 的绘图函数将分析结果可视化。
2. 图像处理：（1）使用MATLAB 的读取图像函数将图像导入MATLAB 环境中。
（2）利用MATLAB 的图像处理工具包进行滤波操作，如均值滤波、中值滤波等。
（3）使用MATLAB 的图像增强函数对图像进行锐化、变换等操作。
（4）利用MATLAB 的图像分割算法对图像进行分割。
3. 算法实现：（1）使用MATLAB 编写排序算法，如冒泡排序、快速排序等。
（2）使用MATLAB 编写搜索算法，如二分查找、线性搜索等。
（3）使用MATLAB 实现图像压缩算法，如离散余弦变换（DCT）。
五、实验结果实验中，我们使用MATLAB 完成了数据分析、图像处理和算法实现的一系列实验操作。
通过数据分析，我们成功导入了数据并对其进行了清洗和预处理，最后得到了数据的统计结果。
在图像处理方面，我们对图像进行了滤波、增强和分割等操作，最终得到了处理后的图像。
程序设计实验报告(matlab)
程序设计实验报告(matlab)实验一: 程序设计基础实验目的：初步掌握机器人编程语言Matlab。
实验内容：运用Matlab进行简单的程序设计。
实验方法：基于Matlab环境下的简单程序设计。
实验结果：成功掌握简单的程序设计和Matlab基本编程语法。
实验二：多项式拟合与插值实验目的：学习多项式拟合和插值的方法，并能进行相关计算。
实验内容：在Matlab环境下进行多项式拟合和插值的计算。
实验方法：结合Matlab的插值工具箱，进行相关的计算。
实验结果：深入理解多项式拟合和插值的实现原理，成功掌握Matlab的插值工具箱。
实验三：最小二乘法实验目的：了解最小二乘法的基本原理和算法，并能够通过Matlab进行计算。
实验内容：利用Matlab进行最小二乘法计算。
实验方法：基于Matlab的线性代数计算库，进行最小二乘法的计算。
实验结果：成功掌握最小二乘法的计算方法，并了解其在实际应用中的作用。
实验六：常微分方程实验目的：了解ODE的基本概念和解法，并通过Matlab进行计算。
实验内容：利用Matlab求解ODE的一阶微分方程组、变系数ODE、高阶ODE等问题。
实验方法：基于Matlab的ODE工具箱，进行ODE求解。
实验结果：深入理解ODE的基本概念和解法，掌握多种ODE求解方法，熟练掌握Matlab的ODE求解工具箱的使用方法。
总结在Matlab环境下进行程序设计实验，使我对Matlab有了更深刻的认识和了解，也使我对计算机科学在实践中的应用有了更加深入的了解。
通过这些实验的学习，我能够灵活应用Matlab进行各种计算和数值分析，同时也能够深入理解相关的数学原理和算法。
这些知识和技能对我未来的学习和工作都将有着重要的帮助。
MATLAB实验报告
v1.0 可编辑可修改实验一 MATLAB 环境的熟悉与基本运算一、实验目的及要求1．熟悉MATLAB 的开发环境； 2．掌握MATLAB 的一些常用命令；3．掌握矩阵、变量、表达式的输入方法及各种基本运算。
二、实验内容1.熟悉MATLAB 的开发环境： ① MATLAB 的各种窗口：命令窗口、命令历史窗口、工作空间窗口、当前路径窗口。
②路径的设置：建立自己的文件夹，加入到MATLAB 路径中，并保存。
设置当前路径，以方便文件管理。
2.学习使用clc 、clear ，了解其功能和作用。
3.矩阵运算：已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A.*B ，并比较结果。
4.使用冒号选出指定元素:已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素；A 中所有列第2，3行的元素； 5.在MATLAB 的命令窗口计算: 1) )2sin(π2) 5.4)4.05589(÷⨯+ 6.关系及逻辑运算1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15]，求: y=a==b ，并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0]，求: x&y+x>y ，并分析结果 7.文件操作1)将0到1000的所有整数，写入到D 盘下的文件 2)读入D 盘下的文件，并赋给变量num8.符号运算1)对表达式f=x 3-1 进行因式分解2)对表达式f=（2x 2*(x+3)-10)*t ，分别将自变量x 和t 的同类项合并 3）求3(1)xdz z +⎰三、实验报告要求完成实验内容的3、4、5、6、7、8，写出相应的程序、结果实验二 MATLAB 语言的程序设计一、实验目的1、熟悉 MATLAB 程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计方法3、函数文件的编写和设计4、了解和熟悉变量传递和赋值二、实验内容1．编写程序，计算1+3+5+7+…+(2n+1)的值（用input 语句输入n 值）。
matlab实验报告
matlab实验报告实验1 熟悉matlab 的开发环境及矩阵操作⼀、实验的教学⽬标通过本次实验使学⽣熟悉MATLAB7.0的开发环境，熟悉MA TLAB ⼯作界⾯的多个常⽤窗⼝包括命令窗⼝、历史命令窗⼝、当前⼯作⽬录窗⼝、⼯作空间浏览器窗⼝等。
掌握建⽴表达式书写规则及常⽤函数的使⽤，建⽴矩阵的⼏种⽅法。
⼆、实验环境计算机、MATLAB7.0集成环境三、实验内容1、熟悉命令窗⼝的使⽤，⼯作空间窗⼝的使⽤，⼯作⽬录、搜索路径的设置。
命令历史记录窗⼝的使⽤，帮助系统的使⽤。
2、在当前命令窗⼝中输⼊以下命令：x=0:2:10 y=sqrt(x)，并理解其含义。
3、求下列表达式的值（1）w=)1034245.01(26-?+?（2）x=ac b e abc cb a ++-+++)tan(22ππ,其中a=3.5,b=5,c=-9.8 四、实验总结1、熟悉了命令窗⼝的使⽤，⼯作空间窗⼝的使⽤。
2、了解了⼯作⽬录、搜索路径的设置⽅法。
---5317383399351542实验2 MATLAB 基本运算⼀、实验的教学⽬标通过本次实验使学⽣掌握向量和矩阵的创建⽅法；掌握矩阵和数组的算术运算、逻辑运算和关系运算；掌握字符数组的创建和运算；了解创建元胞数组和结构体的⽅法。
⼆、实验环境计算机、MATLAB7.0集成环境三、实验内容1、要求在闭区间]2,0[π上产⽣具有10个等距采样点的⼀维数组。
试⽤两种不同的指令实现。
（提⽰：冒号⽣成法，定点⽣成法）2、由指令rng('default'),A=rand(3,5)⽣成⼆维数组A ，试求该数组中所有⼤于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。
（提⽰：find 和sub2ind ）3、创建3阶魔⽅矩阵a 和3阶对⾓阵b ，c=a(1:3,1:3)（1）计算矩阵a,b 和c 的⾏列式、逆矩阵并进⾏最⼤值的统计。
（2）⽐较矩阵和数组的算术运算：b 和c 的*、/、^和.*、./、.^。
matlab实验报告总结
matlab实验报告总结1.求一份matlab的试验报告计算方法试验报告3【实验目的】检查各种数值计算方法的长期行为【内容】给定方程组x'(t)=ay(t),y'(t)=bx(t), x(0)=0, y(0)=b的解是x-y 平面上的一个椭圆，利用你已经知道的算法，取足够小的步长，计算上述方程的轨道，看看那种算法能够保持椭圆轨道不变。
（计算的时间步长要足够多）【实验设计】用一下四种方法来计算：1. Euler法2. 梯形法3. 4阶RK法4. 多步法Adams公式【实验过程】1. Euler法具体的代码如下：clear;a=2;b=1;A=[0 a; -b0];U=[];u(:,1)=[0;b];n=1000000;h=6*pi/n;fori=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5; u(:,i+1)=u(:,i)+h*A*u(:,i);endt=1:n+1;subplot(1, 2,1);plot(1:n,delta);gridon;subplot(1,2,2);plot(u(1,:),u(2,:));gridon;max(abs(delta-ones(1,length(delta)）））；结果如下：2. 梯形法具体的代码如下：clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=300;h=6*pi/n;for i=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;v1=u(:,i)+h*A*u(:,i);v2=u(:,i)+h*A*(u(:,i)+v1)/2;1u(:,i+1)=u(:,i)+h*A*(u(:,i)+v2)/2;endt=1:n+1;sub plot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下 3. 4阶RK法clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=70;h=6*pi/n;for i=1:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;k1=A*u(:,i); k2=A*(u(:,i)+h/2*k2); k3=A*(u(:,i)+h*k3); k4=A*(u(:,i)+h*k3); u(:,i+1)=u(:,i)+h/6*(k1+2*k2+2*k3+k4);endt=1:n+1 ;subplot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下：4. 多步法Adams公式clear;a=2;b=1;A=[0 a; -b 0];U=[];u(:,1)=[0;b];n=200;h=6*pi/n;u(:;2)=u(u,1)+h*A*u(:,1);u(:;3)=u(u,2)+h/2*A*(3*u(:,2)-u(:,1));u(:;4)=u(u,3)+h/12*A*(23*u(:,3)-16*u(:,2)+5*u(:, 1)); delta(1)=((u(1,1)/a)^2+(u(2,1)/b^2)^0.5 delta(2)=((u(1,2)/a)^2+(u(2,2)/b^2)^0.5delta(3)=((u(1,3)/a)^2+(u(2,3)/b^2)^0.5for i=4:n delta(i)=((u(1,i)/a)^2+(u(2,i)/b)^2)^0.5;u(:,i+1)=u(:,i)+h/24*A*(55*u(:,i)-59*u(:,i-1)+37 *u(:,i-1)+37*u(:,i-2)-9*u(:,i-3));endt=1:n+1;sub plot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2)；结果如下：【实验分析】通过这几种方法对比，发现最为稳定的是多步法Adams公式和4阶RK法，其次是梯形法，而欧拉法最为不稳定。
MATLAB综合性实验报告8
一、实验目的：1.学会运用Matlab解决实际问题2.进一步掌握Matlab的一些基本操作3.通过范例体会Matlab的初步建模过程4学会用做动画二．实验仪器、设备或软件：电脑、MATLAB软件三．实验内容：1.已知前两年的猪肉的产量和猪肉的价格分别为：39吨，28吨，12元/公斤，17元/公斤，根据前一年的猪肉价格影响后一年猪肉产量当年猪肉产量影响当年猪肉价格的线性关系，编写程序，利用动画原理画出前十年猪肉的产量——价格的动态图形。
（参数设置为，c1=39,c2=28,c3=36,r1=12,r2=17,k=16)。
2．在地面上建有一座圆柱形水塔,水塔内部的直径为d,并且在地面处开了一个高为H的小门.现在要对水塔内部进行维修施工,施工方案要求把一根长为l(l>d)的水管运到水塔内部.请问水塔的门高H多高时,才有可能成功地把水管搬进水塔内。
四．实验记录：1．蛛网模型：实验代码：clear%c1为产量1, c2为产量2, c3为产量3, r1为%肉价1, r2为肉价2, k为K年后产量与肉价%是否稳定c1=39;c2=28;c3=36;hold offr1=12;r2=17;k=16;a1=[c1 1;c2 1];%系数矩阵b1=[r1,r2]';%列向量a2=[r1 1;r2 1];b2=[c2,c3]';a=a1\b1;b=a2\b2;%x0(1)=c1;for n=1:30y0(n)=a(1)*x0(n)+a(2); %a(1)a(2)为矩阵元素的引用，下行类似x0(n+1)=b(1)*y0(n)+b(2);x(n)=x0(n);y(n)=x0(n+1);endplot(x,y0,'-g',y,y0,'-b')hold onm=moviein(100);for n=1:kfor j=1:30t1=x0(n)+(j-1)*(x0(n+1)-x0(n))/30;t2=x0(n)+j*(x0(n+1)-x0(n))/30;if t2<t1t=t1;t1=t2;t2=t;%这样做方便画图elseendt=t1:0.01:t2;plot(t,y0(n),'.r')%划横线t1=y0(n)+(j-1)*(y0(n+1)-y0(n))/30;t2=y0(n)+j*(y0(n+1)-y0(n))/30;if t2<t1t=t1;t1=t2;t2=t;elseendt=t1:0.01:t2;plot(x(n+1),t,'.r')%划竖线endm(:,n)=getframe;endmovie(m,20)实验结果：ans =Columns 1 through 639.0000 28.0000 36.0000 30.1818 34.4132 31.3358 Columns 7 through 1233.5739 31.9462 33.1300 32.2691 32.8952 32.4398 Columns 13 through 1832.7710 32.5302 32.7053 32.5779 32.6706 32.6032 Columns 19 through 2432.6522 32.6166 32.6425 32.6236 32.6374 32.6274 Columns 25 through 3032.6346 32.6294 32.6332 32.6304 32.6324 32.6310 Column 3132.6320实验图像：2．水塔模型：实验代码：Clearfunction x=lt523(l,d,h) k1=0;a=0l=7;d=3ymax=0;k2=0;b=pi/2;h=3;while (b-a)>10^(-2)k2=k2+1;m=0;a=0;if k1==0n=ceil(b/h)-1 elsen=ceil(b/h);endfor i=1:n+1;x(i)=a+(i-1)*(b-a)/n;y(i)=l*sin(x(i))-d*tan(x(i));endfor i=1:nif y(i)>ymaxymax=y(i);a=x(i-1);else y(i)<ymaxb=x(i);endendm=m+1;k1=k1+1;h=(b-a)/5;endabx=(a+b)/2t=0:0.01:1.35y1=l*sin(t)-d*tan(t);plot(t,y1,'-b')方法：2syms h xh=8*sin(x)-2*tan(x);d1=0;d2=pi/2;d0=(d1+d2)/2;d=0.03;a=subs(h,x,d1);b=subs(h,x,d2);c=subs(h,x,d0);for n=d1:d:d2if (((c-a)/(d0-d1))*((b-c)/(d2-d0)))<0 d1=d1+d; d2=d2-d; d0=(d1+d2)/2;a=subs(h,x,d1);b=subs(h,x,d2);c=subs(h,x,d0);endendd0hold onaxis([0,pi/2,-5,5])x=0:0.05:pi/2;plot(x,8*sin(x)-2*tan(x),'-b')实验运行结果：五、实验总结：通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。
试验一熟悉matlab环境实验报告
试验一熟悉matlab环境(实验报告) 实验报告：熟悉MATLAB环境一、实验目的1.掌握MATLAB的基本操作方法，熟悉MATLAB的界面和工具栏。
2.了解MATLAB的变量类型、数据结构、矩阵运算等基本概念。
3.学会使用MATLAB进行数值计算、绘图和数据分析。
二、实验内容1.MATLAB界面与工具栏：熟悉MATLAB的启动界面，了解MATLAB的菜单栏、工具栏和命令窗口。
掌握如何创建、保存和打开MATLAB脚本文件。
2.变量与数据类型：学习MATLAB的变量命名规则，了解变量类型（如标量、向量、矩阵等）及其在内存中的存储方式。
掌握如何声明变量并赋值，理解MATLAB中的数据结构（如结构体、细胞数组等）。
3.矩阵运算：学习MATLAB的矩阵运算规则，了解矩阵的基本操作（如转置、逆等）。
掌握矩阵运算的常用函数（如求行列式、矩阵的秩等）。
4.数值计算与绘图：学习使用MATLAB进行数值计算和绘图。
掌握MATLAB的数值计算函数（如sin、cos、exp等）及其应用，了解如何绘制二维和三维图形。
5.数据分析：学习使用MATLAB进行数据分析。
了解如何读取和导出数据文件，掌握基本的统计方法和数据分析技术（如求平均值、中位数、标准差等）。
三、实验步骤与操作方法1.打开MATLAB软件，进入命令窗口。
2.创建并保存脚本文件：使用“edit”命令创建一个新的脚本文件，保存文件并命名（如“熟悉MATLAB环境.m”）。
3.变量与数据类型操作：在脚本文件中输入不同类型的变量并赋值，观察变量的输出结果。
例如，声明一个标量变量a和一个向量变量b，然后分别对它们进行赋值和输出。
4.矩阵运算：在脚本文件中输入一个矩阵并进行基本操作。
例如，声明一个3x3的矩阵A，对其进行转置、逆等操作，并输出结果。
5.数值计算与绘图：在脚本文件中使用数值计算函数并进行绘图操作。
例如，创建一个包含0到2π范围内100个点的向量x，计算sin(x)和cos(x)的值，绘制它们的图形并保存结果。
matlab实验报告总结精选
matlab实验报告总结电气工程学院自动化102班 2012年12月21日实验一 MATLAB环境的熟悉与基本运算一、实验目的1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识1.熟悉MATLAB环境MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令变量与运算符变量命名规则如下：变量名可以由英语字母、数字和下划线组成变量名应以英文字母开头长度不大于31个区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式的基本运算表7 两种运算指令形式和实质内涵的异同表的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、新建一个文件夹2、启动，将该文件夹添加到MATLAB路径管理器中。
3、保存，关闭对话框4、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye5、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
6、初步程序的编写练习，新建M-file，保存，学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
注意：每一次M-file的修改后，都要存盘。
练习A：help rand，然后随机生成一个2×6的数组，观察command window、command history和workspace等窗口的变化结果。
学习使用clc、clear，了解其功能和作用。
答：clc是清除命令窗体内容 clear是清除工作区间输入C=1:2:20，则C表示什么？其中i=1,2,3,?,10。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一：Matlab操作环境熟悉
一、实验目的
1．初步了解Matlab操作环境。
2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容
熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：
1．单函数运算操作。
求下列函数的符号导数
(1)y=sin(x); (2) y=(1+x)^3*(2-x);
求下列函数的符号积分
(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1-x^2);(4)y=(x1)/(x+1)/(x+2)
求反函数
(1)y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2));
代数式的化简
(1)(x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x-3*cos(
x)+4*x*sin(x);
2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化
(1)y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6)
y6=x^2/2
3．两个函数之间的操作
求和
(1)sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5
乘积
(1)exp(-x)*sin(x) (2) sin(x)*x
商
(1)sin(x)/cos(x); (2) x/(1+x^2); (3) 1/(x-1)/(x-2);
求复合函数
(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)
(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x
实验二：MATLAB基本操作与用法
一、实验目的
1.掌握用MATLAB命令窗口进行简单数学运算。
2.掌握常用的操作命令和快捷键。
3.了解MATLAB的数据类型。
4.了解MATLAB的操作符。
二、实验内容
1. 在命令窗口依次输入下列命令，根据执行结果分析其功能；
help
lookfor magic
help magic
doc magic
which magic。
2. 以自己姓名拼音来建立自己的工作目录，再讲自己的工作目录设置到
Matlab搜索路径下。
用cd命令查询自己的工作目录；
3. 创建变量，并计算：
1)创建single类型变量a=98，b=168，求：a+b，a-b，b×a3，b×a×a×a
2)创建uint8类型变量m，n，值与（1）同，进行相同计算。
记录命令和结果，并解释（1）与（2）结果为何不同
因uint8类型变量的最大值为255,当超过最大值时结果就等于255
4.先求下列表达式的值，然后记录Matlab工作空间的使用情况和相关变量
Z 3=2
3.0ln )3.0sin(23.3.0+++--a a e e a o a ，a=-3.0,-2.9,-2.8,……，2.8，2.9，3.0 实验三：矩阵运算与元素群运算
一、实验目的
1．掌握数组与矩阵的创建。
2．掌握矩阵运算与数组运算。
3．掌握基本元素群运算。
4．掌握向量与矩阵的特殊处理。
二、实验内容
1.”:”号的用法.用”:”号生成行向量a=[1 2 3 4 5 6 7 8 9 10]、b=[5 3 1 -1 -3 -5];
2．用线性等分命令linspace 重新生成上述的a 和b 向量。
3．在100和10000之间用对数等分命令logspace 生成10维的向量c 。
4．生成范围在[0，10]、均值为5的3×5维的均匀分布随机数矩阵D. 5．利用magic函数生成5×5维的魔方矩阵，取其对角向量e，并根据向量e 生成一个对角矩阵E。
（所谓魔方矩阵就是各行、各列、各对角线元素之和相等。
）
6．另AA是3×3维魔方矩阵，BB是由A旋转180°得到。
CC是一个复数矩阵，其实部为AA，虚部为BB。
DD是CC的转置，EE是CC的共轭。
分别计算CC和EE的模和幅角。
7.f是一个首项为20，公比为0.5的10维等比数列；g是一个首项为1，差为
3的10维等差数列。
试计算向量f和g的内积s。
8．生成一个9×9维的魔方矩阵，提取其中心的3×3维子矩阵M，利用sum 函数检验其各行和各列的和是否相等。
9．已知T=[12
23
34
45
34
45
56
67
],利用函数生成左上三角举证T1=[
12
23
34
40
34
40
00
00
].
*不存在2X2魔方矩阵*
实验四：函数编写与程序设计
一、实验目的
1．掌握函数的编写规则。
2．掌握函数的调用。
3．会用Matlab程序设计实现一些工程算法问题。
二、实验内容
1．编写一个[y,y1,y2]=mwave(f1,m1,f2,m2)函数，实现以下功能，并绘出y1、y2、y在t∈[0，2π]区间500个样点的图形。
（其中调用参数2 ≤ f1、f2 ≤
20 Hz；0.5 ≤ m1、m2 ≤ 2）
y1=m1sin(2πf1t);
y2=m2sin(2πf2t);
y=y1+y2
当输入的参数超过范围时
输入正确时的图形
2．程序设计：相传古代印度国王要褒奖他的聪明能干的宰相达依尔（国际象棋发明者），问他要什么？达依尔回答：“陛下只要在国际象棋棋盘的
第一个格子上放一粒麦子，第二个格子上放二粒麦子，以后每个格子的
麦子数都按前一格的两倍计算。
如果陛下按此法给我64格的麦子，就感激不尽，其他什么也不要了。
”国王想：“这还不容易！”让人扛了一
袋麦子，但很快用光了，再扛出一袋还不够，请你为国王算一下共要给
达依尔多少小麦？（1袋小麦约1.4×108粒）。
共要给达依尔TTTTT∗TT TT粒小麦
3．程序设计：公元前五世纪我国古代数学家张丘建在《算经》一书中提出了“百鸡问题”：鸡翁一值钱五，鸡母一值钱三，鸡雏三值钱一。
百钱买百鸡，问鸡翁、母、雏各几何？
实验五：二维图形和三维图形的创建
一、实验目的
1．掌握二维图形的绘制。
2．掌握图形的标注
3．了解三维曲线和曲面图形的绘制。
二、实验内容
1．生成1×10维的随机数向量a，分别用红、黄、蓝、绿色绘出其连图脉冲图、阶梯图和条形图，并分别标出标题“连线图”、“脉冲图”、“阶
梯图”、“条形图”。
2．在同一个图形窗口中，绘制两条曲线T1=2T、T2=(12)T；并分别在靠近相应的曲线处标注其函数表达式。
3．编写一个mcircle(r)函数，调用该函数时，根据给定的半径r，以原点为圆心画一个如图所示的红色空心圆。
（图例半径r=5）
4．（1）绘一个圆柱螺旋线（形似弹簧）图。
圆柱截面直径为10，高度为5，每圈上升高度为1。
如左图所示。
（2）利用（1）的结果，对程序做少许修改，得到如右图所示图形。
实验六：Matlab多项式和符号运算
一、实验目的
1．掌握Matlab多项式的运算。
2．了解符号运算。
二、实验内容
1．将多项式P(x)=(x−2)(x+3)(x−7)(x+1)化为x的降幂排列.
结果为T T−TT T−TTT T+TTT+TT
2．求一元高次方程的根。
T9−5T8−30T7+150T6+273T5−1365T4−820T3+4100T2+576T−2880=0 3．求一元高次方程的根,并画出左边多项式函数在x∈[2,2]区间内的曲线x4−2x2+1=0
4．求多项式f1(x)=x3+3x2+5x+7和f2(x)=8x3−6x2+4x−2的乘积f(x);并求T(T)−T1(T)
的商和余式
T2(T)
乘积:
商和余式:
5 . 求y=T5+tan(4T2)+3的符号导数. 6．用符号运算求实验内容4中的T(x)的表达式。