Matlab实验第一次实验答案
matlab、实验一答案
在 command 窗口上输入： >> a=@cos; x=[0,2*pi]; >> fplot(a,x) Enter 后弹出图形窗口：
（ 2 ） 在 同 一 坐 标 系 中 绘 制 余 弦 曲 线 y cos(t 0.25) 和 正 弦 曲 线
y sin(t 0.5), t [0,2 ]
6.基本程序设计 （1）编写命令文件：计算 1 2 n 2000 时的最大 n 值。 function[]=program_6() clear all; sum=0; n=1; i=1;
for i=1:100; sum=sum+i; if sum<2000 n=i; else n=i-1; break; end end
在 command 窗口输入： >> t=0:0.01:2*pi; >> plot(t,cos(t-0.25),'-g') >> hold on; >> plot(t,sin(t-0.5)) Enter 后出现下图结果，其中绿色线条是余弦曲线，蓝色线条是正弦曲线：
5.基本绘图控制，绘制 [0,4 * ] 区间上 x1 10 sin t 曲线，并要求 （1）线形为点划线，颜色为红色，数据点标记为加号。 （2）坐标轴控制：显示范围、刻度线、比例、网格线。 （3）标注控制：坐标轴名称、标题、相应文本。 >> t=0:0.1:4*pi; >> x=10*sin(t); >> plot(t,x,'-.+r'); %设置线形为点划线，数据点为加号，红色 >> axis([0,13,-11,11]) %设置坐标轴显示范围，横轴[0,13]，纵轴[-11,11] >> set(gca,'xtick',[0,0.5*pi,pi,1.5*pi,2*pi,2.5*pi,3*pi,3.5*pi,4*pi]) %设置横轴刻度线 >> set(gca,'yminortick','on'); %设置纵轴刻度线 >> grid on %显示网格线 >> xlabel('t 轴');ylabel('x 轴');title('x=10sint'); %设置坐标轴名称、标题 >> text(2*pi,0,'\leftarrow 10sin(2\pi)','fontsize'2 3;4 5 6;7 8 9] >> a=A(1:2,3) a= 3 6 >> b=A(2:3,:) b= 4 5 6 7 8 9 3.多项式
MATLAB实验一：参考--答案
save wenjian.mat A a x load wenjian.mat A a x 小结及思考 在 matlab 中，A*B 与 A .*B 以及 B./A 与 B.\A 之间有什么区别？
A*B 是 A 与 B 进行通常意义上的矩阵乘法得到的矩阵，而 A.*B 是矩阵 A 的每个 元素与 B 的对应位置的元素相乘得到的矩阵。 B./A 是矩阵 B 的每一个元素除以 A 对应位置的元素得到的矩阵， B.\A 是矩阵 A 的每一个元素除以 B 对应位置的元素得到的矩阵。
（1）A=2015*eye(10) （2）A=8*(ones(8)-eye(8))
0 8 （2） 8
8 0 8
8 8 0 88
9. 设有分块矩阵 A
E33 O23
R32 ，其中 E,R,O,S 分别为单位阵、随机阵、零阵和对角阵， S 22
7. 写出下列函数的作用： clc Clear ones det triu fliplr zeros cat
清除命令窗口的记录 清除已声明的变量和函数的记录 生成全 1 矩阵 求一个方阵的行列式 取一个矩阵的上三角 对矩阵进行左右旋转 产生一个零矩阵 拼接两个数组或者矩阵
8. 建立如下矩阵
0 0 2015 2015 0 0 （1） 0 0 2015 1010
3．
计算 1.369 sin
2
7 26.48 2.9 的值（输入程序） 10
1.369^2+sin(7/10*pi)*26.48^(1/2)/2.9
4．
1 3 3 1 1 2 输入矩阵 A 4 4 6 , B 2 2 2 ，指出下列命令的含义 6 8 9 3 3 3
matlab上机实验答案-整理版
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clcclearx=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 自定数据）clccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi;rho=a*cos(b+n*theta);polar(theta,rho)3. 列出求下列空间曲面交线的程序clcclearx=[-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.^2hold onplot3(x1,y1,z3,'*')⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2212y x z -=y x z 322-=4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clcclearx=-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 x clccleara=[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数学实验答案（全）
MATLAB数学实验答案（全）第⼀次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作⼆维、三维⼏何图形，能够⽤Matlab 软件解决微积分、线性代数与解析⼏何中的计算问题。
补充命令vpa(x,n) 显⽰x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形在下⾯的题⽬中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//⾼阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项（要求将结果⽤向量的形式给出）。
同济matlab选修课第一次实验内容参考答案
z1=poly2str(z,'x'); w=polyval(z,x); plot(x,y,'o',x,w,'r')
实验体会与总结
-3-Βιβλιοθήκη 的值，8 3 2 8, b2 3, r2 15
1 10, b1 , r1 28
在三维空间上画出解轨线。 解：
equation.m equation_run.m
function dx=equation(t,x,a,b,c); dx=zeros(3,1); dx=[-a*x(1)+a*x(2);c*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-b*x(3)];
实 验 报 告
班级 实 验 名 称 学号 姓名 Matlab 第一次试验 1．熟悉 Matlab 的操作环境和基本操作； 实 验 目 的 2. 掌握随机矩阵的构建； 3．掌握微分方程数值求解； 4. 掌握曲线拟合 参加实验时间 实 验 地 点 年 月 日 7、8 节课 电信楼 313
实验内容及其结果（写出源程序及其输出结果）
>> x0=[0,1e-10,1]; a1=10;b1=8/3;c1=28; [t,x]=ode45(@equation,[0,100],x0,[],a1,b1,c1); a2=8;b2=3;c2=15; [t,y]=ode45(@equation,[0,100],x0,[],a2,b2,c2); subplot(1,2,1);plot3(x(:,1),x(:,2),x(:,3)); subplot(1,2,2);plot3(y(:,1),y(:,2),y(:,3))
%
[i,j]= find(abs(A)>3);
MATLAB)课后实验答案
实验一 MATLAB 运算基础1、 先求下列表达式得值，然后显示MATLAB 工作空间得使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0、5:2、5 解：4、 完成下列操作：(1) 求[100,999]之间能被21整除得数得个数。
(2) 建立一个字符串向量，删除其中得大写字母。
解：(1) 结果：(2)、 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果就是：实验二 MATLAB 矩阵分析与处理1、 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵与对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5、 下面就是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。
(2) 将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。
(3) 计算系数矩阵A 得条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1、 求分段函数得值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y 值。
（完整版）MATLAB）课后实验答案[1]
![（完整版）MATLAB）课后实验答案[1]](https://img.taocdn.com/s3/m/306ff29b690203d8ce2f0066f5335a8102d266cb.png)
（完整版）MATLAB）课后实验答案[1]实验⼀ MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。
解：(1) 结果：(2). 建⽴⼀个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验⼆ MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证2 2E R RS A O S +??=。
解: M ⽂件如下；5. 下⾯是⼀个线性⽅程组：1231112340.951110.673450.52111456x x x ??=???(1) 求⽅程的解。
(2) 将⽅程右边向量元素b 3改为0.53再求解，并⽐较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M ⽂件如下：123d4e56g9实验三选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ?+-<≠-?=-+≤<≠≠??--?且且及其他⽤if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解：M ⽂件如下：2. 输⼊⼀个百分制成绩，要求输出成绩等级A、B、C、D、E。
Matlab实验一答案
MATLAB 实验一答案1.3 先自定义一个变量，然后分别用8种不同的数字显示格式显示查看。
>> a=pia =3.1416>> format long>> aa =3.14159265358979>> format short>> aa =3.1416>> format short e>> aa =3.1416e+000>> format long e>> aa =3.141592653589793e+000>> format hex>> aa =400921fb54442d18>> format bank>> aa =3.14>> format +>> aa =+>> format rat>> aa =355/1131.4 下面的语句用于画出函数()0.22x y x e -=在[0,10]区间的值x = 0:0.1:10;y = 2*exp(-0.2*x);plot(x,y)1.5 用Matlab 编辑器创建一个m 文件，把上述语句写入这个m 文件并命名为“test1.m ”，保存在当前路径中，然后在命令窗中键入test1，观察结果和运行程序后工作空间的变化.工作空间和结果均与1.4一样1.6 如何清空工作区间数据？键入 clear ；如何关闭图像窗口？键入close ；除了在命令窗输入文件名，还可以怎样运行一个m 文件程序？点击file ，打开m 文件，点击Run 按钮，运行m 文件程序。
1.7 通过以下两种方式得到关于exp 函数的帮助：(1) 在命令窗中输入help exp 命令；(2) 运用帮助空间窗口。
思考，用什么指令可以直接打开帮助空间中关于exp 函数的说明？键入doc exp1.8 假设x =3，y = 4，用Matlab 计算下列表达式：(1) ()232x y x y - (2) 43x y (3) 24x x π- (4) 33x x x y- >> x=3,y=4;>> x^2*(y^3)/(x-y)^2ans =576>> 4*x/(3*y)ans =1>> 4/x*(pi*x^(-2))ans =0.4654>> x^3/(x^3-y^x)ans =-0.72971.9 在当前目录下创建一个m文件，键入以下程序并保存，先把文件保存为“2.m”，运行后观察结果，总结m文件的文件名（包括Matlab标识符）命名规则。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一Matlab使用方法和程序设计
一、实验目的
1、掌握Matlab软件使用的基本方法；
2、熟悉Matlab的数据表示、基本运算和程序控制语句
3、熟悉Matlab绘图命令及基本绘图控制
4、熟悉Matlab程序设计的基本方法
二、实验内容：
1、帮助命令
使用help命令，查找sqrt（开方）函数的使用方法；
解：sqrt
Square root
Syntax
B = sqrt(X)
Description
B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.
Remarks
See sqrtm for the matrix square root.
Examples
sqrt((-2:2)')
ans =
0 + 1.4142i
0 + 1.0000i
1.0000
1.4142
2、矩阵运算
（1）矩阵的乘法
已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B
解：A=[1 2;3 4 ];
B=[5 5;7 8 ];
A^2*B
（2）矩阵除法
已知A=[1 2 3;4 5 6;7 8 9];
B=[1 0 0;0 2 0;0 0 3];
A\B,A/B
解：A=[1 2 3;4 5 6;7 8 9 ];
B=[1 0 0;0 2 0;0 0 3 ];
A\B,A/B
（3）矩阵的转置及共轭转置
已知A=[5+i,2-i,1;6*i,4,9-i];
求A.', A'
解：A=[5+1i,2-1i,1;6*1i,4,9-1i ];
A1=A.',A2=A'
（4）使用冒号表达式选出指定元素
已知：A=[1 2 3;4 5 6;7 8 9];
求A 中第3列前2个元素；A 中所有列第2，3行的元素；
方括号[]
解：A=[1 2 3;4 5 6;7 8 9 ];
B=A([1,2],[3]),C=A(2:end, : )
用magic 函数生成一个4阶魔术矩阵，删除该矩阵的第四列
3、多项式
（1）求多项式42)(3--=x x x p 的根
解：A=[1 0 -2 -4];
B=roots(A)
（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征多项式;
解：
A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4];
A=poly(A);
A=poly2sym(A)
把矩阵A作为未知数代入到多项式中；
4、基本绘图命令
（1）绘制余弦曲线y=cos(t)，t∈[0，2π]
解：t=0:.1:2*pi;
y=cos(t);
plot(t,y),grid
（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]t=0:.1:2*pi;
y1=cos(t-0.25);
y2=sin(t-0.25);
plot(t,y1，t,y2)
grid
5、基本绘图控制
绘制[0，4π]区间上的x1=10sint曲线，并要求：
（1）线形为点划线、颜色为红色、数据点标记为加号；
（2）坐标轴控制：显示范围、刻度线、比例、网络线
（3）标注控制：坐标轴名称、标题、相应文本；
t=0:.01:4*pi;
y=10*sin(t);
plot(t,y,'-.',t,y,'r')
grid
6、基本程序设计
（1）编写命令文件：计算1+2+…+n<2000时的最大n值；
解法1：
s=0;i=0;
while(s<2000) i=i+1;s=s+i; end
s=s-i,i=i-1
解法2：s=0;
for i=1:1000; s=s+i;if(s>2000) ,break;
end
end
s=s-i,i=i-1
（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
解法1：function [ s i ] = mitifang1(n)
s=0;
for i=0:n
s=s+2^i; if(i>n), break;
end
end
s,n
解法2：function [ s i ] = mitifang2(n)
s=0;i=0;
while(i<=n) ,s=s+2^i;i=i+1;
end
s,n
（3）如果想对一个变量x自动赋值。
当从键盘输入y或Y时（表示是），x自动赋为1；
当从键盘输入n或N时（表示否），x自动赋为0；输入其他字符时终止程序。
reply=input('ÇëÊäÈë×Öĸ: ','s');
while reply=='y'||reply=='Y'||reply=='n'||reply=='N'
if (reply=='y'||reply=='Y')
x=1;
disp('x=')
disp(x)
elseif (reply=='n'||reply=='N')
x=0;
disp('x=')
disp(x)
else
disp('shu ru you wu')
end
reply=input('ÇëÊäÈë×Öĸ£º ','s');
end
disp('shu ru you wu ')
解：
三、实验报告要求：
编写实验内容中的相关程序在计算机中运行，程序、运行结果及相关图形一并写在报告上。