Matlab第三章题

合集下载

MATLAB基础及其应用教程-周开利-邓春晖课后答案,第三章

MATLAB基础及其应用教程-周开利-邓春晖课后答案,第三章

第三章习题答案1.代码:a=[1 -1 -1]; roots(a)结果:ans =-0.61801.61802.代码:x=0:10;y=sin(x);xi=0:0.15:10; %选取了67个插值点,要增加n,只需减小步长即可y0=sin(xi); %算精确值y1=interp1(x,y,xi); %分段线性插值y2=interp1(x,y,xi,'spline'); %三次样条插值plot(xi,y0,'o',xi,y1,xi,y2,'-.')legend('精确值','分段线性插值','三次样条插值')结果:3.理论公式为:p=1.0332*exp(-(x+500)/7756),所以拟合模型可写为:p=a*exp(-k*x+b) 式中,a, k, b为常数,两边同时取自然对数,得:log(p)=-k*x+b+log(a)问题转化为线性模型。

注意:自然对数是log(x), 以10为底的对数是log10(x)代码:clear;x=[0 300 600 1000 1500 2000];p=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491];lnp=log(p); %转化为 p 的自然对数值,模型转化为线性模型pk=polyfit(x,lnp,1); % 线性拟合,得到模型的斜率pk(1)和常数pk(2) 模型为: p=exp(pk(1)*x)*exp(pk(2))xi=0:50:2000;p0=1.0332*exp(-(xi+500)/7756); %理论值p1=exp(pk(1)*xi+pk(2)); %拟合模型值p2=interp1(x,p,xi,'spline'); %三次样条插值plot(x,p,'p',xi,p0,xi,p1,'--',xi,p2,'-.');legend('测量值','理论值','拟合值','三次样条值');format long % 数据显示格式为15位有效数字x2=0:200:2000 % 取10个点,比较差异pp1=1.0332*exp(-(x2+500)/7756) %理论值pp2=exp(pk(1)*x2+pk(2)) % 拟合值pp3=interp1(x,p,x2,'spline') % 样条插值err1=sum(abs(pp2-pp1).^2) % 拟合值的误差绝对值总和err2=sum(abs(pp3-pp1).^2) % 样条值的误差绝对值总和结果:0200400600800100012001400160018002000从图像上,都符合得很好,但很难看出差异。

matlab课后习题答案第三章

matlab课后习题答案第三章

