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. 12

B. 9

C. 7

D. 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. a

B. z

C. w

D. v

(5)运行命令>>a=double(sym('sin(pi/2)'))，则变量a是 C。

A. 符号变量

B. 字符串'1'

C. double型的数值1

D. 出错

(6)运行命令>>y=dsolve('x*D2y-3*Dy=x^2','t')求解微分方程，则D。

A. Dy是指dy/dx

B. 得出的y是通解有一个常数符C1

C. D2y是指d2y/dx

D. 得出的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门课成绩，采用的命令中达不到要求的是。

A.student(1,2).scores(2)

B.student(4).scores(2)

C.两种方式均可

D.两种方式均不可

(10)构架数组结构同上，如果增添一个字段ave用来存放每个同学的平均成绩，则要查看所有同学的平均成绩，应该采用的命令是D。

A.[student.ave]

B.{student.ave}

C.两种方式均可

D.两种方式均不可

(11)运行以下命令：

>>x=[1 2 3;4 5 6];

>>y=x+x*i

>>plot(y)

则在图形窗口中绘制A条曲线。

A. 3

B. 2

C. 6

D. 4

(12)运行以下命令：

>>x=[1 2 3;4 5 6];

>>plot(x,x,x,2*x)

则在图形窗口中绘制B条曲线。

A. 4

B. 6

C. 3

D. 5

(13)subplot(2,1,1)是指A的子图。

A. 两行一列的上图

B. 两行一列的下图

C. 两列一行的左图

D. 两列一行的右图

(14)运行命令>>figure(3)，则执行D 。

A. 打开三个图形窗口

B. 打开一个图形窗口

C. 打开图形文件名为“3.fig”

D. 打开图形文件名为“figure 3.fig”

(15)运行以下命令：

>>x=0:0.1:2*pi;

>>y=sin(x);

>>plot(x,y)

则如果要使正弦曲线充满坐标轴则以下A命令不能使用。

A. axis image

B. axis([0,2*pi,-1,1])

C. axis fill

D. axis tight

(16)如果要显示向量中各元素占和的百分比，则使用B 函数绘图。

A. hist

B. pie

C. bar

D. stairs

(17)极坐标图是使用B来绘制的。

A. 原点和半径

B. 相角和距离

C. 纵横坐标值

D. 实部和虚部

(18)meshc函数是D 。

A. 绘制三维曲线图

B. 绘制三维网线图并添加平行于z轴的边框线

C. 绘制三维表面图

D. 绘制三维网线图并添加等高线

(19)三维图形中的默认视角是C。

A. 方位角＝0°俯仰角＝90°

B. 方位角＝90°俯仰角＝0°

C. 方位角＝37.5°俯仰角＝30°

D. 方位角＝0°俯仰角＝180°

(20)二维图形中的colorbar命令运行后，颜色条显示D。

A. 无色

B. 黑色

C. 白色

D. 有颜色但无意义

2.用构架数组来存储5名学生的基本情况数据，每名学生的数据包括学号、姓名、专业和6门课程的成绩。

解：建立学生构架数组过程如下：

student(1).studno='001';

student(1).name='张三';

student(1).major='机械';

student(1).score=[75 82 95 90 87 83];

student

student.score

student =

studno: '001'

name: '张三'

major: '机械'

score: [75 82 95 90 87 83]

ans =

75 82 95 90 87 83

student(5).name='李四';

student

student =

1x5 struct array with fields:

studno

name

major

score

3.建立单元矩阵B并回答有关问题。

B{1,1}＝1；

B{1,2}＝’Brenden’；

B{2,1}＝reshape(1:9，3，3)；

B{2,2}={12，34，2；54，21，3；4，23，67}；

(1)size(B)和ndims(B)的值分别是多少?

(2)B(2)和B(4)的值分别是多少?

(3)B(3)＝[]和B{3}＝[]执行后，B的值分别是多少?

解：建立单元矩阵B的过程如下：

B{1,1}=1;

B{1,2}='Brenden';

B{2,1}=reshape(1:9,3,3);

B{2,2}={12,34,2;54,21,3;4,23,67};

这是一个2×2的单元数组，其中第2行第2列的单元又是一个单元数组，因此是一个嵌套的单元数组。

size(B)

ndims(B)

B(2),B(4)

ans =

2 2

ans =

2

ans =

[3x3 double]

ans =

{3x3 cell}

B(3)=[];B

B =

[1] [3x3 double] {3x3 cell}

此时B变为1×3的单元数组。

B{3}=[];B

B =

[1] [3x3 double] []

此时B仍然是1×3的单元数组，但第3个单元的内容变为空数组。

4.确定下面各符号表达式中的自由符号变量：

1/(log(t)+log10(w*t)) sqrt(t)/y 10*i+x*j exp(-a*result)

解：这四个符号表达式中的自由符号变量分别为：第一个：w、t；第二个：y、t；第三个：x （注意：i和j作为虚数单位看待，不能作为自由符号变量）；第四个：result、a。也可以采用函数findsym来获取自由符号变量：

