matlab习题第八章
第八章
实验指导
1、
>> syms x y;
>> A=x^4-y^4;
>> factor(A)
ans =
(x - y)*(x + y)*(x^2 + y^2)
>> a=5135;
>> factor(a)
ans =
5 13 79 2、(1)
>> syms x
>> f=(x-2)/(x^2-4);
>> limit(f,x,2)
ans =
1/4
（2）
>> syms x
f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1); limit(f,x,-1,'right')
ans =
-Inf
3、
（1）
>> syms x y
>> f=sin(1/x);
>> diff(f)
ans =
-cos(1/x)/x^2
>> diff(f,x,2)
ans =
(2*cos(1/x))/x^3 - sin(1/x)/x^4
>> f=(1-cos(2x))/x;
(2)
>> f=(1-cos(2*x))/x
f =
-(cos(2*x) - 1)/x
>> diff(f)
ans =
(2*sin(2*x))/x + (cos(2*x) - 1)/x^2 >> diff(f,x,2)
ans =
(4*cos(2*x))/x - (4*sin(2*x))/x^2 - (2*(cos(2*x) - 1))/x^3 4、
(1)、
>> x=sym('x');
>> f=sqrt(exp(x)+1);
>> int(f)
ans =
atan((exp(x) + 1)^(1/2)*i)*2*i + 2*(exp(x) + 1)^(1/2) (2)
>> x=sym('x');
y=sym('y');
f=x/(x+y);
>> int(f,y)
ans =
x*log(x + y)
(3)
>> syms x;
>> f=exp(x)*(exp(x)+1)^2;
>> int(f,0,log10(2))
ans =
(exp(5422874305198591/18014398509481984)*(3*exp(5422874305198 591/18014398509481984) + exp(5422874305198591/9007199254740992) + 3))/3 - 7/3
(4)
>> sym('x');
f=x*log10(x);
int(f,1,exp(1))
ans =
(18733482797859000068490812234738*log(3060513257434037/112589 9906842624) -
8099090798701270632748702911993)/(50706024009129176059868128 21504*log(10))
5、(1)
>> s=symsum((-1)^(n+1)*1/n,1,inf)
s =
log(2)
(2)
>> syms n
>> s=symsum(x^(2*n-1)/(2*n-1),1,inf)
s =
piecewise([real(n) < 0, zeta(1 - 2*n)/(2*n - 1)]) 6、(1)
>> syms x;
>> f=(exp(x)+exp(-1*x))/2;
>> taylor(f,x,5,0)
ans =
x^4/24 + x^2/2 + 1
(2)
>> syms x;
>> f=sqrt(x^3-2*x+1);
>> taylor(f,x,6,0)
ans =
- x^5/8 - x^4/8 - x^2/2 - x + 1
7、(1)
>> syms x a;
x=solve('x^3+a*x+1','x')
x =
((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3) - a/(3*((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3))
(3^(1/2)*(a/(3*((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3)) + ((a^3/27 + 1/4)^(1/2) -
1/2)^(1/3))*i)/2 + a/(6*((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3)) - ((a^3/27 +
1/4)^(1/2) - 1/2)^(1/3)/2
a/(6*((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3)) - (3^(1/2)*(a/(3*((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3)) + ((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3))*i)/2 - ((a^3/27 + 1/4)^(1/2) - 1/2)^(1/3)/2
(2)
>> syms x
>> x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')
x =
- 227.29614717392215108027554049443 + 2.607041909691930054078218073756*i （3）
>> syms x y;
>> [x y]=solve('log10(x/y)=9','exp(x+y)=3','x','y') x =
(1000000000*log(3))/1000000001
y =
log(3)/1000000001
8、
思考练习
1、(1)数值计算
>> 1/2+2/3
ans =
1.1667
符号运算
>> p1=sym('1/2'); >> a=sym('2/3'); >> p1+a
ans =
7/6
(2)数值计算>> pi+sqrt(5)
ans =
5.3777 >> a=sym('36'); >> y=36;
>> (1+sqrt(a))/2
ans =
7/2>> c=sym('pi'); >> d=sym('sqrt(5)'); >> c+d
ans =
pi + 5^(1/2)
(3)
数值计算
>> (1+sqrt(36))/2
ans =
3.5000
符号计算
>> a=sym('36'); >> y=36;
>> (1+sqrt(a))/2
ans =
2、
>> syms x B1 B2 a b
s1=2*((cos(x))^2)-(sin(x))^2;
s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);
s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1);
h1=simplify(s1)
h2=simplify(s2)
h3=simplify(s3)
h4=simplify(s4)
h1 =
3*cos(x)^2 - 1
h2 =
sin(B1 - B2)
(2^(1/2)*((a + (a^2 - b)^(1/2))^(1/2) + (a - (a^2 - b)^(1/2))^(1/2)))/2 h4 =
2*x + 3
3、
>> syms x a;
f=abs(x)/x;
limit(f,x,0,'left')
ans =
-1
(2)
>> f=(x+a/x)^x;
limit(f,x,inf)
ans =
4、(1)
>> syms x y m
f=sqrt(x+sqrt(x+sqrt(x)));
m=diff(f,'x')
diff(m,'x')
m =
((1/(2*x^(1/2)) + 1)/(2*(x + x^(1/2))^(1/2)) + 1)/(2*(x + (x +
x^(1/2))^(1/2))^(1/2))
ans =
- ((1/(2*x^(1/2)) + 1)/(2*(x + x^(1/2))^(1/2)) + 1)^2/(4*(x + (x +
x^(1/2))^(1/2))^(3/2)) - (1/(8*x^(3/2)*(x + x^(1/2))^(1/2)) + (1/(2*x^(1/2)) +
1)^2/(4*(x + x^(1/2))^(3/2)))/(2*(x + (x + x^(1/2))^(1/2))^(1/2))
(2)
>> syms x y;
f=x+y-sqrt(x^2+y^2);
z1=diff(f,x)
z2=diff(f,y)
z1 =
1 - x/(x^
2 + y^2)^(1/2)
z2 =
1 - y/(x^
2 + y^2)^(1/2)
5、
>> syms x;
>> f=1/sin(x);
>> int(f)
ans =
log(tan(x/2))
(2)
>> sym x;
f=1/(asin(x)^2*(1-x^2)^(1/2)); int(f)
ans =
-1/asin(x)
5、(1)
>> syms x
>> f=1/(1+sqrt(x));
>> int(f,0,4)
ans =
4 - log(9)
(2) sym x;
f=x^3*sin(x)^2/(x^6+2*x^4+1);
int(f,-1,1)
Warning: Explicit integral could not be found. ans =
int((x^3*sin(x)^2)/(x^6 + 2*x^4 + 1), x == -1..1) 7、(1)
>> syms n
s=symsum(1/4^n,1,inf)
s =
1/3
(2)、
>> sym n;
s=symsum(((n+1)/n)^(1/2),1,inf) eval(y)
s =
sum(((n + 1)/n)^(1/2), n == 1..Inf) ans =
y
8、（1）
>> syms x
f=tan(x);
taylor(f,x,3,0)
x
(2)、
>> syms x
f=sin(x)^2;
taylor(f,x,5,0)
ans =
- x^4/3 + x^2
9、(1)
>> syms x
x=solve('log(1+x)-5/(1+sin(x))=2','x')
x =
521.67926389905839979437366649258
(2)
syms x y z
>> [x y
z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x' ,'y','z')
1/2
- 20.0*z^5 - 21.617977528089887640449438202247*z^4 - 1.2359550561797752808988764044944*z^3 +
3.6741573033707865168539325842697*z^2 -
0.058988764044943820224719101123596*z -
0.011235955056179775280898876404494
y =
1/2
156.0*z^5 + 110.82022471910112359550561797753*z^4 + 10.404494382022471910112359550562*z^3 +
2.5337078651685393258426966292135*z^2 +
0.30056179775280898876404494382022*z +
0.097612359550561797752808988764045
z =
1/2
Z
10、
>> [x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')
x =
-(exp(-t*(3*5^(1/2) + 2))*((3*5^(1/2))/2 - (5^(1/2)*(3*5^(1/2) + 5))/6 +
5*exp(t*(3*5^(1/2) - 2))*exp(t*(3*5^(1/2) + 2))*(5^(1/2)/6 - 1/2) +
3*5^(1/2)*exp(t*(3*5^(1/2) - 2))*exp(t*(3*5^(1/2) + 2))*(5^(1/2)/6 - 1/2) + 5/2))/5
y =
exp(-t*(3*5^(1/2) + 2))*((5^(1/2)*(3*5^(1/2) + 5))/30 - exp(t*(3*5^(1/2) - 2))*exp(t*(3*5^(1/2) + 2))*(5^(1/2)/6 - 1/2))
matlab课后习题解答第二章
第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)
ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t
matlab基础练习题带答案
Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。
matlab课后答案完整版
ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z = 2、a=::; >> b=exp*a).*sin(a+ 3、x=[2 4; 5]; y=log(x+sqrt(1+x.^2))/2 y =4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素 A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = + d=c+a*b/(a+b) d = + 第二章 二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans = 第三章 二1（1）
matlab课后习题及答案详解
第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？
matlab习题及答案
2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来， 赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的 形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。 （1）)/1sin(t ，其中)1,1(-∈t ； （2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三 视图 7. 试求出如下极限。 （1）x x x x 1)93(lim +∞ →； （2）1 1lim 0-+→→xy xy y x ； （3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。 （1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ； （2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。 （1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。
matlab第二章习题答案
第一大题： （1） a = 7/3 b = sym(7/3) c = sym(7/3,'d') d = sym('7/3') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d)) a = 2.3333 b = 7/3 c = 2.3333333333333334813630699500209 d = 7/3 v1 = 0.0 v2 = 0.0 v3 = 0.00000000000000014802973661668756666666667788716（2） a = pi/3 b = sym(pi/3) c = sym(pi/3,'d') d = sym('pi/3') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d)) a = 1.0472 b = pi/3 c = 1.047197551196597631317786181171 d = pi/3 v1 = 0.0 v2 = 0.0 v3 = 0.00000000000000011483642827992216762806615818554 (3) a = pi*3^(1/3) b = sym(pi*3^(1/3)) c = sym(pi*3^(1/3),'d') d = sym('pi*3^(1/3)') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d))
a = 4.5310 b = 1275352044764433/281474976710656 c = 4.5309606547207899041040946030989 d = pi*3^(1/3) v1 = 0.00000000000000026601114166290944374842393221638 v2 = 0.00000000000000026601114166290944374842393221638 v3 = 0.0000000000000002660111416629094726767991785515 第二大题： （1） c1=3/7+0.1 c1 = 0.5286 双精度 （2） c2=sym(3/7+0.1) c2 = 37/70 符号 （3） c3=vpa(sym(3/7+0.1)) c3 = 0.52857142857142857142857142857143 完整显示精度 第三大题： （1） findsym(sym('sin(w*t)'),1) ans = w (2) findsym(sym('a*exp(-X)' ) ,1) ans = a (3) findsym(sym('z*exp(j*theta)'),1) ans = z
matlab习题及答案2
MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？ 进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？ 在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？ 在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。
7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。 相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？ 在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。
MATLAB习题参考答案
第1章 MATLAB简介 1、MA TLAB的主要特点有：①语言简洁，编程效率高。②人机界面友善，交互性好。③绘图功能强大，便于数据可视化。④学科众多、领域广泛的MATLAB工具箱。⑤源程序的开放性。 MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。 2、填空题 ⑴命令窗口、命令历史窗口、当前目录窗口 ⑵查阅、保存、编辑 ⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。 3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。 4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。） 6+7+8+9 1+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。） 第2章矩阵与数值数组 1、填空题： ⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。 ⑵全下标、单下标。 2、阅读程序题： （本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果） ⑴ans = 2 3 7 Sa = 10 20 30 A = 1 20 5 30 9 10 4 6 8 10 ⑵ A = 1 3 5 7 9 11 13 15
2 4 6 8 10 12 14 16 ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16 A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 0 ⑶ ans = -1 -4 6 4 ans = 3 0 5 -2 ans = 3 6 9 12 ans = 3 6 9 12 3、A=magic(4)； L=A<10 L = 0 1 1 0 1 0 0 1 1 1 1 0 1 0 0 1 4、省略。 第3章字符串、元胞和构架数组 1. 直接创建法：S1=['Where there is life,' 'there is hope. '] %注意第2行要加入6个空格使其与第一 行字符数相等，否则报错。 S2=str2mat('Where there is life,','there is hope.') S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行 字符的数目 2. 填空题：A（2，3）；A{2，3} 3. 阅读程序题： ⑴subch = ABc12
MATLAB运算基础(第2章)答案
实验01讲评、参考答案 讲评 未交实验报告的同学名单 批改情况： 问题1： 不仔细，式子中出错。 问题2： 提交的过程不完整。 问题3： 使用语句尾分号(;)不当，提交的过程中不该显示的结果显示。 问题4： 截屏窗口没有调整大小。
附参考答案： 实验01 MATLAB 运算基础 （第2章 MATLAB 数据及其运算） 一、实验目的 1. 熟悉启动和退出MATLAB 的方法。 2. 熟悉MATLAB 命令窗口的组成。 3. 掌握建立矩阵的方法。 4. 掌握MATLAB 各种表达式的书写规则以及常用函数的使用。 二、实验内容 1. 数学表达式计算 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 1.1 计算三角函数 12 2sin 851z e =+（注意：度要转换成弧度，e 2如何给出） 示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适 当调整窗口大小。 命令窗口中的执行过程： 《MATLAB 软件》课内实验 王平
1.2 计算自然对数 221 ln(1)2z x x =++，其中2120.45 5i x +??=??-??（提示：clc 命令擦除命令窗口，clear 则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。 应用点乘方） 命令窗口中的执行过程： 1.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 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。 命令窗口中的执行过程：
MATLAB习题及答案 (入门)
习题 1， 计算?? ?? ??=572396a 与??? ???=864142b 的数组乘积。 a.*b ans = 12 36 3 8 42 40 2， 对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。 inv(a)*b ans = -0.5118 4.0427 1.3318 3， 已知：??? ? ? ?????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 a.*a ans = 1 4 9 16 25 36 49 64 81 >> a^2 ans = 30 36 42 66 81 96 102 126 150
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) >> sin(x./pi) ans = -0.1242 0.9826 0.2465 >> cos(x./pi) ans = -0.9923 -0.1857 0.9692 >> tan(x./pi) ans = 0.1252 -5.2915 0.2543 >> cot(x./pi) ans = 7.9894 -0.1890 3.9321 5， 将矩阵?? ????=7524a 、??????=3817b 和??? ???=2695c 组合成两个新矩阵： （1）组合成一个4?3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵 元素，第三列为按列顺序排列的c 矩阵元素，即 ?? ??? ???? ???237 912685 574 [reshape(a,4,1),reshape(b,4,1),reshape(c,4,1)] ans = 4 7 5 5 8 6 2 1 9 7 3 2 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即
matlab答案
第一章 2.在命令窗口输入‘’w=3+2‘’，然后依次使用clear和clc命令，分别观察命令窗口、工作空间窗口和历史命令窗口的变化。 使用clear命令时，命令窗口无变化，工作空间窗口中的内容被删除，历史命令窗口多出一条命令记录。 使用clc命令时，命令窗口中的内容被删除，工作空间窗口无变化，历史命令窗口中多出一条命令记录 3.将硬盘上一已有目录，加入到搜索路径，并将其设置为当前工作目录。 File-set path-add folder-save 第二章 1.计算复数3+4i与5-6i的乘积。 a=3+4i b=5-6i c=a*b 2.构建结构体Students，属性包含Name、age和Email，数据包括｛’Zhang’，18，*‘Zhang@16 https://www.360docs.net/doc/7a9628955.html,’，’Zhang@https://www.360docs.net/doc/7a9628955.html,’+｝、｛’Wang’，21，[]｝和｛’Li’，[],[]｝，构建后读取所有Name属性值，并且修改’Zhang’的Age属性值为19。 Students(1).Name='Zhang' Students(1).Age=18 Students(1).Email='Zhang@https://www.360docs.net/doc/7a9628955.html,','Zhang@https://www.360docs.net/doc/7a9628955.html,' Students(2).Name='Wang' Students(2).Age=21 Students(2).Email=[] Students(3).Name='Li' Students(3).Age=[] Students(3).Email=[] https://www.360docs.net/doc/7a9628955.html, Student(1).Age=19
Student(1).Age 3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵： A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0] ; F=full(A) S=sparse(A) S=sparse([2,1,4],[1,2,4],[1,1,1],4,5) 4.采用向量构造符得到向量[1,5,9....,41]. A=1:4:41 5.按水平和竖直方向分别合并下述两个矩阵： A=[1 0 0;1 1 0;0 0 1],B=[2 3 4;5 6 7;8 9 10] A=[1 0 0;1 1 0;0 0 1] ; B=[2 3 4;5 6 7;8 9 10] ; C=[A B] D=[A;B] 6.分别删除第五题两个结果的第2行。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] C(2,:)=[] D(2,:)=[] 7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] C(2,4:6)=[11 12 13] D(2,:)=[11 12 13] 8.分别查看第5题两个结果的各方向长度
matlab程序设计例题及答案
1.编写程序：计算1/3+2/5+3/7+……+10/21 法一： s=0; for i=1:10 s=s+i/(2*i+1); end s s = 4.4096 法二： sum((1:10)./(3:2:21)) ans = 4.4096 2.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。 s=0; for i=1:100 if mod(i,3)==0&&mod(i,7)==0 s=s+i; end,end s s = 210 3.画出y=n！的图（1<=n<=10）,阶乘的函数自己编写，禁用MATLAB自带的阶乘函数。 x=1:10; for i=1:10 try y(i)=y(i-1)*i; catch y(i)=1; end,end plot(x,y)
12345678910 0.511.522.533.54x 10 6 4.一个数恰好等于它的因子之和，这个数就称为完数。例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。编程找出2000以内的所有完数。 g=[]; for n=2:2000 s=0; for r=1:n-1 if mod(n,r)==0 s=s+r; end end if s==n g=[g n]; end end g g =6 28 496
5.编写一个函数，模拟numel函数的功能，函数中调用size函数。 function y=numelnumel(x) m=size(x); y=m(1)*m(2); numelnumel([1 2 3;4 5 6]) ans = 6 6. 编写一个函数，模拟length函数的功能，函数中调用size函数。 function y=lengthlength(x) m=size(x); y=max(m(1),m(2)); lengthlength([1 2 3;4 5 6]) ans = 3 7.求矩阵rand（5）的所有元素和及各行平均值，各列平均值。 s=rand(5); sum=sum(sum(s)) mean2=mean(s,2) mean1=mean(s) sum = 13.8469
第2章MATLAB数据及其运算习题答案.doc
第2章 MATLAB数据及其运算 习题2 一、选择题 1．下列可作为MA TLAB合法变量名的是（）。D A．合计B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。C A．+10 B．1.2e-5 C．2e D．2i 3．使用语句t=0:7生成的是（）个元素的向量。A A．8 B．7 C．6 D．5 4．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。B A．1 B．2 C．3 D．4 5．已知a为3×3矩阵，则a(:,end)是指（）。D A．所有元素B．第一行元素 C．第三行元素D．第三列元素 6．已知a为3×3矩阵，则运行a (1)=[]后（）。A A．a变成行向量B．a变为2行2列 C．a变为3行2列D．a变为2行3列 7．在命令行窗口输入下列命令后，x的值是（）。B >> clear >> x=i*j A．不确定B．-1 C．1D．i*j 8．fix(354/100)+mod(354,10)*10的值是（）。D A．34 B．354 C．453D．43 9．下列语句中错误的是（）。B A．x==y==3 B．x=y=3 C．x=y==3 D．y=3,x=y 10．find(1:2:20>15)的结果是（）。C A．19 20 B．17 19 C．9 10 D．8 9 11．输入字符串时，要用（）将字符括起来。C A．[ ] B．{ } C．' ' D．" " 12．已知s='显示"hello"'，则s的元素个数是（）。A A．9 B．11 C．7 D．18
13．eval('sqrt(4)+2')的值是（）。B A．sqrt(4)+2 B．4 C．2 D．2，2 14．有3×4的结构矩阵student，每个结构有name（姓名）、scores（分数）两个成员，其中scores是以1×5矩阵表示的5门课的成绩，那么要删除第4个学生的第2门课成绩，应采用的正确命令是（）。D A．rmfield(student(1,2).scores(2)) B．rmfield(student(4).scores) C．student(4).scores(2)=0 D．student(1,2).scores(2)=[] 15．有一个2行2列的单元矩阵c，则c(2)是指（）。B A．第一行第二列的元素内容B．第二行第一列的元素内容 C．第一行第二列的元素D．第二行第一列的元素 二、填空题 1．从键盘直接输入矩阵元素来建立矩阵时，将矩阵的元素用括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用分隔，不同行的元素之间用分隔。中括号，逗号或空格，分号 2．设A=[1,2;3,4]，B=[5,6;7,8]，则A*B= ，A.*B= 。 A*B=[19,22;43,50]，A.*B=[5,12;21,32] 3．有3×3矩阵，求其第5个元素的下标的命令是，求其第三行、第三列元素的序号的命令是。[i,j]=ind2sub([3 3],5)，ind=sub2ind([3 3],3,3) 4．下列命令执行后的输出结果是。20 >> ans=5; >> 10; >> ans+10 5．下列命令执行后，new_claim的值是。This is a great example. claim= 'This is a good example.'; new_claim=strrep(claim,'good','great') 三、应用题 1．命令X=[]与clear X有何不同？请上机验证结论。 Clear X是将X从工作空间中删除，而X=[]是给X赋空矩阵。空矩阵存在于工作空间中，只是没有任何元素。 2．在一个MA TLAB命令中，6 + 7i和6 + 7*i有何区别？i和I有何区别？ 3．设A和B是两个同大小的矩阵，试分析A*B和A.*B、A./B和B.\A、A/B和B\A 的区别？如果A和B是两个标量数据，结论又如何？ 4．写出完成下列操作的命令。 （1）将矩阵A第2～5行中第1，3，5列元素赋给矩阵B。 （2）删除矩阵A的第5号元素。 （3）求矩阵A的大小和维数。 2
matlab练习题和答案
matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];
计算方法及其MATLAB实现第二章作业
作者：夏云木子 1、 >> syms re(x) re(y) re(z) >> input('计算相对误差：'),re(x)=10/1991,re(y)=0.0001/1.991,re(y)=0.0000001/0.0001991 所以可知re(y)最小，即y精度最高 2、 >> format short,A=sqrt(2) >> format short e,B=sqrt(2) >> format short g,C=sqrt(2)
>> format long,D=sqrt(2) >> format long e,E=sqrt(2) >> format long g,F=sqrt(2) >> format bank,H=sqrt(2) >> format hex,I=sqrt(2) >> format +,J=sqrt(2) >> format,K=sqrt(2)
3、 >> syms A >> A=[sqrt(3) exp(7);sin(5) log(4)];vpa(pi*A,6) 4、1/6251-1/6252=1/6251*6252 5、（1）1/（1+3x）-(1-x)/(1+x)=x*(3*x-1)/[(1+3*x)*(1+x)] (2) sqrt(x+1/x)-sqrt(x-1/x)=2/x/[sqrt(x-1/x)+sqrt(x+1/x)] (3) log10(x1)-log(x2)=log10(x1/x2) (4) [1-cos(2*x)]/x =x^2/factorial(2)-x^4/factorial(4)+x^6/factorial(6)-…
MATLAB复习题(答案仅供参考)1
一、填空题 1. MATLAB命令窗口中可用clear 命令清除工作区中的变量；用clc 命令清除命令窗口中的内容。 2. MATLAB中的运算符包括数学运算符、矩阵运算符、数组运算符。 3 MATLAB的含义为矩阵运算。 4. MATLAB中的函数类型有嵌套函数、主函数、子函数、私有函数和匿名函数。 5. MATLAB中用系统预定义的特殊变量nargin 表示函数输入参数个数。 6. size命令的功能确定数组大小。 7. MATLAB中用linspace 符号定义单元数组。 8. MATLAB中用hold on 命令打开图形保持功能。 9 MATLAB中如果一个变量名为123，该变量名是否合法否（填是或否）。 10. MATLAB中用inv 函数完成矩阵的求逆运算，用det 函数求矩阵的行列式。 11.已知C={2+5i,’Hi’;1:4,[1 2]} 输入s=C(3) 即s=__'_Hi__'___. 12.先输入P=struct(‘name’,’Jack‘,’num’,43) ，再输入s= 即s=___Jack___. 13.当x=[ 1 2 3; 4 5 6] , y=[ 5 8; 5 8 ] 时，__不能__(能﹨不能)使用plot 命令绘图。 14．sym 或syms _能_______(能﹨不能)创建符号矩阵。 15.已知x=3:2:24,则x有__11_个元素。 结构的开始是“if”命令，结束是_end____命令。 17.已知x为一个向量，计算ln（x）的MATLAB命令是计算_log(x)___。 18.命令：>>subplot(3,3,x)中有__9_个子图。 19.已知s=‘显示“hey”’，则s的元素个数是___9__。 y=f() global W W=2 y=f1(5) function y=f() global W y=w+W 创建以上函数文件，在命令窗口中运行“y=f”命令则显示__出错___。 21. MATLAB中用find 函数查找数组中的非零元素，用sort 函数对数组进行排序。 22. MATLAB中的M文件包括M脚本文件和M函数文件。 23. MATLAB中用系统预定义的特殊变量eps 表示浮点数的精度。 24. MATLAB中8位有符号整数用int8 定义。 25. MATLAB中小括号、转置运算符和冒号运算符，具有最高优先级的为转置运算符。 26. MATLAB中用系统预定义的特殊变量nargout 表示函数输出参数个数。 27. MATLAB中用% 符号表示其后为程序注释。 28. 已知矩阵A,则矩阵A的共轭转置可表示为A’。 29. MATLAB中用grid on 命令显示坐标轴下的网格线。 30. MATLAB中如果一个函数名为test，该函数名是否合法是（填是或否）。
MATLAB及其在大学物理中的应用——第二章习题答案
MATLAB 及其在大学物理中的应用——第二章习题答案 作者：荆楚理工吴世华 2.1试求下列极限： （1）x x x 1)93(lim +∞→（2）5 232)5()3()2(lim +++∞→+++x x x x x x x 2.2求下列函数的导数： （1） x e x x x y -=1sin )(（2）) 4)(3()2)(1()(----=x x x x x y （3）)ln(tan 22y x x y a +=（4）0,ln 1)(>+-=n x a x na x y n n 2.3已知参数方程22,sin cos cos ln dx y d dx dy t t t y t x 和求? ??-==。2.4设22,00y u xv yu yv xu ?????=+=+求。2.5设已知函数矩阵 ?? ????++=z y x z e x y x f y sin 3),(23，试求其雅可比矩阵。2.6求下列不定积分：
（1）dx a x x a x x f ?++=2222) (3)(（2） dx x x x x x f ?+++=1)1()(（3） dx bx xe x f ax ?=cos )(（4）cxdx bx e x f ax sin sin )(?=2.7求下列函数的泰勒幂级数展开。（1）dt t t x ?0sin （2）x x -+11ln （3） )3/3sin(5π+-x e x 分别关于x=0,x=a 的幂级数展开。2.8分别用roots 函数和多项式伴随矩阵的特征值求根法求解方程 6251234=++-x x x 的所有根。 2.9分别用矩阵除法和linsolve 函数法求解下列方程组。 ???????=++-+=+--=+-+--=+-+-0 2660835502492633092458432142143214321x x x x x x x x x x x x x x x 2.10用符号法求解方程组：? ??=-=+151926628y x y x