第3章数值数组及其运算习题3及解答1 要求在闭区间]2,0[ 上产生具有10个等距采样点的一维数组。

试用两种不同的指令实现。

〖目的〗●数值计算中产生自变量采样点的两个常用指令的异同。

〖解答〗%方法一t1=linspace(0,2*pi,10)%方法二t2=0:2*pi/9:2*pi %要注意采样间距的选择,如这里的2*pi/9.t1 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.2832t2 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.28322 由指令rng('default'),A=rand(3,5)生成二维数组A,试求该数组中所有大于0.5的元素的位置,分别求出它们的“全下标”和“单下标”。

〖目的〗●数组下标的不同描述:全下标和单下标。

●sub2ind, int2str, disp的使用。

●随机发生器的状态控制:保证随机数的可复现性。

〖解答〗rng('default')A=rand(3,5)[ri,cj]=find(A>0.5);id=sub2ind(size(A),ri,cj);ri=ri';cj=cj';disp(' ')disp('大于0.5的元素的全下标')disp(['行号 ',int2str(ri)])disp(['列号 ',int2str(cj)])disp(' ')disp('大于0.5的元素的单下标')disp(id')A =0.8147 0.9134 0.2785 0.9649 0.95720.9058 0.6324 0.5469 0.1576 0.48540.1270 0.0975 0.9575 0.9706 0.8003大于0.5的元素的全下标行号 1 2 1 2 2 3 1 3 1 3列号 1 1 2 2 3 3 4 4 5 5大于0.5的元素的单下标1 2 4 5 8 9 10 12 13 153 采用默认全局随机流,写出产生长度为1000的“等概率双位(即取-1,+1)取值的随机码”程序指令,并给出 -1码的数目。

matlab1-8章课后作业

matlab1-8章课后作业

MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件,是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境,已被广泛应用于不同领域。

Matlab的基本功能包括:数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。

Matlab第二章1.创建double的变量,并进行计算。

(1)a=87,b=190,计算a+b、a-b、a*b。

(2)创建uint8 类型的变量,数值与(1)中相同,进行相同的计算。

>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2,v=3,计算:(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u)ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4.计算如下表达式:(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。

matlab 教程 第三章Matlab 绘图与例题

matlab 教程 第三章Matlab 绘图与例题
例1、已知向量x=[1 2 3],y=[4 7 9 0],生成它们对应的 格点矩阵。注意:输出的X、Y都是4*3矩阵,X的行 向量都是向量x,Y的列向量都是向量y
机动 目录 上页 下页 返回 结束
二、三维网格图命令mesh
• mesh(X,Y,Z) 生成网格曲面,X,Y,Z是同维数的矩阵 • mesh(x,y,Z) x,y是向量,而Z是矩阵。等价于
机动 目录 上页 下页 返回 结束
三、三维表面图命令surf
• surf的调用格式与mesh相同,不同之处是surf绘的是曲面 而不是网格。
• 三维表面图可以用shading命令修饰其显式形式。
Matlab还提供了waterfall(x,y,z)及contour3(x,y,z)等命 令绘制三维图形。
3
2
在[0,4π]间的图形。
解:syms t

ezplot(‘2/3*exp(-t/2)*cos(3/2*t)’,[0,4*pi])
第二节 Matlab二维特殊图形
Matlab提供了许多其他的二维绘图指令, 大大扩充了Matlab的曲线作图指令,可以满足 用户的不同需要。
注 目录 上页 下页 返回 结束
axis equal 使坐标轴在三个方向上刻度增量相同
axis square 使坐标轴在三个方向上长度相同
axis
返回表示当前图形坐标轴的范围
机动 目录 上页 下页 返回 结束
第四节 Matlab空间曲面绘图
二元函数 z f (x, y) 的图形是三维空间曲面,函
数图形在了解二元函数的特性上帮助很大。
第3章 Matlab 绘图与例题
第一节
第3章
Matlab二维曲线绘图
一、基本绘图指令plot

MATLAB第三章

MATLAB第三章

第三章微积分问题的计算机求解一、实验内容:题目1.试求出如下极限。

①limx→∞(3x +9x )1/ x,②lim x→∞[(x+2)x+2(x+3)x+3 ]/(x+5)2x+5【分析】:该题为单变量函数的极限。

极限问题可以用limit()函数直接求出。

要注意该函数的调用格式为:L=limit(fun,x,x0)(求极限),L=limit(fun,x,x0,’left’或’right’)(求极限)。

还需注意一开始要对函数的字符进行申明。

【解答】:(1)输入如下语句:>> syms x;f=(3^x+9^x)^(1/x);L=limit(f,x,inf)语句运行后显示如下:L =9(2)输入如下语句:>>syms x;f=(x+2)^(x+2)*(x+3)^(x+3)/(x+5)^(2*x+5);>> L=limit(f,x,inf)语句运行后显示如下:L =exp(-5)题目2.试求下面的双重极限。

①lim x→−1y→2 (x2y+xy3)/(x+y) 3,②limx→0 y→0 xy /√(xy+1)−1,③limx→0y→0 [1−cos(x2+y2)]/(x2+y2)e x2+y2。

【分析】:该题为多变量函数的极限问题。

他可以用嵌套使用limit()函数来解决。

在MATLAB上可以用L=limit(limit(f,x,x0),y,y0)或者L=limit(f,y,y0),x,x0)来解决。

其思想是所有的先关于X求导,再所有的关于y求导。

【解答】:(1)输入如下语句:>> syms x y>> f=(x^2*y+x*y^3)/(x+y)^3;>> L=limit(limit(f,x,-1),y,2)语句运行后显示如下:L =-6(2)输入如下语句:>> syms x yf=(x*y)/(sqrt(x*y+1)-1);L=limit(limit(f,x,0),y,0)按ENTER键,语句运行后显示如下:L =2(3)输入如下语句:>> syms x yf=(1-cos(x^2+y^2))/(sqrt(x^2+y^2)*exp(x^2+y^2));L=limit(limit(f,x,0),y,0)按ENTER键,语句运行后显示如下:L =题目3.求出下面函数的导数。

MATLAB语言与控制系统仿真_参考题答案_第3章

MATLAB语言与控制系统仿真_参考题答案_第3章

3.5 MATLAB 绘图实训3.5.1 实训目的1.学会MATLAB 绘图的基本知识;2.掌握MATLAB 子图绘制、图形注释、图形编辑等基本方法;3.学会通过MATLAB 绘图解决一些实际问题;4.练习二维、三维绘图的多种绘图方式,了解图形的修饰方法;5.学会制作简单的MATLAB 动画。

图3-46 炮弹发射示意图3.5.2 实训内容1. 炮弹发射问题〔1炮弹发射的基础知识炮弹以角度α射出的行程是时间的函数,可以分解为水平距离)(t x 和垂直距离)(t y 。

)cos()(0αtv t x = %水平方向的行程; 205.0)sin()(gt tv t y -=α %垂直方向的行程;其中,0v 是初速度;g 是重力加速度,为9.82m/s ;t 是时间。

〔2炮弹发射程序举例:分析以下程序以及图3-47各个图形的实际意义。

a=pi/4; v0=300; g=9.8;t=0:0.01:50; x=t*v0*cos<a>;y=t*v0*sin<a>-0.5*g*t.^2;subplot<221>;plot<t,x>;grid;title<‘时间-水平位移曲线'>; subplot<222>;plot<t,y>;grid;title<‘时间-垂直位移曲线'>; subplot<223>;plot<x,y>;grid;title<‘水平位移-垂直位移曲线'>; subplot<224>;plot<y,x>;grid;title<‘垂直位移-水平位移曲线'>; 图3-4745角发射曲线 〔3编程解决炮弹发射问题①假设在水平地面上以垂直于水平面的角度向上发射炮弹,即发射角90=α,假设初速度分别为[310,290,270]m/s,试绘制时间-垂直位移曲线,编程求取最高射程;绘图要求:◆ 标题设为"炮弹垂直发射问题";◆ 在图上通过添加文本的方式表明初速度; ◆ 在x 轴标注"时间";◆ 在y 轴上标注"垂直距离"; ◆ 添加网格线;◆ 将310m/s 的曲线改为线粗为2的红色实线; ◆ 将290m/s 的曲线改为线粗为3的绿色点划线;◆ 将270m/s 的曲线改为线粗为2的蓝色长点划线;a=pi/2; v1=310; g=9.8;t=0:0.01:50; x1=t*v1*cos<a>;y1=t*v1*sin<a>-0.5*g*t.^2;plot<t,y1>;grid; title<'炮弹垂直发射问题'>; xlabel<'时间'>; ylabel<'垂直距离'>; hold on; v2=290;x2=t*v2*cos<a>;y2=t*v2*sin<a>-0.5*g*t.^2; plot<t,y2>; v3=270;x3=t*v3*cos<a>;y3=t*v3*sin<a>-0.5*g*t.^2; plot<t,y3>;zgsc=[max<y1>; max<y2>; max<y3>] %三次发射的最高射程 运行结果如下: zgsc =1.0e+003 * 4.9031 4.29083.7194最高射程分别为:4903.1米,4290.8米,3719.4米。

matlab第二版课后习题答案

matlab第二版课后习题答案

matlab第二版课后习题答案
《MATLAB第二版课后习题答案》
MATLAB是一种强大的数学软件,被广泛应用于工程、科学和金融等领域。

《MATLAB第二版》是一本经典的教材,为了帮助学生更好地掌握MATLAB的使用,书中提供了大量的课后习题。

下面我们将为大家总结一些MATLAB第二版课后习题的答案,希望能对大家的学习有所帮助。

1. 第一章课后习题答案
第一章主要介绍了MATLAB的基本操作,包括变量的定义、矩阵的运算、函数的使用等。

在课后习题中,有一道题目是要求计算一个矩阵的逆矩阵。

答案是使用MATLAB中的inv函数,将原矩阵作为参数传入即可得到逆矩阵。

2. 第二章课后习题答案
第二章介绍了MATLAB中的绘图功能,包括二维和三维图形的绘制。

有一道课后习题是要求绘制一个正弦曲线和余弦曲线,并在同一张图上显示。

答案是使用MATLAB中的plot函数,分别绘制正弦曲线和余弦曲线,并使用legend函数添加图例。

3. 第三章课后习题答案
第三章介绍了MATLAB中的控制流程,包括if语句、for循环和while循环等。

有一道课后习题是要求编写一个程序,计算1到100之间所有偶数的和。

答案是使用for循环遍历1到100之间的所有数,判断是否为偶数并累加。

通过以上几个例子,我们可以看到MATLAB第二版课后习题的答案涵盖了各种基本和高级的操作,对于学习MATLAB是非常有帮助的。

希望大家在学习MATLAB的过程中能够多加练习,掌握更多的技巧和方法。

matlab1-8章课后作业

matlab1-8章课后作业

MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件,是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境,已被广泛应用于不同领域。

Matlab的基本功能包括:数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。

Matlab第二章1.创建double的变量,并进行计算。

(1)a=87,b=190,计算a+b、a-b、a*b。

(2)创建uint8 类型的变量,数值与(1)中相同,进行相同的计算。

>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =ans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2,v=3,计算:(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans =15.4189>> sqrt(u-3*v)/(u*v) ans =0 + 0.4410i 4.计算如下表达式:(1)(2)>> (3-5*i)*(4+2*i)22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。

第三章matlab矩阵运算

第三章matlab矩阵运算
例3-22 计算向量a=-4:2:6每个元素的符号。 a=-4:2:6; B=sign(a)
3.2.3 坐标变换函数(P52)
例3-23 将迪卡尔坐标系中(1,1,1)分别转换到球坐 标系和极坐标中。 [THETA,PHI,R]=cart2sph(1,1,1) P= [THETA,PHI,R] [THETA,PHI,Z]=cart2pol(1,1,1) Q= [THETA,PHI,Z] R=[P;Q]
ankn
例3-3 求向量eye(4),magic(4)和A=[1,2,3 ;4,5,6;7,8,9]的行列式。 det(eye(4)); det(magic(4)); det(A);
4.矩阵的行列迹: 矩阵的迹定义为对角元素之和。Matlab中用函数trace( )来计算矩阵的行列式。 例3-4 求向量eye(4),magic(4)和A=[1,2,3 ;4,5,6;7,8,9]的行列式。 trace(eye(4)); trace(magic(4)); trace(A);
2.LU分解:
LU分解是将任意一个方正A分解成为一个交换下三角矩阵 L(或是排列(permuted) 的上三角形矩阵)和一个上三角矩 阵U的乘积,A=LU,在Matlab中用函数lu来计算LU分解
例3-14 求矩阵A=[1,4,2;5,6,9;4,1,8]的LU分解,
[L1,U1]=lu(A)
L1*U1
R=rref(A2)
9.矩阵空间之间的角度:
矩阵空间之间的角度代表具有相同行数的两个矩阵线性 相关程度,夹角越小代表线性相关度越高。Matlab中用函 数subspace()来计算矩阵空间之间的角度。
例3-9 求矩阵A1=[1,2,3;4,5,6;7,8,9]和A2=[1,2;3 ,4;5,6]之间的夹角Q。

matlab程序设计第三章课后习题答案资料

matlab程序设计第三章课后习题答案资料

m a t l a b程序设计第三章课后习题答案1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。

>> t=0:0.01:pi;>> x1=t;>> y1=2*x1-0.5;>> x2=sin(3*t).*cos(t);>> y2=sin(3*t).*sin(t);>> plot(x1,y1,'r-',x2,y2,'g-')>> axis([-1,2,-1.5,1])>> hold on>> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)');>> plot(double(s.x),double(s.y),'*');截图:p366 第4题绘制极坐标曲线,并分析对曲线形状的影响。

function [ output_args ] = Untitled2( input_args )%UNTITLED2 Summary of this function goes here% Detailed explanation goes heretheta=0:0.01:2*pi;a=input('请输入a的值:');b=input('请输入b的值:');n=input('请输入n的值:');rho=a*sin(b+n*theta);polar(theta,rho,'k');end下面以a=1,b=1,n=1的极坐标图形为基础来分析a、b、n的影响。

对a的值进行改变:对比发现a只影响半径值的整倍变化对b的值进行改变:对比发现b的值使这个圆转换了一定的角度对n的值进行改变:对比发现当n>=2时有如下规律1、当n为整数时,图形变为2n个花瓣状的图形2、当n为奇数时,图形变为n个花瓣状的图形分别让n为2、3、4、5同时改变b和n的值,对比发现图形变为2n个花瓣状,并且还旋转了一定角度1 绘制sin(x)曲线,通过句柄修改曲线的颜色,要求:a)生成曲线的同时,获取句柄,再修改颜色;>> x=0:pi/50:2*pi;>> y=sin(x);>> h=plot(x,y)h =174.0101>> set(h,'color','r');截图如下:b)先生成曲线,然后再获取句柄,再完成颜色修改。