syms t w x y a result

y1=1/(log(t)+log10(w*t))

y2=sqrt(t)/y

y3=10*i+x*j

y4=exp(-a*result)

findsym(y1)

findsym(y2)

findsym(y3)

findsym(y4)

y1 =

1/(log(t)+log(w*t)/log(10))

y2 =

t^(1/2)/y

y3 =

10*i+i*x

y4 =

exp(-a*result)

ans =

t, w

ans =

t, y ans = x

ans =

a, result

5.对符号表达式cos f x =+，分别使用collect 、expand 和simplify 函数化简，并与

simple 函数的结果比较。

解：运算如下： syms x;

f=cos(x)+sqrt(-sin(x)^2); collect(f) ans =

cos(x)+(-sin(x)^2)^(1/2) expand(f) ans =

cos(x)+(-sin(x)^2)^(1/2) simplify(f) ans =

cos(x)+(-sin(x)^2)^(1/2) simple(f) simplify:

cos(x)+(-sin(x)^2)^(1/2) radsimp:

cos(x)+i*sin(x) combine(trig):

cos(x)+1/2*(-2+2*cos(2*x))^(1/2) factor:

cos(x)+(-sin(x)^2)^(1/2) expand:

cos(x)+(-sin(x)^2)^(1/2) combine:

cos(x)+1/2*(-2+2*cos(2*x))^(1/2) convert(exp):

1/2*exp(i*x)+1/2/exp(i*x)+1/2*((exp(i*x)-1/exp(i*x))^2)^(1/2) convert(sincos):

cos(x)+(-sin(x)^2)^(1/2) convert(tan):

(1-tan(1/2*x)^2)/(1+tan(1/2*x)^2)+(-4*tan(1/2*x)^2/(1+tan(1/2*x)^2)^2)^(1/2)

collect(x):

cos(x)+(-sin(x)^2)^(1/2) mwcos2sin:

cos(x)+(-sin(x)^2)^(1/2) ans =

cos(x)+i*sin(x)

6.已知符号表达式3

2

541f x x x =+++，x

g e -=，求复合函数(())f g x ，并将f 转换为多

项式系数。

解：运算如下： syms x;

f=x^3+5*x^2+4*x+1;

g=exp(-x);

h=compose(f,g,x)

simplify(h)

h =

exp(-x)^3+5*exp(-x)^2+4*exp(-x)+1

ans =

exp(-3*x)+5*exp(-2*x)+4*exp(-x)+1

c=sym2poly(f)

c =

1 5 4 1

7.将符号表达式21

=-中的x-1用a或5替换，并求y的反函数。

y x

解：本题首先应将表达式y进行因式分解，然后再进行相应的替换。运算如下：

syms x;

y=x^2-1;

y=factor(y)

y =

(x-1)*(x+1)

subs(y,'x-1','a')

subs(y,'x-1',5)

ans =

((a))*(x+1)

ans =

5*x+5

g=finverse(y)

Warning: finverse((x-1)*(x+1)) is not unique.

> In sym.finverse at 48

g =

(x+1)^(1/2)

8.对符号表达式2sin(/4)

=+求t趋向1时的极限，并使用级数和求前10项。

y t tπ

解：运算如下：

syms t;

y=2*t*sin(t+pi/4);

limit(y,t,1)

ans =

sin(1)*2^(1/2)+cos(1)*2^(1/2)

y1=taylor(y,t,10,1)

y1 =

2*sin(1+1/4*pi)+(2*cos(1+1/4*pi)+2*sin(1+1/4*pi))*(t-1)+(-

sin(1+1/4*pi)+2*cos(1+1/4*pi))*(t-1)^2+(-sin(1+1/4*pi)-

1/3*cos(1+1/4*pi))*(t-1)^3+(-1/3*cos(1+1/4*pi)+1/12*sin(1+1/4*pi))*(t-1)^4+(1/60*cos(1+1/4*pi)+1/12*sin(1+1/4*pi))*(t-1)^5+(-

1/360*sin(1+1/4*pi)+1/60*cos(1+1/4*pi))*(t-1)^6+(-1/360*sin(1+1/4*pi)-1/2520*cos(1+1/4*pi))*(t-1)^7+(-1/2520*cos(1+1/4*pi)+1/20160*sin(1+1/4*pi))*(t-1)^8+(1/181440*cos(1+1/4*pi)+1/20160*sin(1+1/4*pi))*(t-1)^9

9.求解符号方程组123412412341

234232836

877225

x x x x x x x x x x x x x x x -++=??

++=??-++=??+-+=?

解：首先解系数矩阵和右端向量定义为符号对象，然后求解。运算如下：

A=sym([2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2])

b=sym([8;6;7;5]) X=A\b

A =

[ 2, -3, 1, 2] [ 1, 3, 0, 1] [ 1, -1, 1, 8] [ 7, 1, -2, 2] b = 8 6 7 5 X =

