MATLAB教程练习2
Matlab的操作练习
(5)
5 s,t 5
2、编写程序，绘制右 边图形并加上相应标 注。
3、编写程序，绘制右 边螺旋线图形并加上 相应标注。
4、编写程序，绘制右 边图形并加上相应标 注。
5、绘制下列图形 (1)分别利用ezplot函数和plot函数绘制曲线
x2/3 y2/3 1
(2)分别利用ezplot函数和plot函数绘制曲线
练习2、阅读程序 if A then x= 1 else x=0 end
(a)当A取下面旳数组时，程序输出旳成 果是什么？ (1)A=ones(3) (2)A=eye(3) (3)A=inf (4)A=rand(3)>0.5
(b)假如程序中旳A换成any(A)或者all(A) 呢?
在matalb旳程序编写中，假如能够利用数组运算、矩阵运算替代循环 旳点运算，计算效率会大大提升。
x2 y4 0
(3)分别利用ezplot3函数和plot3函数绘制曲线
x cos(t)2 sin(t); y sin(t)2 cos(t); 0 t 2
z t2 *(2 t)2
6、选择合适旳函数绘制图形 (1) 双扭线 = co（s 2） (2)在复平面上绘制函数图形
z cos t i sin(2t)
(3) z1 cos t i sin(2t) z 2 sin(t) icos(2 t)
绘制图形 z=z1*z2
7、分别利用mesh和ezmesh函数绘制网线图
(1)
z x2 y
(2)
z x2 y2
(3)
z2 x2 2y2
符号运算旳练习
1、利用符号运算解下列方程和方程组
(1)x4-1=0旳实根
(1) 100阶矩阵A=(aij), 其中aij=ij
matlab习题二带答案
matlab习题二带答案Matlab习题二带答案Matlab是一种强大的数值计算和科学编程软件，广泛应用于各个领域的科学研究和工程实践中。
掌握Matlab的使用对于科学家和工程师来说至关重要。
为了帮助大家更好地掌握Matlab，下面将介绍一些常见的Matlab习题，并提供相应的答案。
1. 习题一：计算平均值编写一个Matlab函数，输入一个向量，输出该向量的平均值。
答案：```matlabfunction avg = computeAverage(vector)avg = sum(vector) / length(vector);end```2. 习题二：矩阵操作编写一个Matlab函数，输入两个矩阵A和B，输出它们的乘积C。
答案：```matlabfunction C = matrixMultiplication(A, B)[m, n] = size(A);[~, p] = size(B);C = zeros(m, p);for i = 1:mfor j = 1:pfor k = 1:nC(i, j) = C(i, j) + A(i, k) * B(k, j);endendendend```3. 习题三：图像处理编写一个Matlab函数，输入一张彩色图像，输出该图像的灰度图像。
答案：```matlabfunction grayImage = convertToGray(image)grayImage = rgb2gray(image);end```4. 习题四：数据拟合给定一组数据点(x, y)，编写一个Matlab函数，拟合这些数据点为一条直线，并返回拟合直线的斜率和截距。
答案：```matlabfunction [slope, intercept] = fitLine(x, y)n = length(x);sx = sum(x);sy = sum(y);sxy = sum(x .* y);sxx = sum(x .* x);slope = (n * sxy - sx * sy) / (n * sxx - sx^2);intercept = (sy - slope * sx) / n;end```5. 习题五：数值积分编写一个Matlab函数，输入一个函数f(x)和积分区间[a, b]，输出该函数在该区间上的数值积分结果。
MATLAB教程练习2
1矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=728365219a ，分别对a 进行特征值分解、奇异值分解、LU 分解、QR 分解及Chollesky 分解。
2 将矩阵⎥⎦⎤⎢⎣⎡=5724a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772543 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
4 求解多项式x 3-7x 2+2x +40的根。
5 计算多项式9514124234++--x x x x 的微分和积分。
6 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
7 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。
82的变化趋势，用最小二乘法求解y 。
9矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
10 y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
11 []54321=x ，[]108642=y ，计算x 的协方差、y 的协方差、x 与y 的互协方差。
12 有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi ，用三次样条法进行插值。
因式分解：6555234-++-x x x x >> syms x;>> f=x^4-5*x^3+5*x^2+5*x-6; >> factor(f)ans =(x-1)*(x-2)*(x-3)*(x+1)13 ⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x ex x af ax ，用符号微分求df/dx 。
[整理版]matlab课后习题解答第二章
![[整理版]matlab课后习题解答第二章](https://img.taocdn.com/s3/m/e02c0759f6ec4afe04a1b0717fd5360cba1a8daf.png)
第2章0000符号运算0000习题2及解答00001 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？00003/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1))0000〖目的〗0000不能从显示形式判断数据类型，而必须依靠class指令。
0000〖解答〗0000c1=3/7+0.10000c2=sym(3/7+0.1)0000c3=sym('3/7+0.1')0000c4=vpa(sym(3/7+0.1))0000Cs1=class(c1)0000Cs2=class(c2)0000Cs3=class(c3)0000Cs4=class(c4) 0000c1 =00000.52860000c2 =000037/700000c3 =00000.528571428571428571428571428571430000c4 =00000.528571428571428571428571428571430000Cs1 =0000double0000Cs2 =0000sym0000Cs3 =0000sym0000Cs4 =0000sym00002 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量.0000sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')0000〖目的〗0000● 理解自由符号变量的确认规则。
0000〖解答〗0000symvar(sym('sin(w*t)'),1) 0 0 0 0ans = 0 0 0 0 w 0 0 0 0symvar(sym('a*exp(-X)'),1) 0 0 0 0ans = 0 0 0 0 a 0 0 0 0symvar(sym('z*exp(j*th)'),1) 0 0 0 0ans = 0 0 0 0 z 0 0 0 05求符号矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。
matlab课后习题答案（附图）
matlab课后习题答案（附图）习题2.1画出下列常见曲线的图形y （1）⽴⽅抛物线3x命令：syms x y;ezplot('x.^(1/3)')(2)⾼斯曲线y=e^(-X^2);命令：clearsyms x y;ezplot('exp(-x*x)')(3)笛卡尔曲线命令：>> clear>> syms x y;>> a=1;>> ezplot(x^3+y^3-3*a*x*y)(4)蔓叶线命令：>> clear>> syms x y;>> a=1ezplot(y^2-(x^3)/(a-x))(5)摆线:()()tsin-=,=-by1命令：>> clear>> t=0:0.1:2*pi;>> x=t-sin(t);>>y=2*(1-cos(t)); >> plot(x,y)7螺旋线命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t;>>plot3(x,y,z)(8)阿基⽶德螺线命令：clear>> theta=0:0.1:2*pi;>> rho1=(theta);>> subplot(1,2,1),polar(theta,rho1)(9) 对数螺线命令：cleartheta=0:0.1:2*pi;rho1=exp(theta);subplot(1,2,1),polar(theta,rho1)(12)⼼形线命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nnn n3→命令：>>syms n>>limit((n^3+3^n)^(1/n)) ans = 3（2）)121(lim n n n n ++-+∞→命令：>>syms n>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0) ans = 1/2 （4）)(coslimcm xx ∞→命令：syms x m ； limit((cos(m/x))^x,x,inf) ans = 1（5）)111(lim 1--→exx x命令：syms x>> limit(1/x-1/(exp(x)-1),x,1) ans =(exp(1)-2)/(exp(1)-1) （6）)(2lim x x xx -+∞>> limit((x^2+x)^(1/2)-x,x,inf)ans = 1/2练习2.41. 求下列不定积分，并⽤diff 验证：（1）+x dxcos 1>>Clear >> syms x y >> y=1/(1+cos(x)); >> f=int(y,x) f =tan(1/2*x) >> y=tan(1/2*x); >> yx=diff(y ,x); >> y1=simple(yx) y1 =1/2+1/2*tan(1/2*x)^2 (2)+exdx1clear syms x yy=1/(1+exp(x));f=int(y,x) f =-log(1+exp(x))+log(exp(x)) syms x yy=-log(1+exp(x))+log(exp(x)); yx=diff(y,x); y1=simple(yx) y1 = 1/(1+exp(x)) (3)dx x x ?sin 2syms x yy=x*sin(x)^2; >> f=int(y,x) f =x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2 clearsyms x y y=x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2; yx=diff(y,x); >> y1=simple(yx) y1 = x*sin(x)^2 (4)xdx ?sec3syms x y y=sec(x)^3;f=int(y,x) f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)) clear syms x yy=1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)); yx=diff(y,x); y1=simple(yx) y1 =1/cos(x)^32. 求下列积分的数值解 1）dx x-10clearsyms xy=int(x^(-x),x,0,1) y =int(x^(-x),x = 0 .. 1) vpa(y,10) ans =1.291285997 2）xdx e x cos3202?πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clear syms xy=int((1/(2*pi)^(1/2))*exp(-x^2/2),x,0,1) y =7186705221432913/36028797018963968*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/0,2*pi) y =22/65*exp(pi)^4-22/65vpa(ans,10)(3)dx xe21221-π>> clear >> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1); >> vpa(y,14) ans =.341344746068552（4）>> clear >> syms x>> y=int(x*log(x^4)*asin(1/x^2),1,3); Warning: Explicit integral could not be found. > In sym.int at 58 >> vpa(y,14) ans = 2.45977212823752（5） >> clear >> syms x1判断下列级数的收敛性，若收敛，求出其收敛值。
Matlab的操作练习课件
6、选择合适的函数绘制图形 (1) 双扭线 =co（ s 2） (2)在复平面上绘制函数图形
zcostisin(2t)
(3) z 1 c o s t i s i n ( 2 t )z 2 s i n ( t ) i c o s ( 2 t )
绘制图形 z=z1*z2
7、分别利用mesh和ezmesh函数绘制网线图
Matlab的操作练习
Matlab的操作练习
练习
1、利用赋值语句输入数组
1、A=[1,1,1,1；2,2,2,2；3,3,3,3；…,10,10,10,10];
2、B=[21,22.23,24,…,210];
3、 3 2
0
C
0
3
O
O O 2
0
0
3
100100
4、
0
1 1
D
N
2
99 N N
x2/3 y2/3 1
(2)分别利用ezplot函数和plot函数绘制曲线
x2 y4 0
(3)分别利用ezplot3函数和plot3函数绘制曲线
x cos(t)2 sin(t);
y s i n ( t ) 2 c o s ( t ) ; 0t2
z t2 * (2 t)2
Matlab的操作练习
6
7
8
9
y
1.2 2.5 4.1 4.3 4.6 5.2 7.8 10 12.2
Matlab的操作练习
5、已知
5 3 1
A
3
5
3
1 3 5
6、用help命令找出矩阵函数工具箱 查找并使用你熟悉的矩阵操作的对 函数。
1、求A的行列式、特征值和特征向量； 7、绘制曲线
matlab练习题
matlab练习题Matlab是一种高级的技术计算软件，被广泛应用于科学与工程领域。
通过编写Matlab代码，我们可以解决各种数学问题、进行数据分析和可视化、进行仿真等。
本文将提供一些针对初学者的Matlab练习题，帮助读者熟悉Matlab的基本语法和常用操作。
练习一：矩阵运算在Matlab中，矩阵运算是一项重要的操作。
请编写Matlab代码完成以下任务：1. 创建两个3×3的矩阵A和B，元素的值自行选择。
2. 计算矩阵A和B的和，并将结果存储在变量C中。
3. 计算矩阵A和B的乘积，并将结果存储在变量D中。
4. 计算矩阵A的转置，并将结果存储在变量E中。
提示：使用Matlab的矩阵运算符和内置函数可以简化代码编写过程。
练习二：函数编写与图形绘制Matlab是一个功能强大的数学软件，可以用来编写自定义的函数并进行图形绘制。
1. 编写一个Matlab函数，该函数可以计算任意输入整数n的阶乘，并将结果输出。
2. 调用该函数，计算并输出5的阶乘。
3. 使用Matlab绘制函数y=2x+1的图像，其中x的取值范围为-10到10。
提示：使用Matlab的plot函数可以实现图形的绘制，并可使用grid函数添加刻度线。
练习三：数据分析Matlab提供了丰富的数据分析工具和函数，可以用来处理和分析实验数据。
1. 编写一个Matlab脚本，实现对给定实验数据进行统计分析的功能。
2. 实验数据保存在文本文件data.txt中，每行包含一个数据点。
读取数据并将其存储在一个向量中。
3. 计算数据的均值、方差和标准差，并将结果输出。
4. 使用Matlab的直方图函数绘制数据的频率分布直方图。
提示：使用Matlab的readmatrix函数可以读取文本文件中的数据。
调用mean、var和std函数可以计算均值、方差和标准差。
利用histogram函数可以绘制直方图。
练习四：数值解与绘图Matlab具有强大的数值计算和绘图功能，可以用来求解方程和绘制函数图像。
matlab习题——2
matlab习题——2第二章1。
%3X3矩阵a1=eye(3)a2=ones(3)a3=zeros(3)a4=1-2*rand(3)a5=2*randn(3)+1%15X8矩阵b1=eye(15,8)b2=ones(15,8)b3=zeros(15,8)b4=1-2*rand(15,8) b5=2*randn(15,8)+1 2。
%生成a矩阵x=diag([2,-7,8]); xx=diag([3,5],-1); a=rot90(x+xx)%生成b矩阵x=diag([2 5 8]);xx=rot90(diag([7 0 4])); b=x+xx%变成行向量ah=reshape(a,1,9)bh=reshape(b,1,9)3。
a=5-10*rand(50,2);b=round(10*a)/10;c=num2str(b,'%+10.1f') 4。
alpha=-pi:pi/180:pi;x=sin(alpha);y=cos(alpha);%画曲线plot(alpha,x,'b-')hold onlot(alpha,y,'r-'）5。
a=rand(10);b=a>=0.6;total=sum(sum(b))6。
a=randn(10);b=and(a<.5,a>-.5);total=sum(sum(b))7。
%输入a和ba=input('please input a:'); b=input('please input b:'); %分支结构if and(a<1,b<=.5)disp('语句1')elseif and(a<1,b>.5)disp('语句2')elseif and(a>1,b<=.5)disp('语句3')elsedisp('语句4')end8。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1矩阵⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡=728365219a ，分别对a 进行特征值分解、奇异值分解、LU 分解、QR 分解及Chollesky 分解。
2 将矩阵⎥⎦⎤⎢
⎣⎡=5724a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦
⎤
⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺
序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即
⎥⎥
⎥⎥⎦
⎤⎢⎢⎢
⎢⎣⎡237
912685
574 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []296531877254
3 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
4 求解多项式x 3-7x 2+2x +40的根。
5 计算多项式9514124234++--x x x x 的微分和积分。
6 解方程组⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
7 求欠定方程组⎥⎦
⎤
⎢⎣⎡=⎥⎦⎤⎢
⎣⎡5865394742x 的最小范数解。
8
2的变化趋势，用最小二乘法求解y 。
9矩阵⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
10 y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
11 []54321=x ，[]108642=y ，计算x 的协方差、y 的协方差、x 与y 的互协方差。
12 有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi ，用三次样条法进行插值。
因式分解：6555234-++-x x x x >> syms x;
>> f=x^4-5*x^3+5*x^2+5*x-6; >> factor(f)
ans =
(x-1)*(x-2)*(x-3)*(x+1)
13 ⎥
⎥⎦⎤
⎢⎢
⎣⎡=)sin()log(12
x x e
x x a
f ax ，用符号微分求df/dx 。
14 求代数方程组⎪⎩
⎪⎨⎧=+=++00
2y x c by ax 关于x,y 的解。
15 符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t 的变化范围为[0,2π]。
16 绘制极坐标下sin(3*t)*cos(t)的图形。
17有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。
18 在17题结果图中添加标题-at e =y ，并用箭头线标识出各曲线a 的取值。
19 在17题结果图中添加标题-at e =y 和图例框。
20表中列出了4个观测点的6次测量数据，将数据绘制成为分组形式和堆叠形式的条形图。
21 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。
22 2
2
y x
xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图
形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图。
23用sphere 函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。
24编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。
25 有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等级的程序。
1, 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x 轴,y 轴,标题为“正弦余弦函数图象”.
2, 画出2
2
22)
sin(y
x y x z ++=所表示的三维曲面。
y x ,的取值范围是]8,8[-
3,绘制椭圆44922
2=++z y x 表面图,x 在[-2,2],y 在[-2,2]范围.
4,绘制环面:⎪⎩
⎪
⎨⎧=∈+=∈+=u z v v
u y u v u x sin )2,0(sin )cos 1()2,0(cos )cos 1(ππ
1,求下面的优化问题: min -5x1+4x2+2x3 6x1-x2+x3<=8 x1+2x2+4x3<=10 3>=x1>=-1; 2>=x2>=0; x3>=0;
2,某车间生产A 和B 两种产品.为了生产A 和B,所需的原料分别为2个和3个,而所需要的工时分别为4个和2个单位,现在可以应用的原料为100个单位,工时为120个单位,每生产一台A 和B 分别可获得利润6元和4元,应当安排生产A,B 各多少台,才能获得最大的利润?
3,min )12424()(2212
22
11++++=x x x x x e x f x
x 1+x 2=0;
1.5+x 1x 2-x 1-x 2<=0; -x 1x 2-10<=0;
4,对以下数据分别作二次,三次多项式拟合,并画出图形. x=1:16;
y=[4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5, 10.55, 10.58, 10.6];。