Matlab第一次上机练习
matlab上机习题详细讲解-试题答案解析
学习指导参考P 第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3) y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
第一次上机作业 matlab
2011《工程最优化与（Matlab ）实现》课程第一次Matlab 上机实验实验目的：1. 熟悉Matlab 软件界面及菜单常用操作2. Matlab 命令窗口常用命令熟悉3. Matlab 向量、矩阵操作命令熟悉4. Matlab 帮助使用熟悉上机作业：1. 建立若干变量、向量、矩阵，熟悉who,whos ，what ，which ，save ，load ，clear ，clc 等命令。
2. 生成以4阶魔方阵为实部，4阶帕斯卡矩阵为虚部的复数矩阵，进行求行列式值，求逆，lu 分解，qr 分解（验证q 正交性），转置（包括直接转置和共轭转置），再生成四阶随机方阵与之进行加、减、点乘，矩阵乘，点除（包括左除和右除）等运算。
3. 生成6阶魔方矩阵，寻找15,16所在的行、列，将15,16改为90,91，在矩阵右侧填一列全1向量，将新得的矩阵为单元扩张为一行两列的新矩阵。
4. 求解下列线形方程组的解，并验证，用至少三种方法5. 求以下函数积分，近似到1e-7。
f(x)=sin(x)*x 在（0，inf ）上的积分f(x,y)=sin(x)*exp(-y^2)在x 属于[0,1]，y 属于[0,3]上的积分6.对上述散点使用interp1函数求出0.25:0.1:0.75点函数值7. 试绘制曲线8. 将下列程序用if-elseif-end 编程实现9. 利用inline 法，匿名函数法，函数文件法编制函数f=exp(-x^2-y^2)*sin(x+y^3)/(x+4*y)，并求出在点（2，3），（6，7）处的值 ⎪⎪⎩⎪⎪⎨⎧=++=+-+=+-+=++33222125w z y w z y x w z y x z y x。
matlab入门习题答案
matlab入门习题答案MATLAB入门习题答案MATLAB是一种强大的数学软件工具，广泛应用于工程、科学和计算领域。
它具有丰富的功能和灵活的编程语言，使得用户可以轻松地进行数据分析、可视化、模拟和算法开发。
在学习MATLAB的过程中，习题是一个很好的学习方法，通过练习可以加深对MATLAB的理解和掌握。
下面是一些MATLAB入门习题的答案，希望对初学者有所帮助。
1. 编写一个MATLAB程序，计算1到100之间所有奇数的和。
```matlabsum = 0;for i = 1:2:100sum = sum + i;enddisp(sum);```2. 编写一个MATLAB程序，计算一个矩阵的转置矩阵。
```matlabA = [1 2 3; 4 5 6; 7 8 9];B = A';disp(B);```3. 编写一个MATLAB程序，生成一个包含1到10的整数的向量，并计算向量中所有元素的平均值。
```matlabvec = 1:10;avg = mean(vec);disp(avg);```4. 编写一个MATLAB程序，绘制正弦函数的图像。
```matlabx = 0:0.1:2*pi;y = sin(x);plot(x, y);```通过以上习题的练习，相信大家对MATLAB的基本语法和功能有了更深入的了解。
当然，MATLAB的学习是一个持续的过程，希望大家能够坚持练习，不断提升自己的编程能力。
同时，也欢迎大家分享更多的MATLAB习题答案，让我们一起学习进步。
实验一 MATLAB基础训练
实验一 MATLAB基础训练一、实验目的本次上机实验主要练习使用Matlab的基本操作和基础知识，包括数组（复数、向量、矩阵、结构体数组等）的创建和数组元素的操作和运算、矩阵的运算、Matlab的运算符（尤其是点运算‘.’）、脚本M文件和函数M文件的编写、Matlab文件的编程（基本的流程控制结构）、基本的二维和三维绘图方法以及图形的标注等。
希望通过本次实验使大家尽量在短时间内（4学时）掌握Matlab的基本操作和基础知识，为后面的实验项目奠定基础。
二、实验原理参见PPT中有关内容。
三、实验内容1. 上机练习课件中的例子。
2. 设两个复数a=1+2i, b=3-4i, 计算a+b, a-b, a*b, a/b, a和b的模。
3. 计算下式的结果，其中x=-3.5°, y=6.7°（提示：①应将角度单位由度转换为数学函数所能处理的弧度值；②求根函数sqrt，取绝对值函数abs,具体用法用help查询）4. 对矩阵123456789A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦, 实现下列操作:(1) 左右翻转(fliplr命令)(2) 上下翻转(flipud命令)(3) 利用cat命令分别将A扩展成3×6和6×3的矩阵(3) 分别提取A的第2行,第2列,对角线元素(4) 删除A的第2行2列的元素（提示：将矩阵元素赋空阵[]可以删除元素，注意此时元素的访问只能使用单下标的方式。
观察删除元素后，A中元素的排列方式的变化）5. 创建[0,2π]区间上拥有100个等间隔元素的列向量．．．x．（Matlab默认是行向量），并绘制y=sin(x1/3)的函数图像。
6. 创建如下图所示的单结构体数组。
7. 编程训练：下图所示电路中，R 1=2Ω, R 2=4Ω, R 3=12Ω, R 4=4Ω, R 5=12Ω, R 6=4Ω, R 7=2Ω, u s =10V, 求i 3。
（要求：以脚本M 文件方式建立程序。
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('输入的不是方阵或者矩阵元素不是数!');endfor i=1:afor j=1:bA(i,j)=(-1)^(i+j)*X(i,j)*det(X([1:i-1,i+1:end],[1:j-1,j+1:end]));endendA=A';2．已知表达式y=6*x^5+4*x^3+2*x^2-7x+10，x的范围是[0，100]，使用三阶拟合和五阶的方法得出多项式的表达式，并编程在图中绘制出原曲线、三阶拟合和五段拟合的曲线。
clcclear allx=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 ongrid onplot(x,y1,'b-.',x,y2,'y--*');grid on三条曲线放大图（原曲线与五阶拟合完全重合）：3．读如图像5.tif（450*450），做fourier变换，a)将变换之后所得到的矩阵只保留前10行10列，后10行10列的数据，其他位置都设置为0，然后再作反fourier变换，观察结果；原图：运行程序：clcclear allI=imread('5.tif');figure(1);imshow(I);B=fft2(I);C=ones(size(B));C(11:end-10,11:end-10)=0;B1=B.*C;B1=B1./256;I1=ifft2(B1);figure(2);imshow(I1);图像：b)将变换之后所得到的矩阵将前10行10列，后10行10列的数据置为0，然后再作反fourier变换，观察结果；clcclear allI=imread('5.tif');figure(1);imshow(I);B=fft2(I);C=zeros(size(B));C(11:end-10,11:end-10)=1;B1=B.*C;B1=B1./256;I1=ifft2(B1);figure(2);imshow(I1);c) 从观察的结果你有什么体会。
matlab上机练习（附答案）
1.以下两种说法对吗？（1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
（2）MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗？4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直接键入法？用数组编辑器？用M文件编辑器？5.运用数组算术运算符去掉下面程序里的for/end循环：x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符，重写下面的程序代码：x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7．某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
Matlab第一次上机练习
第一次上机练习：目的：1.Matlab 的安装、卸载：（选作，自己有电脑的可以做） 2.建立自己的工作目录，再将自己的工作目录设置到matlab 的搜索路径下。
3.掌握MATLAB 各种表达式的书写规则 4.运行课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算。
作业：1. 求下列表达式的值，显示MATLAB 工作空间的使用情况并保存全部变量：（1）)1034245.01(26-⨯+⨯=w（2）,)tan(22ac b e abc c b a x ++-+++=ππ 其中a=3.5,b=5,c=9.8。
（3）])48333.0()41[(22απβππα---=y ，其中32.3=α，9.7-=β （4）)1ln(2122t t e z t ++=，其中⎥⎦⎤⎢⎣⎡--=65.05312i t 2. 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=7613870451A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=023352138B ，求下列表达式的值： （1）B A 6+、I B A +-2（其中I 为单位阵）；（2）A*B 、A.*B 、B*A 、B.*A ；（3）A/B 、B\A ；（4）[A, B]、[A([1 3],:)；B^2]。
3. 已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454.9632053256545410778.01023A ，取出其前三行构成矩阵B ，其前两列构成矩阵C ，其右下角3×2子矩阵构成的矩阵D ，B 与C 的乘积构成矩阵E ，分别求E<D 、E&D 、E|D 、~E|~D 。
4. 产生均值为3，方差为1 的5阶正态分布的随机方阵。
5. 使用函数，实现方阵A 左旋 90、右旋 90的功能，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=131074129631185210741A 。
MATLAB上机习题一-2016-解答
4. 在命令窗口输入 ispc，按键盘的 F1 键会弹出帮助，阅读并写出该内置变量的 含义。 （也可通过主菜单中[Help]->[Product Help]菜单打开帮助，并搜索 ispc 来完 成） 解答：
表示 MATLAB 的版本是否是 Windows 系统平台下的版本
5. 观察 MATLAB 命令窗口中关键字、字符串、注释、错误提示的显示颜色。提 示：输入相关表达式来观察，例如字符串可以这样输： ss=’abc’ 解答： 关键字：蓝色 字符串：紫色 注释：绿色需要 MATLAB7.0 以上版本，请确认系统装安装有。第 1、2 机房需要在开机时 选择【APPLY】系统。 截图方法：按键盘上右上角的[Print Screen]键，在作业文档中按 Ctrl+V 键粘贴即可。
1. 打开 MATLAB，观察界面中包含哪些子窗口，写出窗口名字 解答： （不需截图） 命令窗口（ Command Window ） 、工作区窗口（ Workspace ） 、历史命令窗口 （Command History） 、当前文件夹窗口（Current Folder） 、...
u
(e 2）
+ v)
2
v2 − u
3） 解答：
u = 2; v = 3;
u − 3v uv
%1) 4*u*v/log(v) %2) (exp(u)+v)^2/(v^2-u) %3) sqrt(u-3*v)/(u*v) 运行屏幕截图：
9. 计算如下表达式： 1） ( 3 − 5i )( 4 + 2i ) 2） sin ( 2 − 8i ) 解答：
2. 使用 path 命令将计算机的 C 盘添加到系统路径，使用 cd 命令将 C 盘设置为 当前工作路径 解答：(不需截图，仅需要写出相应的指令表达式即可) path('C:\', path) cd('C:\') %或者 path(path, 'C:\')
Matlab第一次上机作业
输入： >>tic, n=9;[u,k]=xsj(n), toc,surf(u)
计算Байду номын сангаас果如下：
u=
0 0 0 0 0 0 0 0 0 0 0 0 0.0256 0.0413 0.0508 0.0560 0.0577 0.0560 0.0508 0.0413 0.0256 0 0 0.0413 0.0686 0.0859 0.0955 0.0986 0.0955 0.0859 0.0686 0.0413 0 0.0508 0.0859 0.1088 0.1216 0.1258 0.1216 0.1088 0.0859 0.0508 0 0 0.0560 0.0955 0.1216 0.1364 0.1412 0.1364 0.1216 0.0955 0.0560 0 0 0.0577 0.0986 0.1258 0.1412 0.1462 0.1412 0.1258 0.0986 0.0577 0 0 0.0560 0.0955 0.1216 0.1364 0.1412 0.1364 0.1216 0.0955 0.0560 0 0 0.0508 0.0859 0.1088 0.1216 0.1258 0.1216 0.1088 0.0859 0.0508 0 0 0.0413 0.0686 0.0859 0.0955 0.0986 0.0955 0.0859 0.0686 0.0413 0 0 0.0256 0.0413 0.0508 0.0560 0.0577 0.0560 0.0508 0.0413 0.0256 0 0 0 0 0 0 0 0 0 0 0 0 0
程序一： function [u,k]=xsbj(n) % xsbj:用块 Jacobi 迭代法求解线性方程组 A*u=f % u:方程组的解； k 迭代次数； n:非边界点数 % a：方程组系数矩阵 Aii 的下对角线元素；b:方程系数矩阵 Aii 的主对角线元素 % c：方程组系数矩阵 Aii 的上对角线元素；d:追赶法所求方程的右端向量 % e:允许误差界；er：迭代误差 f=2*1/(n+1)^(2)*ones(n+2,n+2); a=-1*ones(1,n);b=4*ones(1,n);c=-1*ones(1,n);u=zeros(n+2,n+2);e=0.000000001; for k=1:2000 er=0; ub=u; for j=2:n+1 d(1:n)=f(2:n+1,j)+ub(2:n+1,j-1)+ub(2:n+1,j+1); x=zg(a,b,c,d); u(2:n+1,j)=x'; er=er+norm(ub(:,j)-u(:,j),1); end if er/n^2<e,break; end end 程序二： function x=zg(a,b,c,d) % zg:用追赶法求解线性方程组 n=length(b); % LU 分解 u(1)=b(1); for k=2:n if u(k-1)==0,D=0,return; end l(k)=a(k)/u(k-1); u(k)=b(k)-l(k)*c(k-1); end % 追赶法求解之追过程，求解 Ly=d y(1)=d(1); for k=2:n y(k)=d(k)-l(k)*y(k-1); end % 追赶法求解之追过程，求解 Ux=y if u(n)==0,D=0,return;end x(n)=y(n)/u(n); for k=n-1:-1:1 x(k)=(y(k)-c(k)*x(k+1))/u(k); end 输入：
Matlab基础及应用第一章上机作业
1、新建一个文件夹（自己的名字命名）
答：
2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。
保存，关闭对话框。
使用path命令查看MATLAB搜索路径。
3、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye。
4、使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
5、编写一段简短的脚本文件，保存并显示运行结果。
6、创建两个double型变量a=32、b=5，并计算a+b、a-b、a*b、a/b、a\b；
创建两个int8型变量a、b，取同样数值，并计算a+b、a-b、a*b、a/b、a\b，对于计算结果与前次计算结果不同的情形请给出解释。
7、查看int16数据类型的取值范围(intmin,intmax)；查看单精度数据类型的取值范围和精度(realmin,realmax,eps)。
8、求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量：
(1)sin(60o) (2) e3 (3) cos(3/4π)
9、用两种定义复数的方法计算(直接定义，complex(a,b))
(1)(3-5i)(4+2i) (2) sin(2-8i)。
Matlab上机练习
13sin 25457.0145.175.375.0222-++=---y e x z x y x练习一：1) 以两种方式打开 MA TLAB 工作窗口，进入MATLAB 6.0 的工作环境，并尝试用不同的方式退出。
2) 尝试、熟悉 MA TLAB 6.0 的各栏菜单以及各个工具栏的功能。
3) 重新启动 MATLAB 6.0，进入 MATLAB 工作窗口，用 who 命令查看当前工作空间内有无变量及其值。
4) 绘制函数y=cos(5x+2)/sin(3x+1) 的图像，并求解当 x=2 时的函数值。
5) 此时再次用 who 命令查看工作空间内的变量名及其值，与（ 3）比较，同时用whos 命令查看变量，比较与 who 命令的不同。
6) 熟练掌握 MATLAB 的通用命令。
7) 练习并熟练掌握 MA TLAB 的帮助命令，学会利用 MATLAB 的帮助信息。
8) 用 lookfor 命令查找函数 cos 的信息，并与 help 命令查找的结果相比较，注意采用两种命令之间的差别。
练习二： 1) 创建一个表达式式并求当x =1, y =2时的z 值。
2) 计算i x y 108++=的值，其中451-+-=x 。
3) 创建一个包含“ 她是我们班最好的学生”字符串。
练习三：1) 将如下矩阵 A 进行转置和求逆。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=412303214A 2) 在MATLAB 环境下，用下面三条指令创建矩阵 C ，看输出怎样的结果。
a=2.7358;b=33/79; C =[1,2*a+i*b,b*sqrt(a);sin(pi/4),a+5*b,3.5+i]练习三：3) 用 rand 函数建立一个 3×3 维随机数值矩阵，并将该数值矩阵转换为符号矩阵，比较它们的不同。
4) 求下述矩阵的逆、特征根和矩阵的行列式值。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A ，⎥⎦⎤⎢⎣⎡=22211211a a a a A 练习四： 1) 用曲面图命令 surf 表现函数 22y x z +=的图像。
MATLAB上机练习及作业
MATLAB上机练习及作业第一章1. 熟悉MATLAB环境，练习第一章的命令；2. 练习显示格式（P13）1利用Help命令，搞清楚MATLAB的所有主题，利用help signal和help control 搞清楚这两个工具箱2练习whos,who,disp,clear,which,ver等命令第二章3练习矩阵操作，书上P34～404练习逻辑和关系运算，书上P40～435练习随机数产生，书上P506练习时间和日期函数，P56～577矩阵操作diag,reshape,rot90,fliplr,flipud,tril,triu等8练习MATLAB基本函数P669上交：第二章习题中的No.3, No.6第3章1. 练习产生图3.1～3.5, 3.9～3.12, 3.14, 3.18，3.21，3.22, 3.24～3.312. 第3章习题: (P122) 2, 5--93. 交程序：5；74. 并练习图形旋转、标注、复制等功能5. 交结果：与习题7类似，只是画出你自己上个月的消费饼图第4章1. 上机前好好复习这章内容，练习P137页的两段程序，并将终值扩大到500，体会两者的速度差异。
2. 熟悉MATLAB的调试技术。
3. 第4章习题: (P158)4. 1，3，5，8，95. 交结果：3；5，8（分别求出当I＝54，I＝105时的结果）；9第5章习题(p.222)1, 2, 3, 7, 8, 9, 13, 15, 16, 20交结果: 1, 2, 8, 9, 13, 16第6章P.251 看懂综合设计示例, 设计程序计算每个班的平均成绩, 计算教师的实际学时数。
(交结果)P.269 2, 3第7章1. 在[-1,1]区间上绘制出y=-2x^2+4的曲线, 并在图中标注出最大值点的坐标.2. 任意给定一句英文句子, 请提取其每个单词, 并设计一结构, 其域有Name、No、Length、Value，分别用于存储每个单词的名称、句中序号、单词长度、单词各个字符的ASCII 码。
上机练习专题1：Matlab语言基础
数学物理建模与计算机辅助设计
练习4：Matlab的符号运算（续）
– – – – 运行g1='x+y+z=1';g2='x-y+z=2';g3='2*x-y-z=1'; 再运行f=solve(g1,g2,g3)，观察结果 再运行[f.x, f.y, f.z]，观察结果 运行dsolve('Dx=y','Dy=-x','x(0)=0','y(0)=1')，观察结果
• 求解微分方程
– 运行dsolve('D2y=-a^2*y','y(0)=1','Dy(pi/a)=0') ，观察结果
• 求解微分方程
– 求解方程并绘图
d2 y dy 2 2y 0 2 dx dx dy y 0 1， 0 0 dx
Page 6
数学物理建模与计算机辅助设计
练习3：矩阵的基本运算
– – – – – – – – – – – –
Page 4
在命令窗口中输入A=[1 2 3; 4 5 6; 7 8 9]，回车，观察结果 再输入B=zeros(size(A));C=ones(2,4); D=eye(3);，回车 在命令窗口中输入B=A(1:2,:)，回车，观察结果 再在命令窗口中输入C=repmat(B,2,3)，回车，观察结果 输入A=[1 2 3;2 3 4;3 4 5]; B=[1 1 1;2 2 2;3 3 3]; C=A+B; 输入E=[B,[5 5 5]'];C=A*E;，回车，观察结果 输入A=[5 12 8;6 5 8;9 6 10]; B=[7;11;7]; X=A\B，回车 输入A=[2 1 -3 -1; 3 1 0 7; -1 2 4 -2; 1 0 -1 5]; B=inv(A)回车 输入C=det(A)，回车，观察结果 输入D=[1,2,3;4,5,6;7,8,9]; E=D^3; F=D.^3;，回车观察结果 再输入G=D(1:3,2:3);H=D([1,3],2);，回车观察结果 在命令窗口中输入I=[D;ones(1,3)], J=[D,ones(3,1)]，回车
Matlab上机作业部分参考答案
上机练习二 参考答案
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并 且按照从大到小的顺序排列好！ 【求解】 a=10*rand(1,10)-5; b=sort(a,'descend')
上机练习二 参考答案
2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵A将得出什么 结果？
Matlab 上机课作业
吴梅红 2012.10.15
上机练习一
上机练习一 参考答案
上机练习一 参考答案
上机练习一 参考答案
上机练习二
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并且按 照从大到小的顺序排列好！ 2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵 A将得出什么结果？ 3、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全 部偶数行提取出来，赋给B 矩阵，用A =magic(8) 命令生成A 矩阵，用上述的命令检验一下结果是不是正确。
【求解】用课程介绍的方法可以直接输入这两个矩阵 >> A=[1 2 3 4; 4 3 2 1; 2 3 4 1; 3 2 4 1] A= 1234 4321 2341 3241 若给出A(5,6)=5 命令，虽然这时的行和列数均大于A矩阵当前的维数， 但仍然可以执行该语句，得出 >> A(5,6)=5 A= 123400 432100 234100 324100 000005 复数矩阵也可以用直观的语句输入 3+2i 4+1i; 4+1i 3+2i 2+3i 1+4i; 2+3i 3+2i 4+1i 1+4i; 3+2i 2+3i 4+1i 1+4i]; B= 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i
(完整版)matlab上机题目一
上机一1、（1（2） 查询所有人的成绩。
2、已知：A=[27,24,0,-1;23,35,1,14;e 4,-8,13,20;10,12,19,21;11,18,ln2,0]， 求：（1）求A 中第4、5、2行与第2、4、1列组成的C 阵，并将C 中第2行元素全改成8，得到D 阵；（2） 创建与A 同维的G 阵，使C 阵左右翻转后位于G 中第2行至第4行、第2列至第4列位置，且G 中第1行第1列元素为1，其余元素为0。
（3） 求A 阵除以3的余数阵H 。
3、已知：A=[5,17,-2,3;-1,0,21,11;14,9,-10,8]，B=[2,5,8,-1;12,-2,7,8;11,13,16,4] 求：（1）将A 除以B 阵的对应元素，并将结果朝负无穷大方向取整。
（2）求A./B ， B.\A ，A.*B 和A*B ’的值。
（3）求2*ones(3,4)-eye(3,4)|A-B>B&B-8的值.4、要求在闭区间[0，2π ]上产生50个等距采样的一维数组。
试用两种不同的指令实现。
5、对矩阵，⎥⎦⎤⎢⎣⎡=4321A 问如下MATLAB 指令：)(),(),5.0(^),5.0(.^A sqrtm A sqrt A A 所得的结果相同吗？他们中的哪个结果为复数矩阵？ 6、说出MATLAB 指令A （3，1，2，： ）=1：4所产生数组的维数、大小和长度；然后对A 进行降维处理；最后指出所有非零元素的“全下标”位置。
7、现有五个数组： (2×3)的实数数组A ； (2×3)的复数数组B ； (2×3)的字符串数组C ； (2×3)的异质（元胞）数组D ； (2×3)的结构（构架）数组E 。
问，它们各占多少字节？8、写出使以下这段文字成为字符串的MATLAB 程序。
注意保持这段文字的格式。
在英式用法中，引号通常是单引号，如‘Fire!’。
matlab上机习题详解 试题答案
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
matlab上机练习1
实验1 MATLAB7.0基本编程实验目的：① 掌握脚本和函数；② 掌握matlab 中的变量和M 文件的流控制语句；③ 掌握函数的设计和实现。
实验要求:给出程序和实验结果。
实验内容：一、一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
二、从键盘输入若干个数，当输入10时结束输入，求这些数的平均值和它们之和。
三、求[1000，2000]之间第一个能被17整除的整数。
四、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,10000]之间的全部完数。
五、Fibonacci 数列定义如下：1f =12f =112,(2)n n n f f f n --=+>求Fibonacci 数列的第20项。
六、设计一个猜数游戏。
首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。
根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You w o n”，同时退出游戏。
用户最多可以猜7次。
实验结果：一、 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
程序如下：for m=100:999m1=fix(m/100); %求m 的百位数字m2=rem(fix(m/10),10); %求m 的十位数字m3=rem(m,10); %求m 的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend153370371407二、从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
程序如下：sum=0;cnt=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;cnt=cnt+1;val=input('Enter a number (end in 0):');endif (cnt > 0)summean=sum/cntendEnter a number (end in 0):45Enter a number (end in 0):65Enter a number (end in 0):342Enter a number (end in 0):45Enter a number (end in 0):67Enter a number (end in 0):12Enter a number (end in 0):0sum =576mean =96三、求[100，200]之间第一个能被21整除的整数。
MATLAB上机操作题
实验一 MATLAB 基本操作一、目的和要求（1） 熟练掌握MA TLAB 变量的使用。
（2） 熟练掌握矩阵的创建。
（3） 熟练掌握MA TLAB 的矩阵和数组的运算。
二、内容和步骤1． 创建矩阵矩阵是包括n m ⨯个元素的矩阵结构，矩阵中的元素可以是实数或复数，单个元素构成的标量以及多个元素构成的行向量和列向量都是矩阵的特殊形式。
下面用多种方式创建矩阵。
（1） 直接输入。
例：>> a=[1 2 3;4 5 6;7 8 9] a =1 2 34 5 6 7 8 9练习：从键盘输入一个44⨯的矩阵，矩阵元素自己定义，变量名定义为m1。
（2） 用from:step:to 方式。
例：>> a=[1:3 ; 4:6 ; 7:9]a =1 2 34 5 6 7 8 9练习：从键盘输入一个向量，用from:step:to 方式，要求从零开始，步长为1.5，到20结束，变量名定义为m2。
（3） 用linspace 函数。
例：>> a=[linspace(1,3,3);linspace(4,6,3);linspace(7,9,3)] a =1 2 3 4 5 6 7 8 9 练习：使用logspace 函数创建π4~0的行向量，有20个元素，查看其元素分布情况，变量名定义为m3。
（4） 使用matlab 内部函数创建矩阵。
例：>> a=ones(3)a =1 1 11 1 11 1 1练习：建立一个44⨯的正态分布矩阵，要求其均值为2，方差为2，变量名定义为m4。
2．矩阵的操作（1）矩阵元素的提取。
例：>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9提取元素8，可输入a(3,2)或a(6)。
练习：建立一个44⨯的对角阵，提取它的第三列元素，并将矩阵的第四列元素改为全1。
（2）矩阵的裁剪和拼接。
例：>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9>> b=[2 3 4;5 6 7;8 9 10]b =2 3 45 6 78 9 10>> c=cat(1,a,b)c =1 2 34 5 67 8 92 3 45 6 78 9 10>> c=cat(2,a,b)c =1 2 3 2 3 44 5 6 5 6 77 8 9 8 9 10>> c(:,3:5)=[]c =1 2 44 5 7 7 8 10练习：建立两个任意的44⨯矩阵，将它们沿着第二维拼接成一个新矩阵，并将新矩阵的第3列和第5列删除。
MATLAB上机练习一参考解答
上机练习一参考解答一、实验目的1、 熟悉Matlab 编程2、 体会数学上恒等，算法上不一定恒等二、实验内容1． Using the Taylor polynomial of degree nine and three-digit rounding arithmetic to find an approximationto 5-e by each of the following methods.(A) ∑=--≈905!)5(n n n e , (B) ∑=-≈=9055!5/11n nn e e An approximate value of 5-e correct to three digits is 31074.6-⨯. Which formula, (A) or (B), gives the most accuracy, and why?1） 算法基础利用x e 的Taylor 公式00!!n nk x n n x x e n n ∞===≈∑∑，x -∞<<+∞ （1）及001/1/1/!!n nk x x n n x x e e n n ∞-====≈∑∑，x -∞<<+∞， （2）其中k 是根据精度要求给定的一个参数。
在本题中将k 取为9， x 取为－5或5即可由公式（1）或（2）得到5-e 的近似计算方法（A ）或（B ）。
2） 程序下述程序用公式（A ）及（B ）分别在Matlab 许可精度下及限定在字长为3的算术运算情况下给出5-e的近似计算结果，其中results_1, results_2为用方法（A ）在上述两种情况下的计算结果，err_1, err_2为相应的绝对误差；类似的，results_3, results_4为用方法（B ）在上述两种情况下的计算结果，err_3, err_4为相应的绝对误差；具体程序如下：% Numerical Experiment 1.1 % by Xu Minghua, May 17, 2008 clc; %Initialize the data x=-5; k=9; m=3; %three-digit rounding arithmetic %------------------------------------ % Compute exp(x) by using Method (A) % with the computer precision results_1=1; power_x=1; for i=1:k factor_x=x/i; power_x=power_x*factor_x; results_1=results_1+power_x; end results_1 err_1=abs(exp(x)-results_1)%------------------------------------% Compute exp(x) by using Method (A) % with the 3-digits precisionresults_2=1;power_x=1;for i=1:kfactor_x=digit(x/i,m);power_x=digit(power_x*factor_x,m); results_2=digit(results_2+power_x,m); endresults_2err_2=abs(exp(x)-results_2)%------------------------------------% Compute exp(x) by using Method (B) % with the computer precisiont=-x;results_3=1;power_x=1;for i=1:kfactor_x=t/i;power_x=power_x*factor_x;results_3=results_3+power_x; endresults_3=1/results_3err_3=abs(exp(x)-results_3)%------------------------------------% Compute exp(x) by using Method (B) % with the 3-digits precisiont=-x; results_4=1;power_x=1;for i=1:kfactor_x=digit(t/i,m);power_x=digit(power_x*factor_x,m);results_4=digit(results_4+power_x,m); endresults_4=digit(1/results_4,m)err_4=abs(exp(x)-results_4)%------------------------------------上述主程序用到一个子程序digit.m, digit(x,m)的作用是将x四舍五入成m位数。
Matlab学习教程 第一讲 上机练习
第1章安装及使用前的准备1 本章简介MATLAB是一种功能非常强大的科学计算软件。
在正式使用它之前应对它有一个整体的认识。
本章的基本内容主要包括：（1）MATLAB的简介。
（2）MATLAB的安装。
（3）MATLAB的桌面环境的介绍。
（4）MATLAB的帮助系统。
（5）MATLAB的搜索路径及其扩展的方法。
2 达到的目标（1）熟悉MATLAB的桌面环境极其简单的操作桌面平台包括6个窗口，分别是主窗口、命令窗口、历史窗口、当前目录窗口、发行说明书窗口和工作间管理窗口。
（2）熟悉MATLAB的帮助系统及其操作。
（3）掌握MATLAB的搜索路径及其扩展的方法。
3 具体内容3.1 打开和关闭MATLAB的桌面6个窗口。
3.2 帮助系统的使用（1）联机帮助系统用下面3种方法进入MATLAB联机帮助系统，（a）直接按下MATLAB主窗口中的“？”认按钮。
（b）选中主窗口[Help]下拉菜单的前4项中的任何一项。
（c）在命令窗回中执行helpwin或helpdesk或doc（2）命令窗口查询帮助（a）help系列help plot（b）lookforlookfor plot（c）分析（a）与（b）之间的区别（3）联机演示系统①用下面3种方法进入MATLAB联机演示系统的主演示页面，（a）单击MATLAB主窗口菜单的【Help】→【Demos】选项。
（b）在命令窗口输人demos。
（c）直接在帮助页面上选中Demos选项。
②演示（a）graphics→2D-plots。
（b）graphics→3D-plots。
3.3 MATLAB的搜索路径及其扩展的方法（1）查看搜索路径（a）用搜索路径对话框查看搜索路径（b）用path命令查看搜索路径（c）用genpath命令查看搜索路径（d）用editpath或pathtool命令查看搜索路径（2）扩展搜索路径首先在E：＼路径下建立一个新的目录myfiles，然后将目录myfiles扩展为搜索路径。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Matlab第一次上机练习
一基础知识部分
1打开Matlab，区分工作窗口，历史窗口，工作空间和当前路径
2指定工作路径
二常量与变量（在工作窗口中进行）
1用标量对变量赋值，试试不赋值可不可以定义变量。
学会两个变量值的互换。
比如a=5,b=10。
如何互换。
2查看pi，realmax等常量
3用矢量和矩阵对变量赋值。
输入矢量，用冒号，linespace创建矢量。
输入矩阵，用ones，命令创建矩阵。
4标量，矢量和矩阵的加减乘除，注意要符合线性代数的规律。
比较*和.*的区别。
学会使用关系运算符>, <等等。
5学会访问矩阵中的值。
如矩阵A=[2,4,6,8,10,12,14], 要得到矩阵第一行第三个元素的值是多少，输入A(1,3)即可，要得到第二行第3个元素到第6个元素的值，输入A(2,3:1:6)。
6字符串赋值给变量
三编写脚本文件和函数文件，二者合称M文件。
1计算从1到100所有偶数的和，写出脚本文件
2编写函数文件，计算从1到输入参数之间所有偶数的和。
学会使用function命令，%作注释。
了解输入参数，输出参数。
文件名要和函数名一致。
3函数的调用
3.1调用函数得到函数的返回参数值。
体会下面函数
function avg=average(x)
%求矢量中各数值的平均值
n=length(x) %该函数可以得到矢量含多少个元素。
avg=mean(x,n)
function a=mean(u,d)
a=sum(u)/d %sum命令可以求得矢量各元素的和。
此种方式调用是得到函数的返回值。
3.2调用函数本身。
体会下面函数
Function a=result(x)
a=feval(@zwb,x) %函数句柄。
feval命令的意思是获得zwb这个函数在x处的函数值。
function y=zwb(x)
y=x^3+x^2-3*x-3;
将这种创建函数的方法与匿名函数和内联函数进行对比。
我个人认为有了这种方式，就不需
要学习匿名函数和内联函数了。
匿名函数是以如下方式创建函数：fx=@（x）5*x^2+x-6。
内联函数是以如下方式创建函数：fx=inline(‘5*x^2+x-6’)。
四学会使用流程控制语句
1if语句
2switch语句
3for语句。
注意比较它与c语言的异同。
4while语句
5continue，break，return语句
6灵活运用上面语句对含有10数值的矢量（该矢量为输入参数）进行从小到大的排序。
7一个球从100米的高处落下，每次落地后反跳回原来高度的一半，再又落下。
问，第10次落下着地时，共经过多少米？第十次反跳高度是多少？写成函数文件，可根据输入参数（次数）自动求得上面两问。
8猴子开始有n个桃子，每天吃桃是他有桃个数的一半还多一个，第10天时，猴子想吃桃时就发现只有一个桃子了。
问n是多少？
五绘图(参考教材P223)
1创建两个维数相等的矢量，使用plot命令绘图。
学会使用help。
2使用plot命令时加入更多的参数，对图形进行修饰
3尝试学会使用fplot命令。
（使用help）
四.6.7.8题的M文件内容课后发送一份到我的邮箱，并打印
或抄写一份给我，下次课交。