611/241 290/241 1646/241 -35/241

或者直接采用函数solve 进行求解：

X=solve('2*x1-3*x2+x3+2*x4=8','x1+3*x2+x4=6','x1-x2+x3+8*x4=7','7*x1+x2-2*x3+2*x4=5','x1','x2','x3','x4')

X =

x1: [1x1 sym] x2: [1x1 sym] x3: [1x1 sym] x4: [1x1 sym] 查看方程组的解：

X.x1,X.x2,X.x3,X.x4 ans = 611/241 ans = 290/241 ans = 1646/241 ans = -35/241

10.求符号微分方程tan cos dy

y x x dx

+=的通解和当y(0)=2的特解。 解：运算如下：

eqn='Dy+y*tan(x)=cos(x)'; 求通解：

dsolve(eqn)

ans =

exp(-tan(x)*t)*cos(x)^2/sin(x)*exp(sin(x)/cos(x)*t)+exp(-

tan(x)*t)*C1

求特解：

dsolve(eqn,'y(0)=2')

ans =

exp(-tan(x)*t)*cos(x)^2/sin(x)*exp(sin(x)/cos(x)*t)+exp(-

tan(x)*t)*(2-1/sin(x)*cos(x)^2)

11.在0～10的坐标轴范围内绘制三条曲线，一条水平线，一条垂直线，一条对角斜线。

解：编程如下：

t=0:10;

l=length(t);

x1=t;

x2=5*ones(l);

plot(t,x1,t,x2,'r',x2,x1,'b')

12.绘制一条半径为2的圆，要求在图形窗口中显示的是圆形。

解：编程如下：

t=-pi:0.01*pi:pi;

x=2*cos(t);y=2*sin(t);

plot(x,y)

axis square

13.绘制函数曲线5sin(2)y t t π=，t 的范围为0~2。 解：编程如下： t=0:0.01:2;

y=5*t.*sin(2*pi*t); plot(t,y)

14.在同一图形窗口绘制曲线1sin()y t =，t 的范围为0~2π，2sin(2)y t =，t 的范围为π～4π；要求1y 曲线为黑色点划线，2y 曲线为红色虚线圆圈，使用鼠标将文字标注添加到两条曲线上。

解：编程如下：

t=pi:0.01*pi:4*pi; y1=sin(t);

y2=sin(2*t);

plot(t,y1,'k-.',t,y2,'r--o')

15.在同一图形窗口分别绘制1y x =、2

2y x =、3x

y e -=三条函数曲线，x 的范围为[-2 6]，要

求给整个图形加上标题，给横坐标轴加上标注，图的右上角标注三条曲线的图例，使用文字标注

x=1点，并在x=1处绘制一条[-2,10]的垂直线。

解：编程如下： x=-2:0.01:6; y1=x; y2=x.^2;

y3=exp(-x);

plot(x,y1,x,y2,x,y3) xlabel('X 轴') ylabel('Y 轴')

title('三条函数曲线')

legend('y_1=x','y_2=x^2','y_3=e^{-x}') text(1,0,'x=1') hold on

plot([1,1],[-2,10])

16.已知某班10个同学的成绩为65、98、68、75、88、82、94、85、56，分别统计并绘制60分以下、60-70、70-80、80-90、90-100分数段的人数图；并使用饼形图显示各分数段所占的百分比。

解：编程如下：

score=[65,98,68,75,88,82,94,85,56];

x=[55,65,75,85,95];

hist(score,5)

N=hist(score,x)

N =

1 2 1 3 2

17.已知某班5个同学的三次成绩为

6578869369

7585929570

7280799272

??

??

??

??

??

，使用柱状图和阶梯图显示每个

同学的成绩变化。

解：编程如下：

score=[65,78,86,93,69;75,85,92,95,70;72,80,79,92,72]';

bar(score)

stairs(score)

18.用semilogx 命令绘制传递函数为

10

(0.51)

s s +的对数幅频特性曲线，横坐标为w 纵坐标为

Lw ，w 范围为10-2

~103

按对数分布，并绘制双对数坐标图。

解：传递函数()G s 的对数幅频特性为()20lg ()L G j ωω=，因此首先要利用符号运算求出

()L ω，然后才能绘制图形。具体步骤如下：

syms s;

G=10/(s*(0.5*s+1))

G=subs(G,'s','j*omega')

G =

10/s/(1/2*s+1) G =

10/(j*omega)/(1/2*(j*omega)+1) L=20*log10(abs(G))

L =

20*log(10*abs(1/j/omega/(1/2*j*omega+1)))/log(10) w=logspace(-2,3,10);

Lw=20*log(10*abs(1/j./w./(1/2*j.*w+1)))/log(10); semilogx(w,Lw)

19.绘制sin(2)y x =的曲线，并使用图形窗口的图形面板、绘图浏览器和属性编辑器添加文字和箭头。

解：编程如下： x=0:0.01*pi:2*pi; y=sin(2*x); plot(x,y)