Matlab第三章习题答案

Matlab第三章习题答案

Matlab第三章习题答案第三章3.5计算多项式乘法()()222254x x x x ++++>> a=[1 2 2];>> b=[1 5 4];>> c=conv(a,b)c =1 7 16 18 8poly2sym(c)ans =x^4 + 7*x^3 + 16*x^2 + 18*x + 83.6计算多项式的除法()()3231368/4x x x x ++++ >> a=[3 13 6 8];>> b=[1 4];>> [q r]=deconv(a,b)q =3 1 2r =0 0 0 0>> poly2sym(q)ans =3*x^2 + x + 23.8求多项式4324121459x x x x --++的微分和积分(1)微分:>> a=[4 -12 -14 5 9];>> b=polyder(a)b =16 -36 -28 5>> poly2sym(b)ans =16*x^3 - 36*x^2 - 28*x + 5(2)积分:polyint(a)ans =0.8000 -3.0000 -4.6667 2.5000 9.00000 >> poly2sym(ans)ans =(4*x^5)/5 - 3*x^4 - (14*x^3)/3 + (5*x^2)/2 + 9*x3.9求代数⽅程32349x x x +-+=0的解。

>> [x]=solve('x^3+3*x^2-4*x+9')x =- 7/(3*(15/2 - (108^(1/2)*4703^(1/2))/108)^(1/3)) - (15/2 -1/108*108^(1/2)*4703^(1/2))^(1/3) - 17/(6*(15/2 - (108^(1/2)*4703^(1/2))/108)^(1/3)) + (15/2 -(108^(1/2)*4703^(1/2))/108)^(1/3)/2 - (3^(1/2)*i*(7/(3*(15/2 -(108^(1/2)*4703^(1/2))/108)^(1/3)) - (15/2 - 1/108*108^(1/2)*4703^(1/2))^(1/3)))/2 - 1 7/(6*(15/2 - (108^(1/2)*4703^(1/2))/108)^(1/3)) + (15/2 -(108^(1/2)*4703^(1/2))/108)^(1/3)/2 + (3^(1/2)*i*(7/(3*(15/2 -(108^(1/2)*4703^(1/2))/108)^(1/3)) - (15/2 - 1/108*108^(1/2)*4703^(1/2))^(1/3)))/2 – 1 3.10求线性⽅程组223430x xy y x x ?++=??-+=??的解 eqn1='x^2+x*y+y=3';>> eqn2='x^2-4*x+3=0';[x,y]=solve(eqn1,eqn2)x =13y =1-3/23.11求微分⽅程64dy y x dx=+的通解。

