中南大学matlab习题复习课0401222

习题课Any函数判断元素是否为非零元素any(v),如果v是非零元素返回true(即1)否则返回flase(即0)如果A是矩阵，则把矩阵的列当做向量来处理，函数返回每个列向量的逻辑值2024/8/1222024/8/123Magic函数2024/8/124线性方程组求解试求数值方程组10x-y=9

-x+10y-2x=7

-2x+10z=6的解2024/8/1252024/8/1262024/8/127九、方程符号求解2024/8/128符号常微分方程MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。语法：dsolve('eq','con','v') %求解微分方程dsolve('eq1,eq2…','con1,con2…','v1,v2…') %求解微分方程组（1）当y是因变量时，微分方程'eq'的表述规定为：y的一阶导数

或

表示为Dy；y的n阶导数

或

表示为Dny。（2）微分初始条件'con'应写成'y(a)=b，Dy(c)=d'的格式；当初始条件少于微分方程数时，在所得解中将出现任意常数符C1、C2……，解中任意常数符的数目等于所缺少的初始条件数。2024/8/129【例

】

求微分方程

，y(1)=0，y(0)=0的解。>>y=dsolve('x*D2y

3*Dy=x^2','x') %求微分方程的通解y=C3*x^4-x^3/3+C2>>y=dsolve('x*D2y

3*Dy=x^2','y(1)=0,y(0)=0','x')%求微分方程的特解y=x^4/3-x^3/3【例

】求微分方程组

的解。>>[x,y]=dsolve('Dx=y,Dy=

x')x=C8*cos(t)+C7*sin(t)y=C7*cos(t)-C8*sin(t)程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量，结果相同：>>[x,y]=dsolve('Dx=y,Dy=

x','t')2024/8/1210非线性方程的符号注解语法：fsolve(‘fun’,x0) （1）fun为所要求解的函数名，常以M文件形式给出；（2）X0为初始向量或矩阵Fc.mfunctiony=fc(x)y(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2));y(2)=x(2)-0.7*cos(x(1))+0.2*sin(x(2));y=[y(1)y(2)]x0=[0.50.5]>>fsolve('fc',x0)ans=0.52650.50792024/8/1211九、方程符号求解1求微分方程 ，y(1)=0，y(0)=0的解a)x^4/5-x^3/4b)x^3/3-x^2/3c)x^4/3-x^3/3d)x^4/5-x^3/4（c）2024/8/12122、利用非线性方程组的符号法求解如下方程组：在（0.5，0.5）附近的解。a)x1=0.5265x2=0.5079b)x1=0.5079x2=0.5265c)x1=0.6265x2=0.7079d)x1=0.7079x2=0.6265（a）2024/8/12131对于圆周率值，若取值为3.1415，则其有数字为a、3位b、4位c、5位d、2位（b）2024/8/12142024/8/12152下面的MATLAB语句中正确的有：a)a=sin(x)b)Record-1=3+4ic)a=sin(pi)d)C=1+6j（cd）2024/8/12163符合MATLAB变量命名规则的有a)var!2015b)VAR15c)my-var2015d)2015_VAR（b）2024/8/12174在循环结构中跳出循环，但继续下次循环的命令为

(a)

return;

(b)

break

(c)

continue

(d)

keyboad（c）2024/8/1218二、多项式及其运算2024/8/1219二、多项式及其运算1求多项式在x=3,8处的值

a)-52-1730b)-54-1732c)-52-1732d)-54-1730（c）2024/8/12202计算多项式x3+11x2+55x+125在x=[12;34]的矩阵计算值

a)[192287;416585]b)[190285;406582]c)[194292;422580]d)[192288;420580]（a）2024/8/12213试用一个二次多项式拟合下列数据x 0.5 1.0 1.5 2.0 2.5 3.0y 1.75 2.45 3.81 4.80 8.00 8.60a)0.46001.86010.9730b)0.84002.94010.85608c)0.49001.86010.9730d)0.49001.25010.8560（d）2024/8/1222三、矩阵及数组代数运算2024/8/1223三、矩阵及数组代数运算1、已知A=[0

9

6;1

3

0;357];B=[1

4

3;1

5

0;259],则A.*Ba)[217554;4193;227272]b)[03618;1150;62563]c)[4193;217554;227272]d)[62563;03618;1150]（b）2024/8/12242、求数组b=[123;456;789]的指数运算a)1.0e+06*[1.11891.37481.6307;2.53393.11343.6929;3.94894.85205.7552]b)1.0e+03*[0.00270.00740.0201;0.05460.14840.4034;1.09662.98108.1031]c)[00.69311.0986;1.38631.60941.7918;1.94592.07942.1972]d)[1.00001.41421.7321;2.00002.23612.4495;2.64582.82843.0000]（b）2024/8/12253求矩阵的特征值a)

[1;2;3]b)

[2;3;4]c)

[3;4;5]

d)

[1;3;5]（a）2024/8/1226四、数组抽取2024/8/1227四、数组抽取1、若A为一二维数组，要删除其第34行的元素，可使用命令

a)A(34,:)=0b)A(:,34)=[]c)A(34,:)=[]d)A(:,34)=0（c）2024/8/12282、设矩阵A=[1,2,3,4;3,4,5,6;3,6,-1,0;9,12,1,2],抽取成向量[2;4;6;12]，则使用命令a)A(:,2)b)A(2,:)’

c)A([3,4,5,6])d)A(2,:)（a）2024/8/1229五、数组关系逻辑运算2024/8/1230五、数组关系逻辑运算1、已知A=[0

9

6;1

3

0;357];B=[1

4

3;1

5

0;259],则A<Ba)[011;110;111]b)[100;010;001]c)[101;100;011]d)[001;010;101]（b）2024/8/12312、已知A=magic(5),A(:,3)=zeros(5,1)，则all(A(:,1)<10)a)0b)1c)[00100]d)[11011]（a）2024/8/1232六、符号求极限1.求极限a）-1/6b)Infc)–Infd）-1（a）2024/8/12332a）-1b)1/2c)0d）1（c）2024/8/1234七、符号求导2024/8/12351y=,求x=3，t=2时的dy/dx值a)[2.0,0;-1.97998,20.0855]b)[2.0,0;-1.98898,20.0865]c)[2.0,0;-1.96998,20.0655]d)[2.0,0;-1.98998,20.0855]（a）2024/8/1236七、符号求导2求函数在指定点的导数值。a)62454b)63484c)43654d)43060（a）2024/8/1237八、符号求积1符号求积