20.绘制cos(10)y x =的曲线，保存图形文件为.bmp 和.jpg 格式，并在其他图像软件中查看该图形。

解：编程如下： x=0:0.01*pi:2*pi; y=cos(10*x); plot(x,y)

21.绘制z=维网线图和表面图，x在[-5,5]范围，y在[-5,5]范围，将网线图用gray色图并用颜色条显示色图，改变视角显示二维图形。

解：编程如下：

x=-5:5;

y=x;

[X,Y]=meshgrid(x,y);

Z=sqrt(X.^2+Y.^2);

figure(1);

mesh(X,Y,Z)

colormap gray

colorbar

7

6

5

4

3

2

1 figure(2);

surf(X,Y,Z)

view(2)

matlab期末大作业

电气学科大类 Modern Control Systems Analysis and Design Using Matlab and Simulink Title: Automobile Velocity Control Name: 巫宇智 Student ID: U200811997 Class：电气0811

电气0811 巫宇智 Catalogue Preface (3) The Design Introduction (4) Relative Knowledge (5) Design and Analyze (6) Compare and Conclusion (19) After design (20) Appendix (22) Reference (22)

Automobile Velocity Control 1.Preface: With the high pace of human civilization development, the car has been a common tools for people. However, some problems also arise in such tendency. Among many problems, the velocity control seems to a significant challenge. In a automated highway system, using the velocity control system to maintain the speed of the car can effectively reduce the potential danger of driving a car and also will bring much convenience to drivers. This article aims at the discussion about velocity control system and the compensator to ameliorate the preference of the plant, thus meets the complicated demands from people. The discussion is based on the simulation of MATLAB. Key word: PI controller, root locus

北邮模式识别课堂作业答案(参考)

第一次课堂作业 1.人在识别事物时是否可以避免错识 2.如果错识不可避免，那么你是否怀疑你所看到的、听到的、嗅 到的到底是真是的，还是虚假的 3.如果不是，那么你依靠的是什么呢用学术语言该如何表示。 4.我们是以统计学为基础分析模式识别问题，采用的是错误概率 评价分类器性能。如果不采用统计学，你是否能想到还有什么合理地分类 器性能评价指标来替代错误率 1.知觉的特性为选择性、整体性、理解性、恒常性。错觉是错误的知觉，是在特定条件下产生的对客观事物歪曲的知觉。认知是一个过程，需要大脑的参与.人的认知并不神秘，也符合一定的规律，也会产生错误 2.不是 3.辨别事物的最基本方法是计算 . 从不同事物所具有的不同属性为出发点认识事物. 一种是对事物的属性进行度量，属于定量的表示方法(向量表示法 )。另一种则是对事务所包含的成分进行分析，称为定性的描述(结构性描述方法)。 4.风险 第二次课堂作业 作为学生，你需要判断今天的课是否点名。结合该问题(或者其它你熟悉的识别问题，如”天气预报”)，说明: 先验概率、后验概率和类条件概率 按照最小错误率如何决策 按照最小风险如何决策 ωi为老师点名的事件,x为判断老师点名的概率 1.先验概率: 指根据以往经验和分析得到的该老师点名的概率,即为先验概率 P(ωi ) 后验概率: 在收到某个消息之后，接收端所了解到的该消息发送的概率称为后验概率。 在上过课之后,了解到的老师点名的概率为后验概率P(ωi|x) 类条件概率:在老师点名这个事件发生的条件下,学生判断老师点名的概率p(x| ωi ) 2. 如果P(ω1|X)>P(ω2|X)，则X归为ω1类别 如果P(ω1|X)≤P(ω2|X)，则X归为ω2类别 3.1)计算出后验概率 已知P(ωi)和P(X|ωi)，i=1,…，c，获得观测到的特征向量X 根据贝叶斯公式计算 j=1,…，x

matlab第三章习题解答

3-1 >> N=[1 2 3 4 5] N = 1 2 3 4 5 (1) >> N1=2*N N1 = 2 4 6 8 10 (2) >> N2=N/2 N2 = 0.5000 1.0000 1.5000 2.0000 2.5000 (3) >> N3=1./N N3 = 1.0000 0.5000 0.3333 0.2500 0.2000 (4) >> N4=1./(N.*N) N4 = 1.0000 0.2500 0.1111 0.0625 0.0400 3-2 m=input('请输入一个三位整数'); i=fix(m/100); j=mod(fix(m/10),10); k=mod(m,10); n=100*k+10*j+i; disp('反向输出的数为'); disp(n); 3-3 用if来实现 score=input('请输入成绩'); if score>100|score<0 disp('成绩必须>=0且<=100'); elseif score>=90&score<=100 disp('A');