MATLABSimulink电力系统建模与仿真第三章例题主讲教师

MATLABSimulink电力系统建模与仿真第三章例题主讲教师

MATLAB电力变压器模块
• 选用的各模块名称及提取路径。
模 块 名 三相双绕组变压器 三相可编程交流电源 接地模块 Ground 中性点模块 node 10 万用表模块 Multimeter 信号分离模块 Demux 示波器 Scope 提 取 路 径 SimPowerSystems/Elements SimPowerSystems/Electrical Sources SimPowerSystems/Elements SimPowerSystems/Elements SimPowerSystems/Measurements Simulink/ Signal Routing Simulink/Sinks
Lm*
Lm V12 n Pn 2πf n
1.674
MATLAB电力变压器模块
• 双击万用表模块打开万用表参数设置窗口。 • 选择对变压器上所有变量进行测量,因此可测量的参数包括各绕组端 电压、电流、励磁电流和磁通,选择变压器一次侧的a相电压、a相主 磁通和a相电流为测量变量。
MATLAB电力变压器模块
• 打开菜单Simulation>Configuration Parameters,在“算法选择 ”(Solver options)窗口中选择“变步长”(Variable -step)和“刚性 积分算法(ode15s)”。
MATLAB电力变压器模块
• 仿真及结果 • 图中波形从上至下分别为一次侧相电压、空载电流和主磁通。 为了进行比较,在仿真得到的各个波形上叠加了理想正弦波。 • 空载相电流为正弦波,主磁 通发生了很小的畸变但仍近似 为正弦波。由于磁滞损耗的存 在,主磁通滞后空载相电流一 个铁耗角,由主磁通产生的感 应电动势滞后主磁通90°且波 形发生畸变。此时,尽管电压 源电压波形为理想正弦波,但 变压器一次侧空载相电压并不 是理想正弦波,而只是近似正 弦波。 • 仿真结果与理论分析一致。

