MATLAB实验题答案
matlab上机实验答案-整理版
第一次实验答案
1. 设
要求以0.01秒为间隔,求出y 的151个点,并求出其导数的值和曲线。 clc
clear
x=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
clear
x=[-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 on
surf(X,Y,z2)
k=find(abs(z1-z2)<0.5);
x1=X(k)
y1=Y(k)
z3=x1.^2-2*y1.^2
hold on
plot3(x1,y1,z3,'*')
⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2
212y x z -=y x z 322-=
4、设 ⎥⎦⎤⎢⎣
⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线,要求有图形标注。
MATLAB 课后实验答案
实验一 MATLAB 运算基础
1. 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012
2sin851z e =+ (2) 221ln(1)2z x x =++,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦
(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022
a a e e a z a a --+=++=--
(4) 2242011
122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:0.5:2.5 解: M 文件:
z1=2*sin(85*pi/180)/(1+exp(2))
x=[2 1+2*i;-.45 5];
z2=1/2*log(x+sqrt(1+x^2))
a=-3.0:0.1:3.0;
z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)
t=0:0.5:2.5;
z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 运算结果:
z1=2*sin(85*pi/180)/(1+exp(2))
x=[2 1+2*i;-.45 5];
z2=1/2*log(x+sqrt(1+x^2))
a=-3.0:0.1:3.0;
z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)
MATLAB)课后实验答案
实验一 MATLAB 运算基础
1、 先求下列表达式得值,然后显示MATLAB 工作空间得使用情况并保存全部变量。 (1) 012
2sin 851z e =+
(2) 21
ln(2z x =+,其中2120.45
5i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022
a a e e a z a a --+=++=--L (4) 224201
1
122123t 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 O
S +⎡⎤=⎢
⎥⎣⎦。
解: M 文件如下;
5、 下面就是一个线性方程组:
1231112340.951
110.673450.521
11456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣
⎦⎢⎥⎢⎥⎣⎦
(1) 求方程得解。
(2) 将方程右边向量元素b 3改为0、53再求解,并比较b 3得变化与解得相对变化。
MATLAB)课后实验答案
实验一 MATLAB运算基础 1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0122sin851ze (2) 221ln(1)2zxx,其中2120.455ix (3) 0.30.330.3sin(0.3)ln,3.0,2.9,,2.9,3.022aaeeazaaL (4) 2242011122123ttzttttt,其中t=0:0.5:2.5 解: M文件: z1=2*sin(85*pi/180)/(1+exp(2)) x=[2 1+2*i;-.45 5]; z2=1/2*log(x+sqrt(1+x^2)) a=-3.0:0.1:3.0; z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) t=0:0.5:2.5; z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 4. 完成下列操作: (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量,删除其中的大写字母。 解:(1) 结果: m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量 例如: ch='ABC123d4e56Fg9';则要求结果是: ch='ABC123d4e56Fg9'; k=find(ch>='A'&ch<='Z'); ch(k)=[] 实验二 MATLAB矩阵分析与处理 1. 设有分块矩阵33322322ERAOS,其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵,试通过数值计算验证22ERRSAOS。 解: M文件如下; 5. 下面是一个线性方程组: 1231112340.951110.673450.52111456xxx ch = 123d4e56g9 (1) 求方程的解。 (2) 将方程右边向量元素b3改为0.53再求解,并比较b3的变化和解的相对变化。 (3) 计算系数矩阵A的条件数并分析结论。 解: M文件如下: 实验三 选择结构程序设计 1. 求分段函数的值。 2226035605231xxxxyxxxxxxx且且及其他 用if语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。 解:M文件如下: 2. 输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。其中90分~100分为A,80分~89分为B,79分~79分为C,60分~69分为D,60分以下为E。 要求: (1) 分别用if语句和switch语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。 解:M文件如下 3. 硅谷公司员工的工资计算方法如下: (1) 工作时数超过120小时者,超过部分加发15%。 (2) 工作时数低于60小时者,扣发700元。 (3) 其余按每小时84元计发。 试编程按输入的工号和该号员工的工时数,计算应发工资。 解:M文件下 实验四 循环结构程序设计 1. 根据2222211116123nL,求π的近似值。当n分别取100、1000、10000时,结果是多少? 要求:分别用循环结构和向量运算(使用sum函数)来实现。 解:M文件如下: 运行结果如下: 2. 根据11113521ynL,求: (1) y<3时的最大n值。 (2) 与(1)的n值对应的y值。 解:M—文件如下: 3. 考虑以下迭代公式: 1nnaxbx 其中a、b为正的学数。 (1) 编写程序求迭代的结果,迭代的终止条件为|xn+1-xn|≤10-5,迭代初值x0=1.0,迭代次数不超过500次。 (2) 如果迭代过程收敛于r,那么r的准确值是242bba,当(a,b)的值取(1,1)、(8,3)、
MATLAB程序实验参考答案
MATLAB程序设计与应用(第二版)实验参考答案%实验一MATLAB运算基础
%第一题
%(1)
z1=2*sin(85*pi/180)/(1+exp(2))
%(2)
x=[2,1+2i;-0.45,5];
z2=0.5*log(x+sqrt(1+x.^2))
%(3)
a=-3.0:0.1:3.0;
z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)
%(4)
t=0:0.5:2.5;
z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)
%第二题
A=[12 34 -4;34 7 87;3 65 7];
B=[1 3 -1;2 0 3;3 -2 7];
A+6*B
A-B+eye(size(A))
A*B
A.*B
A^3
A.^3
A/B
B\A
[A,B]
[A([1,3],:);B^2]
%第三题
A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]
C=A*B
F=size(C)
D=C(F(1)-2:F(1),F(2)-1:F(2))
whos
%第四题
%(1):
A=100:999;
B=rem(A,21);
C=length(find(B==0))
%(2):
A='lsdhKSDLKklsdkl';
MATLAB)课后实验答案
实验一 MATLAB 运算基础
1. 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012
2sin851z e =+ (2) 221
ln(1)2z x x =++,其中2120.45
5i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022
a a e e a z a a --+=++=-- (4) 2242011
122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩
,其中t =0:0.5:2.5 解: M 文件:
z1=2*sin(85*pi/180)/(1+exp(2))
x=[2 1+2*i;-.45 5];
z2=1/2*log(x+sqrt(1+x^2))
a=-3.0:0.1:3.0;
z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) t=0:0.5:2.5; z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)
运算结果:
z1=2*sin(85*pi/180)/(1+exp(2))
x=[2 1+2*i;-.45 5];
z2=1/2*log(x+sqrt(1+x^2))
a=-3.0:0.1:3.0;
z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)
MATLAB全部实验及答案
MATLAB全部实验及答案
实验一、MATLAB基本操作
实验内容及步骤
1、命令窗口的简单使用
(1)简单矩阵的输入
(2)求[12+2×(7-4)]÷32的算术运算结果
2、有关向量、矩阵或数组的一些运算
(1)设A=15;B=20;求C=A+B与c=a+b?
(2)设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与
A.*B?
A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘(3)设a=10,b=20;求i=a/b=0.5与j=a\b=2?
(4)设a=[1 -2 3;4 5 -4;5 -6 7]
请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全
下标的形式),并将其单下标转换成全下标。
clear,clc
a=[1 -2 3;4 5 -4;5 -6 7];
[x,y]=find(a<0);
c=[];
for i=1:length(x)
c(i,1)=a(x(i),y(i));
c(i,2)=x(i);
c(i,3)=y(i);
c(i,4)=(y(i)-1)*size(a,2)+x(i);
end
c
(5)在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?前面那
个是虚数矩阵,后面那个出错
(6)请写出完成下列计算的指令:
a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?
a^2= 22 16 16
25 26 23
26 24 28
matlab实验报告答案
matlab实验报告答案
篇一:参考答案Matlab实验报告
实验一Matlab基础知识
一、实验目的:
1. 熟悉启动和退出Matlab的方法。
2. 熟悉Matlab命令窗口的组成。
3. 掌握建立矩阵的方法。
4. 掌握Matlab各种表达式的书写规则以及常用函数的使
用。
二、实验内容:
1. 求[100,999]之间能被21整除的数的个数。(rem)
2. 建立一个字符串向量,删除其中的大写字母。(find)
3. 输入矩阵,并找出其中大于或等于5的元素。(find)
4. 不采用循环的形式求出和式i?1
三、实验步骤:
? 求[100,199]之间能被21整除的数的个数。(rem)
1. 开始→程序→Matlab
2. 输入命令:
?m=100:999;
?p=rem(m,21);
?q=sum(p==0)
ans=43
? 建立一个字符串向量,删除其中的大写字母。(find)
1. 输入命令:
?263i的数值解。(sum)
?k=input( ’,’s’);
Eie48458DHUEI4778
?f=find(k =’A’ k =’Z’);
f=9 10 11 12 13
?k(f)=[ ]
K=eie484584778
? 输入矩阵,并找出其中大于或等于5的元素。(find)
1. 输入命令:
?h=[4 8 10;3 空间。
2. 熟悉M文件与M函数的编写与应用。
3. 熟悉Matlab的控制语句。
4. 掌握if,switch,for等语句的使用。
二、实验内容:
1. 根据y=1+1/3+1/5+??+1/(2n-1),编程求:y 5时最
Matlab上机实验答案
Column 61
0.3388 + 0.0000i
t 2 0 t 1 2 (4) z4 t 1 1 t 2 ,其中 t=0:0.5:2.5 2 t 2t 1 2 t 3
>> t=0:0.5:2.5; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2.*t+1 ) z4 =
>> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2))
z2 =
0.7114 - 0.0253i 0.2139 + 0.9343i
0.8968 + 0.3658i 1.1541 - 0.0044i
e0.3a e0.3a 0.3 a sin(a 0.3) ln , a 3.0, 2.9, , 2.9, 3.0 (3) z3 2 2
0.2000 0.1667 0.1429 0.1250 0.1111
1 1 1 1 1
1 2 3 4 5
1 3 6 10 15
1 4 10 20 35
1 5 15 35 70
MATLAB实验答案
2,
(1)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
(3) 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 ];
B=A([1,2],[3]),C=A(2:end, : )
(5) M=magic(4),M(:,4)=[]
3,
(1) p1=[1 0 2 4];PS1=poly2str(p1,'x')
r=roots(p1)
(2) p=poly(A)
polyval(p,20)
poly2str(p,’A’)
4,
(1) t=0:.1:2*pi;
y=cos(t);
plot(t,y),grid
(2) t=0:.1:2*pi;
y1=cos(t-0.25);
y2=sin(t-0.25);
plot(t,y1,t,y2)
grid
t=0:.01:4*pi;
y=10*sin(t);
plot(t,y,'r:+')
grid
axis([0,10,-15,15])
xlabel('x','FontSize',16)
ylabel('y','FontSize',16)
title('正弦函数')
6,
(1)
sum=0;i=1;
while(sum<2000)
sum=sum+i;
i=i+1;
end
n=i-2
{ n=0;s=0;
while s<2000
n=n+1;
s=s+n;
end
n=n-1
}
(2)
解法1:function [ s i ] = mitifang1(n)
MATLAB-课后实验答案
实验一 MATLAB 运算基础
1. 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012
2sin851z e =+
(2) 21ln(2z x =,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦
(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022
a a e e a z a a --+=++=--
(4) 2242011
122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:: 解:
2. 已知:
1234413134787,2033657327A B --⎡⎤⎡
⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦
求下列表达式的值:
`
(1) A+6*B 和A-B+I (其中I 为单位矩阵)
(2) A*B 和A.*B
(3) A^3和A.^3
(4) A/B 及B\A
(5) [A,B]和[A([1,3],:);B^2] 解:
运算结果:
3. 设有矩阵A 和B
12
3
4
53
166789101769,11
121314150
23416171819209
7021
2223242541311A B ⎡⎤⎡⎤
⎢⎥⎢⎥-⎢⎥⎢
⎥⎢⎥⎢⎥==-⎢⎥⎢
⎥
⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦
(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。
解:. 运算结果:
4. 完成下列操作:
)
(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量,删除其中的大写字母。解:(1) 结果:
MATLAB)课后实验答案-精简版
实验一 MATLAB 运算基础
1. 先求下列表达式的值,然后显示MA TLAB 工作空间的使用情况并保存全部变量。 (1) 0
122sin 851z e
=+
(2) 21ln(2
z x =+
,其中2120.45
5i x +⎡⎤
=⎢
⎥-⎣⎦
(3) 0.30.330.3sin(0.3)ln
,
3.0, 2.9,,2.9,3.02
2
a
a
e
e a z a a --+=
++=--
(4) 2
24201
1
122123
t t z t t t t t ⎧≤<⎪
=-≤<⎨⎪-+≤<⎩
,其中t =0:0.5:2.5
2. 已知:
123441
3134
787,2
033
65
73
2
7A B --⎡⎤
⎡⎤
⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦
⎣⎦
求下列表达式的值:
(1) A+6*B 和A-B+I (其中I 为单位矩阵) (2) A*B 和A.*B
(3) A^3和A.^3
(4) A/B及B\A
(5) [A,B]和[A([1,3],:);B^2]
3. 设有矩阵A 和B
123453
0166789101769,1112
1314150
234
16171819209702122
23
24
254
13
11A B ⎡⎤
⎡⎤
⎢⎥⎢⎥-⎢
⎥⎢⎥⎢⎥⎢⎥==-⎢
⎥
⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦
⎣⎦
(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×
2子矩阵赋给D 。 (3) 查看MA TLAB 工作空间的使用情况。 4. 完成下列操作:
(1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量,删除其中的大写字母。 (2). 建立一个字符串向量 例如:
Matlab实验习题集答案
1:用以上两种形式计算36sin 5e ++π算术运算结果。
>> 5^6+sin(pi)+exp(3)
ans =
1.5645e+004
>> x=5^6+sin(pi)+exp(3)
x = 1.5645e+004
2:已知矩阵
⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=2211,2121B A ,对它们做简单的关系与逻辑运算C=(A<B)&(A= =B)
>> A=[1 2;1 2];
>> B=[1 1;2 2];
>> C=(A<B)&(A==B) C =
0 0
0 0
3:对数7sin 5+=a 用五位定点、十五位定点以及有理数形式表示出来。
>> a=5+sin(7);
format short,a
a =
5.6570
>> a=5+sin(7);
>> format long,a
a =
5.656986598718789
4:直接输入创建矩阵
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=98760154321A
>> A=[1 2 3;4 15 60;7 8 9]
A =
1 2 3
4 1
5 60
7 8 9
5:输入矩阵111111111⎛⎫ ⎪ ⎪ ⎪⎝
⎭。 %利用MATLAB 命令直接输入矩阵
OneMatrix=ones(3,3,1)
OneMatrix =
1 1 1 1 1 1 1 1 1
6:输入矩阵
00000 00000⎛⎫ ⎪⎝⎭
>> OneMatrix=ones(2,5,1);
ZeroMatrix=zeros(size(OneMatrix))
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 lim
x mx mx x →-与3
sin lim
x mx mx
x →∞- syms x
limit((902*x-sin(902*x))/x^3) ans =
366935404/3
limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos
1000
x
mx
y e =,求''y syms x
diff(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 计算
22
11
x y e
dxdy +⎰⎰
dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.1394
1.4 计算4
2
2
4x dx m x +⎰ syms x
int(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 =求//高阶导数
参考答案Matlab实验报告
实验一 Matlab基础知识
一、实验目的:
1.熟悉启动和退出Matlab的方法。
2.熟悉Matlab命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握Matlab各种表达式的书写规则以及常用函数的使
用。
二、实验内容:
1.求[100,999]之间能被21整除的数的个数。(rem)
2.建立一个字符串向量,删除其中的大写字母。(find)
3.输入矩阵,并找出其中大于或等于5的元素。(find)
4.不采用循环的形式求出和式
63
1
2i
i=
∑
的数值解。(sum)
三、实验步骤:
●求[100,199]之间能被21整除的数的个数。(rem)
1.开始→程序→Matlab
2.输入命令:
»m=100:999;
»p=rem(m,21);
»q=sum(p==0)
ans=43
●建立一个字符串向量,删除其中的大写字母。(find)
1.输入命令:
»k=input('’,’s’);
Eie48458DHUEI4778
»f=find(k>=’A’&k<=’Z’);
f=9 10 11 12 13
»k(f)=[ ]
K=eie
●输入矩阵,并找出其中大于或等于5的元素。(find)
1.输入命令:
»h=[4 8 10;3 6 9; 5 7 3];
»[i,j]=find(h>=5)
i=3 j=1
1 2
2 2
3 2
1 3
2 3
●不采用循环的形式求出和式的数值解。(sum)
1.输入命令:
»w=1:63;
»q=s um(2.^w)
q=1.8447e+019
- 1 -
实验二 Matlab 基本程序
一、 实验目的:
1. 熟悉Matlab 的环境与工作空间。
Matlab实验题目与答案
(温馨提示:实验课结束后,请将所有作业(题目、代码、结果)利用word 整
理成一个完整的实验报告,加上封面,打印,纸质档于18周周一交)
第一次上机作业
目的:
1. 掌握MATLAB 各种表达式的书写规则
2. 运行课堂上讲过的例子,熟悉矩阵、表达式的基本操作和运算。
作业:
1. 熟悉matlab 集成环境界面。回答以下问题,并操作相关的指令:
(1) 分别写出清除命令窗口和清除变量的指令。
答: clc 和clear
(2)在命令行输入命令后,matlab 的搜索过程是怎样的?
答: (1)检查该命令对象是不是一个变量。
(2)检查该命令对象是不是一个内部函数。
(3)检查该命令对象是否为当前目录下的程序文件。
(4)检查该命令对象是否为MATLAB 搜索路径中其他目录下的M 文件。
(3)什么是matlab 的当前工作目录?写出两种设置当前工作目录的方法? 答: 就是matlab 当前文件读取和存储的默认路径
(1)在当前目录窗口中更改
(2)在MATLAB 桌面工具栏中更改
(3)使用cd 命令:cd c:\mydir---将c :\mydir 设置为当前目录
(4)什么是matlab 的搜索路径?写出两种设置搜索路径的方法?
答: 指Matlab 运行文件时进行搜索的目录。
(1)用path 命令设置:
(2)用Set Path 对话框设置
(5)help 命令和doc 命令有什么作用,它们有什么区别?
答: help 命令:最基本的帮助命令,查询信息直接显示在命令窗口。
doc 命令:在帮助窗口中显示HTML 帮助文档,显示函数的详细用法及 例子,比help 命令更详细。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
result5 =
( 1 ) a = 1 : 2 : 5
a =
1 3 5
( 2 ) b = [ a' , a' , a' ;a ]
b =
1 1 1
3 3 3
5 5 5
1 3 5
( 3 ) c = a + b ( 2 , : )
c =
4 6 8
2、下列运算是否合法,为什么如合法,
结果是多少
>> result2=a*b
Error using *
Inner matrix dimensions must agree. >> result3=a+b
result3 =
3 6 2
58 11
>> result4=b*d
result4 =
31 22 22
40 49 13
31 22 22
40 49 13
-5 -8 7
>> result6=a.*b result6 =
2 8 -3
415 30
>> result7=a./b result7 =
>> result8=
Attempt to reference field of non-structure array.
>> result9=a.\b result9 =
>> result10=a92
result10 =
1 4 9
16 25 36
>> resultl 1=29a
result11 =
2 4 8
16 32 64
>>
result5=[b;c']*d 3、用MATLAB求解下面的的方程组。
1、求以下变量的值,并在MATLAB^验证。
1 2 x1
3 2 x2
11 5 x3
2 1
3 x4
>> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0]
>> B=B'
>> x=inv(A)*B
>> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5]
>> x2=inv(A1)*B2
7 2 1 2
9 15 3 2
2 2 11 5
1 3
2 13
(1)求矩阵A的秩(rank)
(2)求矩阵 A 的行列式(determinant)
(3)求矩阵 A 的逆(inverse)
(4)求矩阵 A 的特征值及特征向量
(eigenvalue and eigenvector)
>> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3)
>> a=det(A3) >> [V,D]=eig(A3)
10
n 10
查看y 的值)
m1=0;
for m=-10:10 m仁m1+2^m;
end
m1
m1 =
6、求分段函数的值。
用if 语句实现,算出下列表中x 对应的y 值。
x=input('enter x=');
if x<0
y=x A2+x-6;
elseif x>=0&&x<5
y=xA2-5*x+6;
else
y=xA2-x-1;
end
y
7、分别用if 和switch 语句实现,将百分
制成绩转换为成绩等级A、B、C、D、E。
其中90~1 00分为A,80~89 分为B,70~79
分为C,60~69 分为D,60 分以下为E。
对超出百分制范围的成绩,给出错误提示
信息。
if 结构程序:
x=input('please enter score=');
if x>=90&&x<=100
9
2
10
disp('A')
7 2
9 15
(1)
2 2
1 3 4 7 1 0
A
4、已知
2n 2 10 29
elseif x<80&&x>=70 otherwise
disp('C')
disp('error')
elseif x<70&&x>=60
end
disp('D')
8、思考题
elseif x<60&&x>=0
设计程序,完成成两位数的加、减、乘、除
disp('E') 四则运算,即产生两个两位随机整数,再输入一个运算符号,做相应的运算,
else 并显示相应的结果。
disp('error')
x=input(' 请输入运算符')
end
a=num2str(floor(rand(1)*90+10));
switch 结构程序:
a
x=input('please enter sco re=');
b=num2str(floor(rand(1)*90+10));
switch fix(x/10)
b
case{10,9}
if x=='+'
if x>100
y=a+b;
disp('error')
elseif x=='-'
else
y=a-b;
elseif x=='*'
disp('A')
y=a*b;
end
elseif x=='/'
case{8}
y=a/b;
disp('B')
else
case{7}
disp('error')
disp('C')
end
case{6}
y
disp('D') 9、启动MATLAB后,点击File|New|M-File ,
启动MATLAB的程序编辑及调试器
case{0,1,2,3,4,5} ( Editor/Debugger ),编辑以下程序,点
击File|Save 保存程序,注意文件名