elseif score>=80&score<=89 disp('B'); elseif score>=70&score<=79 disp('C'); elseif score>=60&score<=69 disp('D'); else disp('E'); end 用switch来实现 score=input('请输入成绩'); switch score case num2cell(90:100) disp('A'); case num2cell(80:89) disp('B'); case num2cell(70:79) disp('C'); case num2cell(60:69) disp('D'); case num2cell(0:59)%case{0,1,2,3,4,5}也可以 disp('E'); otherwise disp('成绩必须>=0且<=100'); end 注意:以下程序是错的 score=input('请输入成绩'); switch fix(score/10) case {9,10} disp('A'); case {8} disp('B'); case {7} disp('C'); case {6} disp('D'); case num2cell(0:5)%case{0,1,2,3,4,5}也可以

西南大学含弘在线课堂PowerPoint多媒体课程第一次作业答案

一、 单项选择题(本大题共5小题，共25.0分) 1. 传播信息的电缆、电磁波等称为： (5.0分 ) A ：媒介 选择： B B ：媒材 C ：媒体 D ：以上三个均不对 2. 在演示文稿制作时，如果在每一页幻灯片上添加版权信息，则可以在下列哪一种视图模式中进行操作设置： (5.0分 ) A ：讲义母版 选择： B B ：幻灯片母版 C ：备注母版 D ： 备注页视图 3. 在演示文稿中，如果仅需要查看幻灯片内的文字信息，而不需要观看其缩略图内容，则可以选择下列哪一个选项卡： (5.0分 ) A ：大纲选项卡 选择： A B ：幻灯片选项卡

C ：属性选项卡 D ：以上三个均不对 4. 在计算机中用构成一幅图的位置、维数、大小和形状的指令来进行描述的是： (5.0分 ) A ：图形 选择： B B ：图像 C ：图形与图像 D ：以上三个均不对 5. 适用于表现比较细节、层次以及丰富色彩的是 ： (5.0分 ) A ：图形 选择： C B ：图像 C ：图形与图像 D ：以上三个均不对 二、 多项选择题(本大题共2小题，共10.0分) 1. 多媒体技术所涉及到以下哪些形式: (5.0分 ) 选择：A,B,C

A ：媒介 B ：媒材 C ：媒体 D ：以上三个均不对 2. 评价多媒体课件质量的指标有 : (5.0分 ) A ：科学性 选择： A,B,C,D B ：技术性 C ：艺术性 D ：教育性 三、 判断题(本大题共3小题，共15.0分) 1. 多媒体课件在开发时其艺术性必须服务于教育性。 (5.0分) √ × 选择：√ 2. 多媒体课件在开发时要求必须使用多种媒体元素。 (5.0分) √ × 选择：× 3. 多媒体课件脚本要严格按照模板制作，不能有个性化调整。 (5.0分) √ × 选择：× 四、

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仿真期末大作业

MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：

2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。

仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。

第一次作业

数理经济学第一次作业题 董志勇 1. You can write your answer in either English or in Chinese or in both. 2. You have to do this paper independently; giving identical answers will be assumed cheating, which will lead to severe punishment. 3. Answer to the point; be precise and avoid unnecessary length. 1、单个垄断企业分别在两个市场上实现价格歧视，求证其曲线从下面与各曲线的水平之和相交。 2、设，为两种投入要素，共同生产一种产品，该产品的价格是，生产函数是，的价格是，的价格是，，是使得利润最大化的要素投入量，求证：和不会同时小于0。 3、生产函数是，求出要素需求曲线，对根据该生产函数进行生产而且追求利润最大化的企业进行比较静态分析，证明对该企业而言，交叉效应项的符号为负。 4、一个垄断者在两个市场上销售同一种产品，市场1的收益为，市场2的收益为，两个市场的总成本为，假设在市场1中，每销售一个单位的产品要征税t。 (a)证明t增加将会使得市场1销售得产品量降低。 (b)有关t增加对市场2中的产品的影响，最大化假设意味着什么？ (c)即使假设二阶条件总是成立，市场1中的税收的增加也可能导致总产出，请说明在什么情况下会是这样的。 (d)假设根据政府法律，市场2的产出固定在前面所述的利润最大化水平上，与不加管制的情况相比，在受到政府管制的情况下，市场1中的产出对税收增加得反应在绝对量上较小，请加以证明，并对此进行直观的经济解释。 5、垄断者在两个市场中销售其产品，收入函数分别为和。总产出的函数。在两个市场中，每销售一个单位的产品，征收同样的税收t。 (a)求解，和。如果有解的话，利润最大化隐含了哪些偏导的符号？ (b)假设产出固定不变，求解，能够确定该式的符号？

Matlab 第三章 习题 答案

第三章 3.5计算多项式乘法()()222254x x x x ++++ >> a=[1 2 2]; >> b=[1 5 4]; >> c=conv(a,b) c = 1 7 16 18 8 poly2sym(c) ans = x^4 + 7*x^3 + 16*x^2 + 18*x + 8 3.6计算多项式的除法()()3231368/4x x x x ++++ >> a=[3 13 6 8]; >> b=[1 4]; >> [q r]=deconv(a,b) q = 3 1 2 r = 0 0 0 0 >> poly2sym(q) ans = 3*x^2 + x + 2 3.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.0000 0 >> poly2sym(ans) ans = (4*x^5)/5 - 3*x^4 - (14*x^3)/3 + (5*x^2)/2 + 9*x