MATLAB程序设计知到章节答案智慧树2023年广西科技大学

MATLAB程序设计知到章节答案智慧树2023年广西科技大学

MATLAB程序设计知到章节测试答案智慧树2023年最新广西科技大学第一章测试1.当在命令行窗口执行命令时,如果不想立即在命令行窗口中输出结果,可以在命令后加上()。

参考答案:分号(;)2.MATLAB命令行窗口中提示用户输入命令的符号是()。

参考答案:>> ;3.fix(264/100)+mod(264,10)*10的值是()。

参考答案:424.下列可作为MATLAB合法变量名的是()。

参考答案:xyz_2a。

5.在命令行窗口输入下列命令后,x的值是()。

>> clear>> x=i*j参考答案:-1;6.建立矩阵时,不同行的元素之间用()分隔。

参考答案:分号;7.输入字符串时,要用()将字符序列括起来。

参考答案:' ';8.最初的MATLAB核心程序是采用()语言编写的。

参考答案:FORTRAN;9.下列选项中能反应MATLAB特点的是()。

参考答案:编程效率高。

10.如果要重新执行以前输入的命令,可以使用()。

参考答案:上移光标键(↑)。

11.plot(x,y)是一条()命令。

参考答案:绘图;12.以下两个命令行的区别是()。