3.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) - 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 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求线性方程组223430 x xy y x x ?++=??-+=??的解 eqn1='x^2+x*y+y=3'; >> eqn2='x^2-4*x+3=0'; [x,y]=solve(eqn1,eqn2) x = 1 3 y = 1 -3/2 3.11求微分方程64dy y x dx =+的通解。 >> y=dsolve('Dy=6*y+4*x','x') y = (C3*exp(6*x))/9 - (2*x)/3 - 1/9 3.14求微分方程组2543dy x dt dx y dt ?=+????=-??在y(0)=3,x(0)=1时的特解。 >> eqn1='Dy=2*x+5'; >> eqn2='Dx=4*y-3'; >> [x y]=dsolve(eqn1,eqn2,'y(0)=3,x(0)=1') x = (2^(1/2)*(7*2^(1/2) - 9))/(8*exp(2*2^(1/2)*t)) + (2^(1/2)*exp(2*2^(1/2)*t)*(7*2^(1/2) +

matlab课后习题答案第四章

第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分 diff或数值梯度gradient指令计算)(t y'曲线 y'，然后把)(t y和)(t 绘制在同一图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1);

yc=gradient(y)/dt; plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =

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) 请输入一个三位数:489 984 第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<=100 disp('A'); elseif s>=80&s<=89 disp('B'); elseif s>=70&s<=79 disp('C'); elseif s>=60&s<=69 disp('D'); elseif s>=0&s<60 disp('E'); else disp('false'); end 请输入成绩分数:89 B Switch语句: >> 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 请输入成绩分数:88 B 第3题 输入20个数,求其中最大数与最小数。要求分别用循环结构与调用MA TLAB的max函数、min函数来实现。 1循序结构 A=[1 2 3 4 5 6 7 23 32 12 53 1 23 53 5 235 31 8888 33 3 21]; %还可用A=input('ENTER A:'); min=A(1); max=A(1); for i=1:19 if A(i+1)>max max=A(i+1); end if A(i+1)f; f=m; end end end funtion f=min(A) f=A(1); for i=1:20; %还可用A=input('ENTER A:'); m=A(i); if m

人教版小学四年级下册语文课堂作业本答案

语文园地二 1、精湛财富属于朴实器重喉结例外赞赏闷闷不乐狼吞虎咽面呈菜色 2、维持生计得到宽容回忆往事动作熟练神情严肃生活拮据 3、(1)行果 (2)言而有信 (3)勿施于人学而不厌，诲人不倦 4、(1)深信(2)自信(3)守信(4)诚信 (1)寡言少语(2)豪言壮语(3)花言巧语(4)胡言乱语 5、(1)严肃道德难题，把车留下来还是还给库伯 （2）小明这次语文考得很差，试卷发下来后，他战战兢兢地递给了妈妈。妈妈的脸顿时由晴转阴，大声嚷道：“你上课都听了些什么呀？这么简单的题目也会错？” 6、(1)不小心，文中指小男孩不小心掉进了大西洋。肯定，一定文中指船长肯定孩子掉进了大海。(2)肯定/确定用力/尽力/奋力的确/确实 (3)他实在游不动了，准备放弃，加上没人发现他掉进大海里，他觉得自己活不下去了认为船长一定会返航去救他的小男孩相信船长一定会去救他，说明船长被孩子所信任 (4)一个人能被他人相信也是一种幸福。他人在绝望时想起你，相信你会给予拯救更是一种无可替代的幸福，我们要真诚地去对待别人。 九、自然之道 1、沙滩侦察蠢事棒球帽愚不可及无遮无拦气喘吁吁 2、啄逐逐啄彻切彻切 4、(1)这个词突出了向导非常不情愿，生怕惹出麻烦，对同伴们善良的呼喊非常理解，但也感到非常无奈 (2)这个词突出了向导悲伤地叹息，对自

己之前的所作所为感到非常羞愧，表达了“自己”自责、悔恨的心理。5、(1)看着说飞走答道抱起 (2)R 艹 5 ① T 页 7 ③ (3)眼前的一幕是指，当幼龟踌躇不前时，一只嘲鸫突然飞来，它用尖嘴啄幼龟的头，企图把它拉到沙滩上去。 (4) ③ 十一、蝙蝠和雷达 1、清朗蝙蝠飞蛾避开蚊子揭开障碍荧光屏 2、bǔ líng qī zhuàng cāng ruì 3、启发灵敏(灵活) 仿照清晰 4、横七竖八地拉了许多绳子系着许多铃铛 5.(1)√×× (2)①波浪②天上的云/春天的风小羊/妈妈的手飘来飘去/吹在脸上真舒服。 （3）蝙蝠用嘴发射超声波，超声波向四周传播，遇到障碍物就反射回来，反射回来的那部分被蝙蝠用耳朵接收，经过大脑处理，就知道哪里有障碍物了。 6、嘴巴超声波耳朵青蛙－电子蛙眼蜻蜓－飞机蝴蝶－卫星螳螂－镰刀顺风耳－电话

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) 请输入一个三位数:489 984 第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<=100 disp('A'); elseif s>=80&s<=89 disp('B'); elseif s>=70&s<=79 disp('C'); elseif s>=60&s<=69 disp('D'); elseif s>=0&s<60 disp('E'); else disp('false'); end 请输入成绩分数:89 B Switch语句： >> 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 请输入成绩分数:88 B 第3题 输入20个数，求其中最大数和最小数。要求分别用循环结构和调用MATLAB的max函数、min函数来实现。 1循序结构 A=[1 2 3 4 5 6 7 23 32 12 53 1 23 53 5 235 31 8888 33 3 21]; %还可用A=input('ENTER A:'); min=A(1); max=A(1); for i=1:19 if A(i+1)>max max=A(i+1); end if A(i+1)f; f=m; end end end funtion f=min(A) f=A(1); for i=1:20; %还可用A=input('ENTER A:'); m=A(i); if m

《科学计算与MATLAB》期末大作业

杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业

给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t)=sin(-pi/4)));%函数表示 plot(t,y); %画图 end

2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]]

A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵

人教版小学四年级下册语文课堂作业本答案

人教版小学四年级下册语文课堂作业本答案 一、古诗词三首 1、敬亭山洞庭湖青螺 2.×√××√ 3.(1)完了，没有了悠闲孤独众孤因怀才不遇而感到的孤独、寂寞的心情 (2)月光下的洞庭湖洞庭湖的君山翠白青 (3)胜过仿佛，像蓝草，一种植物，叶蓝绿色，可提取青蓝色染料5、江南春唐杜牧千里莺啼绿映红，水村山郭酒旗风。南朝四百八十寺，多少楼台烟雨中。 二、桂林山水 1、无瑕攀登屏障骆驼危峰兀立波澜壮阔奇峰罗列 2、sì juàn dào m?ng 3、波澜壮阔水平如镜峰峦雄伟红叶似火衣服（图画）山峰（山峦） 4、(1)秋天，落叶仿佛像一个会跳舞的精灵，从树上翩翩飞舞而下。 (2)老房子的围墙真危险呐，手一碰，墙壁上的石灰好像要掉下来似的。 5、(1)M 辶 6 ④ D 斤 7 ② (2)静清绿奇秀险 (3) 绿树红花竹筏小舟云雾山水 6、田野上的野花真多啊，多得像天上的星星数也数不清；田野上的野花真美啊，美得仿佛那是一条五彩的锦缎；田野上的野花真香啊，香得引来群群蜜蜂争相采蜜！（我游览过水平如镜的西湖，玩赏过飞泻而下的瀑布，却从没看到过秦淮河的水。秦淮河的水真凉啊，凉得让你想马上和河水来

个亲密的接触；秦淮河的水真清啊，清得就像是一块彩色的大玻璃；秦淮河的水真亮啊，当河边的万家灯火映在水中，像是水底铺满了熠熠生辉的大钻石，使河水变得无比明亮、美丽；秦淮河的水真野啊，刚才还是块大玻璃，现在就波涛起伏，使得河上的画舫也一个劲儿地摇摆，再加上朱雀桥洞内随风而舞的狂草，更体现了秦淮河的粗犷之美；秦淮河的水真神啊，要不，它怎么会吸引来那么多的小雨点儿？想当年大自然在创造世界时，一定挥舞着最大的画笔，而放弃了那些精细的刻刀，毫无雕琢地在南京勾勒出了这样一条粗犷、神奇、美丽的秦淮河吧？） 三、记金华的双龙洞 1、浙江臀部稍微额角擦伤蜿蜒依据石钟乳 2、一簇臀部稍微蜿蜒民族宫殿树梢大碗 3、(1)②① (2)①② 4、（1）小妹妹总是时而笑，时而哭，时而忧，时而喜，真是捉摸不透。（六 月的天气，时而风，时而雨，时而凉，时而热，真是变化多端。）（2）即使这件事情再困难，我们也要尽自己最大的力量去做好它。 5、(1)根据名称欣赏 (2)，，、、， (3)在洞顶的那两条龙，一条黄龙，一条青龙。 (4)蜿蜒在洞顶的双龙变化多端，颜色各异的石钟乳和石笋。 6、怎样美的笔盒呢？让人一看就爱不释手，是这样美的笔盒。（怎样简陋的小屋呢？裂缝的石灰墙，风化的水泥地，一扇破旧的窗，一张油漆脱落的书桌，根本找不到一样值钱的东西，就是这样一间简陋的小屋。） 语文园地一 1、洞庭玩赏扩散泰山拥挤孔隙仰卧油桐水平如镜拔地而起突兀

MATLAB大作业

选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c