>> x=5,y=x+10>> x=5,y=x+10;参考答案:第一个命令行同时显示x和y的值,第二个命令行只显示x的值;13.下列选项中,不是MATLAB帮助命令的是()。

参考答案:search;14.下列数值数据表示中错误的是()。

参考答案:2e;15.下列语句中错误的是()。

参考答案:x=y=3;16.已知s='显示"hello"',则s的元素个数是()。

参考答案:9;17.建立矩阵时,同一行的元素之间用()分隔。

参考答案:空格;;逗号;18.在当前文件夹和搜索路径中都有fpp.m文件,那么在命令行窗口输入fpp时,下列说法错误的是()。

参考答案:先执行当前文件夹的fpp.m文件,再执行搜索路径中的fpp.m文件;;先执行搜索路径中的fpp.m文件,再执行当前文件夹的fpp.m文件;;执行搜索路径中的fpp.m文件;19.内存变量文件是二进制格式文件,扩展名为.mat。

matlab第三章答案

matlab第三章答案

第 1 题: 从键盘输入一个3位数,将它反向输出。

如输入639 输出936、n1=fix(n/100); n2=rem(fix(n/10),10);n3=rem(n,10); m=n3*100+n2*10+n1;disp(m)请输入一个三位数:489984第2题:输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。