(2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵：

（1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2

第4章matlab程序流程控制_习题答案

第4章 MATLAB程序流程控制 习题4 一、选择题 1．下列关于脚本文件和函数文件的描述中不正确的是（）。 A A．函数文件可以在命令行窗口直接运行 B．去掉函数文件第一行的定义行可转变成脚本文件 C．脚本文件可以调用函数文件 D．函数文件中的第一行必须以function开始 2．下列程序的输出结果是（）。D y=10; if y==10 y=20; elseif y>0 y=30 end disp(y) A．1 B．30 C．10 D．20 3．有以下语句： a=eye(5); for n=a(2:end,:) for循环的循环次数是（）。C A．3 B．4 C．5 D．10 4．设有程序段 k=10; while k k=k-1 end 则下面描述中正确的是（）。A A．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段： x=reshape(1:12,3,4); m=0;

n=0; for k=1:4 if x(:,k)<=6 m=m+1; else n=n+1; end end 则m和n的值分别是（）。C A．6 6 B．2 1 C．2 2 D．1 2 6．调用函数时，如果函数文件名与函数名不一致，则使用（）。A A．函数文件名 B．函数名 C．函数文件名或函数名均可 D．@函数名 7．如果有函数声明行为“function [x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是（）。B A．x=f1(a,b,c) B．[x,y,z,w]=f1(a,b,c) C．[x,b,z]=f1(a,y,c) D．[a,b]=f1(x,y,z) 8．执行语句“fn=@(x) 10*x;”，则 fn是（）。A A．匿名函数 B．函数句柄 C．字符串 D．普通函数 9．执行下列语句后，变量A的值是（）。D >> f=@(x,y) log(exp(x+y)); >> A=f(22,3); A．22,3 B．22 C．3 D．25 10．程序调试时用于设置断点的函数是（）。A A．dbstop B．dbclear C．dbcont D．dbstack 二、填空题 1．将有关MATLAB命令编成程序存储在一个扩展名为.m的文件中，该文件称为。M文件 2．有语句“for k=[12;34]”引导的循环结构，其循环体执行的次数为。1 3．MATLAB中用于控制不确定重复次数的循环语句为，若在循环执行过程中需要终止该循环时采用的语句为。while…end，break 4．函数文件由语句引导。在函数定义时，函数的输入输出参数称为参数，简称。在调用函数时，输入输出参数称为参数，简称。 function，形式，形参，实际，实参 5．在MATLAB中，函数文件中的变量是变量。定义变量是函数间传递信息的一种手段，可以用命令定义。局部，全局，global 6．应用程序的错误有两类，一类是错误，另一类是运行时的错误，即错

(完整版)有限元大作业matlab---课程设计例子

有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师：

连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。

[程序原理及实现]： 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下： 问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法

输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

课堂作业本答案

1、看拼音，写词语。 书籍抽屉窝头粗暴魔鬼苦刑拘留解救精兵简政死得其所重于泰山轻于鸿毛 2、多音字组词。 （转移、转身）转（转动、旋转） （兴旺、兴盛、兴起）兴（高兴、兴趣、即兴） 3、体会下列词语的构词特点，并照样子写一写。 乱蓬蓬热腾腾（冷清清）（光秃秃）（静悄悄）（白茫茫） 跳跃陈设（（奔跑）（肥胖）（吵闹）（寻找）两个字意思相近） 4、从“词语盘点”中选择合适的词语填空。 （残暴）的匪徒（尖锐）的问题（含糊）的回答（恐怖）的氛围 （陷入）困境（解救）人质（审阅）文件（拘留）嫌犯 5、阅读回顾。 （1）毛泽东为人民服务《十六年前的回忆》李大钊同志无私无畏、坚贞不屈的革命精神（或灯光献身的情景等）

（2）《》首尾呼应突出了作者一直把父亲的被害牢郝副营长为了后代的幸福和安乐英勇记在心的思想感情，加深了我们的印象。 6、学习“日积月累”，完成练习。 （1）叶挺被国民党关押（2）① 7、默读短文，完成练习。 （1）正确的读音打勾：sìhècánɡ （2）解释词语：矜持不苟——端庄严谨，毫不马虎,态度认真，一丝不苟。齿冷——耻笑，讥笑。指不光彩，不正当的行为让人鄙视耻笑。 （3）摘录文中的四字词语，并选一个写句子。 矜持不苟舍己为公意外之财毫无所得齿冷三天洁白朴素等 例1：舍己为公——周总理工作劳苦，生活俭朴，舍己为公的精神令我们敬佩不已。 例2：毫无所得——三天打鱼，两天晒网，这样的学习将毫无所得的。 （4）描写国民党士兵言行的词句：①哪知道从我上身摸到下身……一个铜板也没搜出。②于是他们激怒起来……不然炸死你！”③“你骗谁……会没有钱！” ④他们又在周围低头搜寻了一番……就带着我走了。他们的言行反衬出方志敏清贫朴素、矜持不苟、舍己为公的精神