其中90~100分为A,80~89分为B、70~79 分为C、60~69 分为D,60 分以下为E。

要求:(1) 分别用if 语句与switch 语句实现。

(2)输入百分制成绩后要判断成绩的合理性,对不合理的成绩应输出错信息。

If 语句:>> s=input(' 请输入成绩分数:');if s>=90&s<=100disp('A');elseif s>=80&s<=89disp('B');elseif s>=70&s<=79disp('C');elseif s>=60&s<=69disp('D');elseif s>=0&s<60disp('E');else disp('false');end请输入成绩分数:89BSwitch 语句:>> s=input(' 请输入成绩分数:');switch fix(s/10)case {9,10}disp('A');case {8}disp('B');case {7}disp('C');case {6}disp('D');case {5,4,3,2,1,0}disp('E');otherwise disp('false');end请输入成绩分数:88B第3题输入20个数,求其中最大数与最小数。

MATLAB第三章数值数组及其运算

MATLAB第三章数值数组及其运算

行向量
如:array=[2, pi/2, sqrt(3), 3+5i]
x=[1,2,3,4,5都已知.如对 少量实验数据的处理可用此种方法.
4
(2) 冒号生成法: array=a: inc: b
<向量>
a---数组的第一个元素
inc---采样点之间的间隔, 即步长. 最后一个元素不一定等于b, 其大小为b’=a +inc*[(b-a)/inc]; 步长可以省略, 默认为 1; inc可以取正数或负数, 但要注意当取正时,要保证b>a, 数 组最后一个元素不超过b, 取负时b<a, 最后一个元素不小于b.
(2) 数值计算解法
delt=0.01; x=0:delt:4;
y=exp(-sin(x));
sx=delt*cumtrapz(y);
plot(x,y, 'r', 'LineWidth', 6); hold on;
plot(x, sx, '.b', 'MarkerSize', 15);
plot(x, ones(size(x)), 'k');
a inc>0 b
b inc>0 a
特点: 等差数列
方便对数据之间的间隔(步长)进行控制.但要注意三个数值之 间的关系,可能得到空数组.另外要注意生成的数组的元素的 个数.如x=a: (b-a)/n :b (b>a)得到n+1个元素的数组.
5
x=1:5x=[1,2,3,4,5]
y=5:-1:1y=[5, 4, 3, 2, 1]
2. 在命令窗中输入MyMatrix
11
3.5 二维数组的标识 (mxn, m>1, n>1)

MATLAB第三章和第四章作业

MATLAB第三章和第四章作业

习题三1.选择题(1)已知a=4,b='4',下面说法错误的为C 。

A. 变量a比b占用的存储空间大B. 变量a和b可以进行加、减、乘、除运算C. 变量a和b的数据类型相同D. 变量b可以用eval命令执行(2)已知s='显示"hello"',则s的元素个数是B 。

A. 12B. 9C. 7D. 18(3)运行命令>>a=sym('pi','d')则对于变量a的描述D是正确的。

A. a是符号变量B. a显示为10位的数值C. a显示为32位的数值D. a不存在(4)对于符号表达式g=sym('sin(a*z)+cos(w*v)'),如果运行命令diff(g),则求导自变量是D。

A. aB. zC. wD. v(5)运行命令>>a=double(sym('sin(pi/2)')),则变量a是 C。

A. 符号变量B. 字符串'1'C. double型的数值1D. 出错(6)运行命令>>y=dsolve('x*D2y-3*Dy=x^2','t')求解微分方程,则D。

A. Dy是指dy/dxB. 得出的y是通解有一个常数符C1C. D2y是指d2y/dxD. 得出的y是通解有两个常数符C1和C2(7)有一个2行2列的元胞数组c,则c(2)是指C 。

A. 第1行第2列的元素内容B. 第2行第1列的元素内容C. 第1行第2列的元素D. 第2行第1列的元素(8)对于一个2行2列的元胞数组c,如果要删除第二列的元素,应该使用的命令是B 。

A.c{:,2}=[]B.c(:,2)=[]C.A、B两种方式均可D.A、B两种方式均不可(9)有一个3行10列的构架数组student,每个构架有name(姓名)、scores两个字段,其中scores是以1×5数组表示的五门课的成绩,那么要查看第4个学生的第2门课成绩,采用的命令中达不到要求的是。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

%T3.3
%绘制一条半径为2的
园,要求在图形窗口中
显示是圆形
theta=0:0.1:2*pi;
r1=sin(theta);
polar(theta,r1,'r')
title('T3.3')
%T3.4
%绘制曲线y=5*t*sin(2*pi*t),t的范围是0:2
t=0:0.1:2;
y=5*t*sin(2*pi*t);
plot(t,y) %T3.6
%在同一图形窗口分别绘制y1=x y2=x^2 y3=e^-x x 范围【-2 6】给整个图形加标题横坐标加标注右上角标注三条曲线的图例文字标注x=1点,并在x=1处绘制一条【-2 10】的垂直线
x=-2:0.1:6;
y1=x ;
plot(x,y1)
hold on
y2=x.^2 ;
plot(x,y2)
y3=exp(-x);
plot(x,y3)
xlabel('x')
legend('y1=x','y2=x^2','y3=e^-x',1)%1:右上角
text(1,0,'x=1')%在(1,0)处写x=1
line(1,-2:.1:10)
title('T3.6')
hold
%T3.5
%在同一窗口绘制y1=sin(t) t=0:2*pi 黑色点画线,y2=sin(2*t) t=pi:4*pi 红色许仙圆圈
t=0:0.1:2*pi;
y1=sin(t);
plot(t,y1,'k-.')
hold on;
t=pi:0.1:4*pi;
y2=sin(2*t);
plot(t,y2,'r--')
gtext({'y1=sin(t)','y2=sin(2*t)'}) title('T3.5')
Hold %3.9
%用semillogx命令绘制传递函数为10/(s*(.5*s+1)) 的对数幅频特性曲线,横坐标为w纵坐标为Lw,w范围是10.^-2 ~10.^3按对数分布,兵绘制双对数坐标图
w=logspace(-2,3,10);
Lw=-20*log10(w.*sqrt((0.5*w). ^2+1));
semilogx(w,Lw)
loglog(w,Lw)
title('T3.9')。

相关文档
最新文档