文档库 最新最全的文档下载
当前位置:文档库 › 西南交大MATLAB编程作业相关复习资料

西南交大MATLAB编程作业相关复习资料

西南交大MATLAB编程作业相关复习资料
西南交大MATLAB编程作业相关复习资料

编程题目

《历年竞赛题目及编程培训》P(135至137)

1.编写一个程序,要求输入五个整数,然后由小到大排序再输出.

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t1中

%然后在命令窗口定义 x 为所输入的5个整数的行向量

形式,再输入 t1(x)

function f=t1(x)

for i=1:4

for j=(i+1):5

if x(i)>x(j)

a=x(i);

x(i)=x(j);

x(j)=a;

end

end

end

disp(x)

2.将一个整型数组的元素按逆序重新存放(如原序为:8,6,5,4改为:4,5,6,8)。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t2中

%然后在命令窗口定义 x 为所输入的数组,再输入

t2(x)

function t2(x)s=length(x);

a=fix(s/2);

for i=1:a

j=s-i+1;

c=x(i);

x(i)=x(j);

x(j)=c;

end

disp(x)

3.输入一个字符,如果是大写字母,则将其转换成小

写并输出;若是小写则直接输出;若是非字母字符则

打印:“Dataerror!”。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t3中

%然后在命令窗口直接输入 t3

function t3()

a=input('输入一个字符:string');%输入的形式为

'x' x为输入的字符

if a<=57&a>=48

fprintf('“Dataerror!”\n')

elseif a<=90&a>=65

fprintf('%c\n',a+32)

elseif a<=122&a>=97

a

end

4.输入一个整数,写一程序输出它是几位数。

%此程序可以计算一万位整数的位数

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t4中

%然后在命令窗口直接输入 t4(x) x为所输入的整数

function t4(x)

i=0;

for j=1:10000

x=abs(x/10);

a=fix(x);

if a>0

i=i+1;

end

if a==0

i=i+1;

break

end

end

disp(i)

5.写一程序求1!+2!+ (10)

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t5中

%然后在命令窗口直接输入t5

function t5()

j=1;s=0;

for i=1:10

j=i*j;

s=s+j;

end

disp(s)

6.从键盘上输入a与n的值,计算

sum=a+aa+aaa+aaaa+…(共n项)的和。

例a=2,n=4,则sum=2+22+222+2222.

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t6中

%然后在命令窗口直接输入 t6(a,n) a和n分别是一个数function t6(a,n)

x=a;t=a;

for i=1:(n-1)

t=t+10^i*a;

sum=x+t;

end

sum

7.编程求数列1,1/2 ,1/3 ,1/4 ,1/5 ,……的所有大于等于0.000001的数据项之和并输出结果

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t7中

%然后在命令窗口直接输入 t7

function t7()

s=0;

for c=1:100000

d=1/c;

s=s+d;

end

s

8.求3*3矩阵的主对角线元素之和。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t8中

%然后在命令窗口直定义 x 为一个三行三列矩阵%最后输入 t8(x)

function t8(x)

s=0;

for i=1:3;

j=i;

s=s+x(i,j);

end

s

9.从键盘输入的10个整数中,找出第一个能被7

整除的数。若找到,打印此数后退出循环;若

未找到,打印“not exist”。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t9中

%然后在命令窗口定义 x为输入的10个数的行向量形式

%最后输入 t9(x)

function t9(x)

for i=1:10

a=x(i)/7;

b=fix(a);

if a-b==0

disp(x(i))

break

end

if i==10

fprintf('not exist\n');

end

end

10. 设计一个用于计算个人所有税的程序。假设

个人所得税的缴纳标准为:月收入少于等于800

元者不纳税;超出800元的部分,纳税5%;超出

2000元的部分,纳税10%;超出5000元的部分,

纳税20%;超出10000元的部分,纳税30%;超出

100000元的部分,纳税40%。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t10中

%然后在命令窗口直接输入 t10(x) x为月收入

function g=t10(x)

b=90000*0.3;

c=5000*0.2;

d=3000*0.1;

e=1200*0.05;

if x>=100000

a=(x-100000)*0.4;

g=a+b+c+d+e;

elseif x>=10000

a=(x-10000)*0.3;

g=a+c+d+e;

elseif x>=5000

a=(x-5000)*0.2;

g=a+d+e;

elseif x>=2000

a=(x-2000)*0.1;

g=a+e;

elseif x>=800

g=(x-800)*0.05;

else

g=0;

end

elseif x>=5000

a=(x-5000)*0.2;

g=a+d+e;

elseif x>=2000

a=(x-2000)*0.1;

g=a+e;

elseif x>=800

g=(x-800)*0.05;

else

g=0;

end

11.编写一个程序,用循环结构输出以下图形。

&&&&&&&&&

&&&&&&&

&&&&&

&&&

&

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t11中

%然后在命令窗口直接输入 t11

function t11()

for i=9:-2:1

for j=1:i

fprintf('&')

end

fprintf('\n');

end

12.编写一个程序求出1000以内的所有素数。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t12中%然后在命令窗口直接输入 t12

function t12()

x=[];

for i=2:1000

t=1;j=0;

for j=2:i-1

a=i/j-fix(i/j);

if a==0

t=0;

break

end

end

if t>0

x=[x,i];

end

end

disp(x)

13. 输入两个整型变量a,b的值,输出

a+b,a-b,a*b,a/b,的结果,要求连同算式一起输出,

每个算式占一行。

例:a等于10,b等于5,a+b的结果输成

10+5=15

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t13中

%然后在命令窗口直接输入 t13(a,b) a,b为要输入的两

个整形变量

function t13(a,b)

c=a+b;

d=a-b;

e=a*b;

f=a/b;

fprintf('%d+%d=%d\n%d-%d=%d\n%d*%d=%d\n%d/%d=%d\

n',a,b,c,a,b,d,a,b,e,a,b,f);

14. 求前驱字符和后继字符。输入一个字符,

找出它的前驱字符和后继字符,并按ASCII码

值,按从大到小的顺序输出这三个字符及其对

应的ASCII码值。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t14中

%然后在命令窗口直接输入 t14('x') x为输入的字母

function t14(x)

if x=='a'

x=='z';

else a=x-1;

end

if x=='z'

x='a';

else b=x+1;

end

fprintf('%c %d,%c %d,%c %d \n',b,b,x,x,a,a);

15. 输入三角形的边长,写一求三角形面积的函数。

(面积=sqrt(s(s-a)(s-b)(s-c)), s=(a+b+c)/2)

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t15中

%然后在命令窗口直接输入 t15(a,b,c) a,b,c分别是三

角形的边长

function t15(a,b,c)

s=(a+b+c)/2;

d=sqrt(s*(s-a)*(s-b)*(s-c))

fprintf('面积=%d\n',d)

16. 输入某个点A的平面坐标(x,y),判断(输出)A 点是在圆内、圆外还是在圆周上,其中圆心坐标为(2,2),半径为1 。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t16中

%然后在命令窗口直接输入 t16(x,y) x,y为点A的横纵坐标

function t16(x,y)

a=(x-2)^2;

b=(y-2)^2;

c=sqrt(a+b);

if c==1

fprintf('A点在圆上\n')

end

if c>1

fprintf('A点在圆外\n')

end

if c<1

fprintf('A点在圆内\n')

end

17. 求爱因斯坦数学题。有一条长阶,若每步跨2阶,则最后剩余1阶;若每步跨3阶,则最后剩2阶;若每步跨5阶,则最后剩4阶;若每步跨6阶,则最后剩5阶;若每步跨7步,最后正好一阶不剩。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t17中

%然后在命令窗口直接输入 t17(x) x为台阶数

function t17(x)

s=[];

for i=7:x

if mod(i,2)==1

if mod(i,3)==2

if mod(i,5)==4

if mod(i,6)==5

if mod(i,7)==0

s=[s,i];

end

end

end

end

end

end

disp(s)

18.输入一串字符,直到输入一个星号(*)为

止,统计(输出)其中的字母个数和数字字符

个数。(要求可输入任意字符)

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t18中

%然后在命令窗口直接输入 t18('x') x为输入的字符串,

必须以*号结束

function t18(x)

a=numel(x);

j=0;k=0;

for i=1:a

if x(i)=='*';

break

elseif x(i)>= 'a'&x(i)<='z'| x(i)>=

'A'&x(i)<='Z'

j=j+1;

elseif x(i)>= '0'&x(i)<='9'

k=k+1;

end

end

fprintf('字母个数为%d\n数字字符个数

为%d\n',j,k);

19. 从键盘输入一个正整数n,编程判断这个数

是否同时含有奇数字和偶数字。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t19中

%然后在命令窗口直接输入t19(x)其中x为正整数

function t19(x)

a=num2str(x);

b=numel(a);

g=0;k=0;

for i=1:b

c=mod(a(i),2);

if c==0

k=1;

break

end

end

for i=1:b

c=mod(a(i),2);

if c>0

g=1;

break

end

end

if g==1&k==1

fprintf('Yes!\n');

end

if g==0|k==0

fprintf('No!\n');

end

20. 设某县2000年工业总产值为200亿元,如果该县预计平均年工业总产值增长率为4.5%,那么多少年后该县年工业总产值将超500亿元?

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t20中

%然后在命令窗口直接输入 t20

function t20()

sum=200;a=0;

for i=1:1000

a=sum*0.045;

sum=sum+a;

if sum-200>=300

break

end

end

fprintf('%d\n',i);

21. (1)输入一个三位数,判断是否是一个“水仙花数”。水仙花数是指三位数的各位数字的立方和等于这个三位数本身。例如: 153=1*1*1+5*5*5+3*3*3%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t211中

%然后在命令窗口直接输入 t211(x) x为一个三位数

function t211(x)

y=x;

a=mod(x,10);

x=x/10;

b=fix(mod(x,10));

x=x/10;

c=fix(x);

if y==a*a*a+b*b*b+c*c*c

fprintf(' %d是水仙花数\n',y)

else fprintf(' %d不是水仙花数\n',y)

end

(2)一个整数等于该数所有因子之和,则称该

数是一个完数。例如:6和28都是完数。因

为:6=1+2+3,28=1+2+4+7+14. 输出三位数中所

有完数。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t212中

%然后在命令窗口直接输入 t212

function t212()

for i=100:999

sum=0;

for j=1:i-1

if mod(i,j)==0

sum=sum+j;

end

end

if sum==i

disp(sum)

continue

end

end

22. 输入10个数,将10个整数按生序排列输出,并且

奇数在前,偶数在后。如果输入的10个数是:10 9 8 7

6 5 4 3 2 1 ,则输出:1 3 5

7 9 2 4 6

8 10。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t22中

%然后在命令窗口将需输入的10个数存入数组 x 中,

%最后直接输入 t22(x)

function t22(x)

y=[];d=[];k=1;a=0;

for i=1:10

if mod(x(i),2)~=0

d(k)=x(i);

k=k+1;

end

end

t(d)

k=1;

for i=1:10

if mod(x(i),2)==0

d(k)=x(i);

k=k+1;

end

end

t(d)

function t(d)

ii=numel(d);

for i=1:ii

y(i)=d(i);

end

for j=1:ii

for r=1:ii-1

if y(r)>y(r+1)

t=y(r);

y(r)=y(r+1);

y(r+1)=t;

end

end

end

fprintf('%d ',y)

23.将数组a[n]中的每一个元素依次循环向后移动一位。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t23中

%然后在命令窗口直接输入 t23

function t23

x=0;

x=input('输入一个数组a[],输入格式为[a1 a2

a3 ……]\n','s');

x=str2num(x);

if x~=0

t(x)

else

fprintf('Not input!');

end

function f=t(x)

a=[];k=0;b=numel(x);

for i=1:b-1

a(i+1)=x(i);

end

a(1)=x(b);

for j=1:b

end

disp(a)

24.将一个正整数n输出成“千分位”形式,即从个

位数起,每三位之间加一个逗号,例如,将7654321

输出成7,654,321.

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t24中

%然后在命令窗口直接输入 t24

function t24

x=input('输入一个整数\n','s');

a=numel(x);

x=str2num(x);

b=ceil(a/3);

for i=1:b

a(i)=mod(x,1000);

x=fix(x/1000);

end

for j=b:-1:2

fprintf('%d,',a(j));

end

fprintf('%d',a(1));

fprintf('\n');

25. 求出矩阵a中每行元素的平均值和平均最

大的行号。

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t25中

%然后在命令窗口定义矩阵a为需要输入的矩阵,再输入

t25(a)

function t25(a)

m=[];n=[];max=0;k=0;p=[];

c=numel(a);

m=a(1,:);

b=numel(m);

d=c/b;

for i=1:d

n(i)=sum(a(i,:));

end

for i=1:d

p(i)=n(i)/b;

if n(i)>max

max=n(i);

k=i;

end

end

fprintf('每行元素的平均值为:\n');

for i=1:d

fprintf('%d ',p(i));

end

fprintf('\n');

fprintf('平均值最大的行号为: %d\n',k);

26. 编写一个函数,找出矩阵A中最大元素和最小元

素所在的行列号.

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t26中

%然后在命令窗口将需要处理的矩阵存入 a 中

%最后输入 t26(a)

function t26(a)

e=1;f=1;

g=1;k=1;

c=numel(a);

m=a(1,:);

b=numel(m);

d=c/b;

min1=a(1,1);

max1=a(1,1);

for i=1:d

for j=1:b

if max1<=a(i,j)

max1=a(i,j);

e=i;

f=j;

end

if min1>=a(i,j)

min1=a(i,j);

g=i;

k=j;

end

end

end

fprintf('最大元素所在的行列号为%d %d\n',e,f); fprintf('最小元素所在的行列号为%d %d\n',g,k);

27.给出年、月、日,计算该日是该年的第几天?%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t27中

%然后在命令窗口将需要处理的年月日按年月日顺

序存入数组x 中

%最后输入 t27(x)

function t27(x)

a=0;

if

mod(x(1),4)==0&mod(x(1),100)~=0|mod(x(1),400)==0

a=1;

end

tian=30*(x(2)-1)+x(3);

for i=1:x(2)-1

if i==3|i==5|i==7|i==8|i==10

tian=tian+1;

end

end

if x(2)~=1|x(2)~=2

if a==0

tian=tian-1;

end

end

tian

28. 输入数组a,b的元素,用数组a和b构造数组

c,使得:

c(i)=a(i)-b(i) 当a(i)>b(i)时;

c(i)= a(i)+b(i) 当a(i)<=b(i)时;

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t28中

%然后在命令窗口直接输入 t28

function t28()

a=input('请输入数组a输入格式为[x1 x2 ……]\n');

b=input('请输入数组b,b与a的规格相同\n');

d=numel(a);c=[];

for i=1:d

if a(i)>b(i)

c(i)=a(i)-b(i);

elseif a(i)<=b(i)

c(i)=a(i)+b(i);

end

end

c

*30.编写程序求所有各位数字的立方和等于420的三位

数.

%首先在MATLAB编辑器中建立如下程序并保存在文件名

为t30中

%然后在命令窗口直接输入 t30

function t30

s=0;

for k=100:999

i=k;a=zeros(1,3);

for j=1:3

a(j)=mod(i,10);

i=fix(i/10);

end

if a(1)^3+a(2)^3+a(3)^3==420

disp(k)

s=1;

end

end

if s==0

fprintf('不存在这样的三位数\n');

end

*31.学校要举行篮球循环赛,共有十个队参加比赛,请你安排一下比赛场次,使每队每次打完比赛后至少能隔一场不比赛,以便得到休息从而保持体力,请用matlab解决这个问题,并给出至少一个比赛日程表(每队休息的场次越多越好,但每队得到的休息应大体相同,不能出现某队刚打完比赛,马上和另一个已经休息了诺干场的球队比赛,这是不公平的)。

%首先在MATLAB编辑器中建立如下程序并保存在文件名为t31中

%然后在命令窗口直接输入 t31

function t31

a=[1:5];

b=[10:-1:6];

for k=1:9

fprintf('第%d轮\n',k);

for i=1:5

fprintf('%2d-----%2d\n',a(i),b(i));

end

for j=5:-1:1

a(j+1)=a(j);

end

a(2)=b(1);

b(1)=[];

b(5)=a(6);

a(6)=[];

end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

知识点总结

了解matlab中常见变量:ans pi Nan eps等的基本含

Pi圆周率 eps浮点运算的相对精度10^-25 Nan不

定值 ans求极值

了解matlab中变量定义的规则。

1、将矩阵

?

?

?

?

?

?

=

5

7

2

4

a

?

?

?

?

?

?

=

3

8

1

7

b

?

?

?

?

?

?

=

2

6

9

5

c

合两个新矩阵:

求a和b的和,b和c的矩阵乘,a和c的数组乘。

组合成一个3

4?的矩阵,第一列为按列顺序排列的a

矩阵元素,第二列为按列顺序排列的b矩阵元素,第三

列为按列顺序排列的c矩阵元素,即

?

?

?

?

?

?

?

?

?

?

?

?

2

3

5

9

1

2

6

8

7

5

7

4

按照a、b、c的列顺序组合成一个行矢量,即

[]2

9

6

5

3

1

8

7

5

2

7

4

a=[4 2;7 5];b=[7 1;8 3];c=[5 9;6 2];...

a+b

ans=11 3

15 8

>>b*c

ans=41 65

58 78

>>a.*c

ans=20 18

42 10

d=[a(:),b(:),c(:)]

d=4 7 5

7 8 6

2 1 9

5 3 2

d(:)' %转置符号“‘”

ans =

4 7 2

5 7 8 1 3

5 6 9 2

法2

e1=cat(1,a(:),b(:),c(:))

>> e2=e1'

A=[4 9 2;7 6 4;3 5 7];

B=[1 2 3;4 5 6;7 8 9];

X=A\B

X=0.1043 0.1374 0.1706

-0.1754 -0.0948 -0.0142

1.0806 1.1517 1.2227

>>det(X) %计算矩阵行列式的值

ans=-3.1825e-017

>>eig(X) %计算特征值

ans=1.3385

0.0000

-0.1062

x=[0:0.02*pi:4*pi];%必须有时间域申明

y1=sin(x);y2=cos(x);...

plot(x,y1.*y2)

>>x1=input('x1='); y1=input('y1='); z1=input('z1=');

x2=input('x2=');...

y2=input('y2='); z2=input('z2=');

d=sqrt((x1-x2)^2+(y1-y2)^2+(z1-z2)^2)

x1=2

y1=3

z1=1

x2=8

y2=-5

z2=-1

d=10.1980

>>

t=[0:0.5:10];y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5* t);...

plot(t,y1,'-.b',t,y2,'-k',t,y3,':r')

>> grid on,title('指数曲线')

>> grid on,xlabel('t'),ylabel('y')

>> gtext('蓝色a=0.1 黑色a=0.2 红色a=0.5')

%以鼠标取点方式标注

3、求出任意矩阵a中每行元素的平均值和平均最大的行号。

m=[];n=[];max=0;k=0;p=[];a=input('请输入一个矩阵');

[d,b]=size(a);

for ii=1:d

n(ii)=sum(a(ii,:));

end

for ii=1:d

p(ii)=n(ii)/b;

if n(ii)>max

max=n(ii);

k=ii;

end

end

fprintf('每行元素的平均值为:\n');

for i=1:d

fprintf('%4.3f ',p(ii));

end

fprintf('\n');

fprintf('平均值最大的行号为: %d\n',k);

编写一个函数,找出任意矩阵A中最大元素和最小元

素所在的行列号.

e=1;f=1;

g=1;k=1;

a=input('请输入一个矩阵');

[d,b]=size(a);

min1=a(1,1);

max1=a(1,1);

for ii=1:d

for jj=1:b

if max1<=a(ii,jj)

max1=a(ii,jj);

e=ii;

f=jj;

end

if min1>=a(ii,jj)

min1=a(ii,jj);

g=ii;

k=jj;

end

end

end

fprintf('最大元素所在的行列号为%d %d\n',e,f);

fprintf('最小元素所在的行列号为%d %d\n',g,k);

飞船的运转半径。一个做匀速圆周运动的物体,其向

心加速度公式为

r

v

a

2

. a代表向心加速度,单位为

m/s*s。v 代表物体运动的速率,单位为m/s,,r 代表

半径,单位为m. 假设这个物体是一个飞机,回答下列

问题:

1、假设飞机的运动速度为0.85 马赫, 即声速的85%.

如果向心加速度2g, 那么飞机的半径为多少?

2、假设飞行员能忍耐的最大加速度为7g. 那么以1. 5

马赫飞行的最小半径为多少?

3、画出以向心加速度为自变量的半径函数, 向心加速

度的取值为[2g,6g], 假设运转速度为0.85 马赫.

程序:

1.

% Define variables:

% g -- Max acceleration (g)

% grav -- Acceleration of gravity

(9.81 m/s2)

% mach1 -- Speed of sound (340 m/s) % radius -- Turning radius (m)

% speed -- Aircraft speed in Mach % Initialise values

grav = 9.81;

mach1 = 340;

% Get speed and max g

speed = input('Enter speed (Mach): ');

g = input('Enter max acceleration (g): '); % Calculate radius

radius = (speed * mach1).^ 2 / ( g * grav ); % Tell user

fprintf('Turning radius = %f m\n',radius);

2.

% max_speed -- Maximum speed in Mach numbers

% min_speed -- Minimum speed in Mach numbers

% radius -- Turning radius (m)

% speed -- Aircraft speed in Mach

% Initialise values

grav = 9.81;

mach1 = 340;

% Get speed and max g

min_speed = input('Enter min speed (Mach): ');

max_speed = input('Enter min speed (Mach): ');

g = input('Enter max acceleration (g): ');% Calculate range of speeds

speed =

min_speed:(max_speed-min_speed)/20:max_s

peed;

% Calculate radius

radius = (speed * mach1).^ 2 / ( g * grav );

% Plot the turning radius versus speed

plot(speed,radius/1000);

title('Plot of turning radius versus

speed');

xlabel('Speed (Mach)');

ylabel('Turning radius (km)');

grid on;

3.

% max_g -- Maximum accleration in

g's

% min_g -- Minimum accleration in

g's

% radius -- Turning radius (m)

% speed -- Aircraft speed in Mach

% Initialise values

grav = 9.81;

mach1 = 340;

% Get speed and max g

speed = input('Enter speed (Mach): ');

min_g = input('Enter min acceleration (g):

');

max_g = input('Enter min acceleration (g):

');

% Calculate range of accelerations

g = min_g:(max_g-min_g)/20:max_g;

% Calculate radius

radius = (speed * mach1).^ 2 ./ ( g * grav );

% Plot the turning radius versus speed

plot(g,radius/1000);

title('Plot of turning radius versus

acceleration');

xlabel('Centripetal acceleration (g)');

ylabel('Turning radius (km)');

grid on;

Decibels. Engineers often measure the ratio of two

power measurements in decibels, or dB. The equation

for the ratio of two power measurements in decibels is

1

2

10

log

10

P

P

dB

Where P2 is the power level being measured, and P1 is

some reference power level.

A.Assume that the reference power level P1 is 1

milliwatt, and write a program that accepts an input

power P2 and converts it into dB with respect to 1 mW

reference level.(Engineers have a special unit for dB

power levels with respect to a 1 mW reference: dBm).

Use good programming practices in your program.

B.Write a program that creates a plot of power in

watts versus power in dBm with respect to a 1 mW

reference level. Create both a linear xy plot and a

log-linear xy plot.

程序1:

% Script file: decibel.m

%

% Define variables:

% dBm -- Power in dBm

% pin -- Power in watts

% Prompt the user for the input power.

pin = input('Enter the power in watts: ');

% Calculate dBm

dBm = 10 * log10( pin / 1.0e-3 );

% Tell user

disp (['Power = ' num2str(dBm) ' dBm']);

程序2:

%

% Define variables:

% dBm -- Power in dBm

% pin -- Power in watts

% Create array of power in watts

pin = 1:2:100;

% Calculate power in dBm

dBm = 10 * log10( pin / 1.0e-3 );

% Plot on linear scale

figure(1);

plot(dBm,pin);

title('Plot of power in watts vs power in dBm');

xlabel('Power (dBm)');ylabel('Power (watts)');

grid on;

% Plot on semilog scale

figure(2);

semilogy(dBm,pin);

title('Plot of power in watts vs power in

dBm');

xlabel('Power (dBm)');

ylabel('Power (watts)');

grid on;

Current Through a Diode. (Suggestion: you

had better use the

function file )

Figure. A semiconductor diode.

The current flowing through the semiconductor

diode(半导体二极管) shown in Figure is given

by the equation

(1)

D

qV

kT

D

i I e

=-

Where

D

V= the voltage across the diode, in volts

D

i= the current flow through the diode ,in

amps

I= the leakage current of the diode, in amps.

q = the charge on an electron, 1.602×

10-19coulombs(库仑)

k = Boltzmann's constant, 1.38×10-23joule/k

T = temperature, in kelvins (K)

The leakage current I0of the diode is

2.0uA (Hint:2×10-6A). Write a program to

calculate the current flowing through this diode

for all voltages from -1.0 V to +0.6V, in

0.1Vsteps. Repeat this process for the following

temperatures: 750F, 1000F and 1250F. Create a

plot of the current as a function of applied

voltage, with the curves for three different

temperatures appearing as different colors.

(Hint:The formula for the conversion

from temperature in degrees Fahrenheit to

absolute temperature in kelvins:

15

.

273

)0.

32

)

(

9

5

(

)

(+

-

=F

in

T

inKelvins

T )

% Define variables:

% i0 -- Leakage current (A)

% id -- Diode current (A)

% k -- Boltzmann's constant

(joule/K)

% q -- Charge on an electron

(coul)

% temp_f -- Temperature (deg F)

% temp_k -- Temperature (K)

% vd -- Diode voltage (V)

% Initial values

i0 = 2.0e-6; % amps

k = 1.38e-23; % joule/K

q = 1.602e-19; % Coulombs

vd = -1.0:0.01:0.6; % Volts

temp_f = [75 100 125]; % def F

for ii = 1:length(temp_f)

% Convert temperature to kelvins.

temp_k = (5/9) * (temp_f(ii) - 32) + 273.15;

% Calculate currents

id = i0 .* ( exp((q*vd)/(k*temp_k)) - 1 );

% Plot line in various colors

if ii == 1

plot(vd,id,'b-','LineWidth',2);

hold on;

elseif ii == 2

plot(vd,id,'k--','LineWidth',2);

elseif ii == 3

plot(vd,id,'r:','LineWidth',2);

hold off;

end

end

legend('75 F', '100F', '125 F')

grid on;title('\bfPlot of diode voltage vs diode

current');

xlabel('\bf\itv_{D}');

ylabel('\bf\iti_{D}');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5%%%

1、Fibonacci numbers. The n th Fibonacci

number is defined by the following recursive

equations:

f(1)=1

f(2)=2

f(n)=f(n-1)+f(n-2)

Therefore, f(3)=f(2)+f(1)=2+1=3, and so forth

for higher numbers. Write an M-file to

calculate and write out the n th Fibonacci

number for n>2, where n is input by the user.

Use a while loop to perform the calculation.

% Define variables:

% fn -- Fibonacci number

% n -- The item in the sequence to

calculate

% Get n

n = input('Enter the Fobonacci number n to

evaluate (n>2): ');

% Check to see that n is an integer greater

than two

if n <= 2

disp('Error--n must greater than two!');

elseif round(n) ~= n

disp('Error--n must be an integer!');

else

% Calculate fn

fn = zeros(1,n);

fn(1) = 1;

fn(2) = 2;

ii = 3;

while ii <= n;

fn(ii) = fn(ii-1) + fn(ii-2);

ii = ii + 1;

end

% Display result

disp(['The 'int2str(n) 'th Fibonacci

number = ' int2str(fn(n))]);

end

2、Geometric Mean. The geometric mean of a

set of numbers x1 through xn is defined as the

n th root of the product of the numbers Geometric mean=)..(321n n x x x x

Write a MATLAB program that will accept an arbitrary number of positive input values and calculate both the arithmetic mean(i.e, the average) and the geometric mean of the numbers. Use a while loop to get the input values and terminate the inputs when a user enters a negative number. Test your program by calculating the average and geometric mean of the four numbers 10, 5, 2, and 5.

% Define variables: % x -- Input values % g -- Geometric mean

% nvals -- Number of input values % prod -- Product of input values

% Initialize product and nvals prod = 1; nvals = 0;

% Read the first number

x = input('Enter first number: ');

% Read the remaining numbers while x > 0

prod = prod * x; nvals = nvals + 1;

x = input('Enter next number: '); end

% Calculate geometric mean g = prod ^ (1/nvals);

% Tell user

fprintf('The geometric mean is %.4f\n',g);

3、Harmonic Mean. The harmonic mean is yet another way of calculating a mean for a set of numbers. The harmonic mean of a set of numbers is given by the equation

ha

N

x x x N 1...1121+++

Write a MATLAB program that will read in an

arbitrary number of positive input values and calculate the harmonic mean of the numbers. Use any method that you desire to read in the input values. Test your program by calculating the harmonic mean of the four numbers 10, 5, 2 and 5.

% Define variables: % hmean -- Harmonic mean

% nvals -- Number of input values

% sumr -- Sum of reciprocals of the input values

% x -- Input values

% Initialize sum sumr = 0;

% Read the number of input values

nvals = input('Enter number of values: ');

% Read the numbers for ii = 1:nvals

x = input('Enter number: '); sumr = sumr + 1/x; end

% Calculate harmonic mean hmean = nvals / sumr;

% Tell user fprintf('The harmonic

mean

is %.4f\n',hmean); aa=input('The number:','s'); b=str2num(aa); n=length(b); prod=1; for ii=1:n

prod=prod*b(ii); end

9、12、13、14 in textbook P57. Reference answer:

chapter3_9

% Define variables:

% dist1AB -- Distance between A point and B point

% dist2BC -- Distance between B point and C point

% dist3AC -- Distance between A point and C point

% Ax1, Ay1, Az1 -- Point A

% Bx2, By2, Bz2 -- Point B

% Cx3, Cy3, Cz3 -- Point C

% Prompt the user for the input points

Ax1=input('Ax1: ');

Ay1=input('Ay1: ');

Az1=input('Az1: ');

Bx2=input('Bx2: ');

By2=input('By2: ');

Bz2=input('Bz2: ');

Cx3=input('Cx3: ');

Cy3=input('Cy3: ');

Cz3=input('Cz3: ');

% Calculate the distance

dist1AB=sqrt((Ax1-Bx2)^2+(Ay1-By2)^2+(Az1-Bz2)^2 );

dist2BC=sqrt((Bx2-Cx3)^2+(By2-Cy3)^2+(Bz2-Cz3)^2 );

dist3AC=sqrt((Ax1-Cx3)^2+(Ay1-Cy3)^2+(Az1-Cz3)^2 );

% Tell user

disp (['The distance AB is ' num2str(dist1AB)]); disp (['The distance BC is ' num2str(dist2BC)]); disp (['The distance AC is ' num2str(dist3AC)]);

chapter3_12 %Define variables:

x=-5:0.1:5;

y1=zeros(size(x));

y2=zeros(size(x));

y3=zeros(size(x));

N=length(x);

for k=1:N

if x(k)>1

y1(k)=1;

elseif x(k)<=1&x(k)>=-1

y2(k)=x(k);

else y3(k)=-1;

end

end

y=y1+y2+y3;

plot(x,y)

c hapter3_13

非循环

clc;

tic,t=cputime;

i=0:1:63;

k=sum(2.^i)

%disp(['k=1+2+2^2+2^3+...+2^62+2^63= '

num2str(k)]);

toc,cputime-t;

for循环

tic,t=cputime;

k=0;

for i=0:63

k=k+2.^i;

end

k

toc,cputime-t;

while 循环

tic,t=cputime;

sum=0;

i=0;

while 1

sum=sum+2.^i;

i= i+1;

if i> 63

break;

end

end

sum

toc,cputime-t;

chapter3_14

%Define variables:

% t--the equation involved

%x1--temp value of x in order to find the value of

x

%x2--temp value of y in order to find the value of x

%y1--temp value of x in order to find the value of y

%y2--temp value of y in order to find the value of y

%initialize t

t=input('please enter the value of the t :')

%get the initial values of x1 and x2

x1=input('please enter the value of the x1 :')

x2=input('please enter the value of the x2 :')

%count the value of x by loops

for i=1:32

y1=log(x1)-cos(x1+t);

y2=log(x2)-cos(x2+t);

x=0.5*(x1+x2);

y=log(x)-cos(x);

if y*y1>0

x1=x;

end

if y*y2>0

x2=x;

end

end

format long;

x,y

function [x,esterr]=equval(phi)

x1=0.2;x2=3;

for ii=1:50

y1=log(x1)-cos(x1+phi); y2=log(x2)-cos(x2+phi);

x=0.5*(x1+x2);

y=log(x)-cos(x+phi);

if y*y1>0, x1=x;end

if y*y2>0, x2=x;end

end

esterr=log(x)-cos(x+phi);

format long;

fprintf('\nWhen PHI is ');

disp(sym(phi));

fprintf('Value of x is %3.4f,',x);

fprintf('Error of the solution is %e.\n',esterr);

>>phi=0:pi/8:3*pi/8;

x=zeros(1,4);

for k=1:4

x(k)=equval(phi(k));

end

The output is:

When PHI is 0

Value of x is 1.3030,Error of the solution is

2.053913e-015.

When PHI is pi/8

Value of x is 1.0909,Error of the solution is

-9.575674e-016.

When PHI is pi/4

Value of x is 0.8957,Error of the solution is

-1.332268e-015.

When PHI is 3*pi/8

Value of x is 0.7230,Error of the solution is

-3.608225e-015.

一matlab基础知识

1 数值的表示

Matlab的数只采用习惯的十进制表示,可以带

小数点和负号;其缺省的数据类型为双精度浮点型

(double)。

例如:3 -10 0.001 1.3e10 1.256e-6

变量命令规则

变量名、函数名对字母的大小写是敏感的。如myVar

与myvar表示两个不同的变量。

变量名第一个字母必须是英文字母。

变量名可以包含下英文字母、划线和数字。

变量名不能包含空格、标点。

变量名最多可包含63个字符(6.5及以后的版本)

2 数组的分类

一维数组,也称为向量(vector) 。

行向量(row vector)、列向量(column vector)。

二维数组(矩阵matrix)。

多维数组。

有效矩阵:每行元素的个数必须相同,也就是说其每

列元素的个数耶也必须相同。

(1)创建一维数组变量

方法1 使用方括号“[ ]”操作符。例如A[1 2 3]

方法2 使用冒号“:”操作符例如B=1:10

方法3利用函数linspace 例如a=linspace(1 ,2,5)

方法4利用函数logspace

方法5 (3) 利用已有的矩阵剪裁方法创建向量。

列向量的创建

使用方括号“[ ]”操作符,使用分号“;”分割行

使用冒号操作符

创建数组变量的一般方法

创建变量的赋值语句的一般格式

var=expression

var为变量名

expression为MA TLAB合法表达式

可以是单独的常数值或数值数组;

也可以由常数值、其他变量(部分或全部)、数值数组和运算符(+、-等)构成。

创建二维数组变量

第一种方法:使用方括号“[ ]”操作符

“[ ]”操作符创建二维数组的使用规则

数组元素必须在“[ ]”内键入;

行与行之间须用分号“;”间隔,也可以在分行处用回车键间隔;

行内元素用空格或逗号“,”间隔。

函数方法

函数ones(生成全1矩阵)、zeros (生成全0矩阵) 、reshape

元素的寻访

a(3:-1:1) %由前3个元素倒序构成的子数组

修改a([2 5])=[1 1]

A(end,:) 表示最后一行;a(:,end)表示最后一列例子:a=[1 2 3 4 ;5 6 7 8]

A(:,end) %最后一列元素

4

8

a(end,:) %最后一行元素

5 6 7 8

a(:,end:-1:3)%从最后一列取整列到第三列

a([1 2])%按列开始从第一个元素取值,取两个元素1 5

a(end,[2:end-1])%最后一行从第2个元素开始取值至

该行倒数第2个

5

3 运算:算数关系逻辑

二matlab二维绘图

多次叠绘

多次调用plot来在一幅图上绘制多条曲线,需要hold

指令的配合。

hold on 保持当前坐标轴和图形,并可以接受下一次

曲线绘制。

hold off 取消当前坐标轴和图形保持,这种状态下,

调用plot绘制完全新的图形,不保留以前的坐标格式、

曲线。

subplot(m, n, k) 使(mxn)幅子图中第k个子图成为

当前图

A、分格线与grid指令

grid on 画出分格线

grid off 不画分格线

MATLAB的缺省设置是不画分格线;分格线的疏密取

决于坐标刻度(改变坐标刻度,可改变分格线的疏密)。

B、坐标框

box on 使当前坐标呈封闭形式(MATLAB的缺

省设置)

box off 使当前坐标呈开启形式

C、刻度设置

指令及格式:

set(gca, ‘xtick’, xs, ‘ytick’, ys)

xs、ys可以使任何合法的实数向量,用于分别设置x、

y轴的刻度。图形标识包括:

图名(title)

坐标轴名(xlabel、ylabel)

图形文本注释(text)

图例(legend)

六绘图

>> t=0:0.5:10;

>> y1=exp(-0.1*t);

>> y2=exp(-0.2*t);

>> y3=exp(-0.5*t);

>> plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')

>> title('\ity\rm=e^{-\itat}')

>> title('\ity\rm=e^{-\itat}','FontSize',12)

>>

text(t(6),y1(6),'\leftarrow\ita\rm=0.1','FontSize',11)

>>

text(t(6),y2(6),'\leftarrow\ita\rm=0.2','FontSize',11)

>>

text(t(6),y3(6),'\leftarrow\ita\rm=0.5','FontSize',11)

>> title('\ity\rm=e^{-\itat}','FontSize',12)

>> legend('a=0.1','a=0.2','a=0.5')

3 简述题

1 matlab 是一种什么样的软件,其主要特点有哪些?答:MA TLAB已经发展成为一个世界顶级的、功能强大的、易学易用的、综合数学工具软件。

它以高效的数组运算为特征;

它集符号运算与数值运算于一体,以数值计算见长;以函数形式内迁了许多优化的数学算法;

具有众多不同领域的专用工具箱;

不领域的系统建模、仿真、分析以及嵌入式系统的快速开发。

应用广泛(工程、科研、大学教学)

包括多个软件产品:

Matlab

Matlab专用工具箱(Toolbox)

Matlab编译器

Simulink

Simulink专用模块库(BlockSet)

语言特点:功能强大

数值运算优势

符号运算优势(Maple)

强大的2D、3D数据可视化功能

许多具有算法自适应能力的功能函数语言简单、内涵丰富

语言其书写形式非常接近于常规数学书写形式;

其操作和功能函数指令就是常用的计算机和数学书上

的一些简单英文单词表达的,如:help、clear等;

完备的帮助系统,易学易用。

扩充能力、可开发能力较强

MATLAB完全成了一个开放的系统

用户可以开发自己的工具箱

可以方便地与Fortran、C等语言接口

编程易、效率高

Matlab以数组为基本计算单元

具有大量的算法优化的功能函数

2 Matlab脚本程序文件的特点?函数程序文件的特

点?

M文件有两种形式:脚本文件(Script File)函数文

件(Function File )

这两种文件的扩展名,均为“. m”。M脚本文件

(Script File)

——实际上是一串MATLAB指令集合。

M脚本文件的特点:

脚本文件的构成比较简单,只是一串按用户意图排列

而成的(包括控制流向指令在内的)MATLAB指令集。

脚本文件运行后,所产生的变量都驻留在MATLAB

基本工作空间(Base workspace)中。只要用户不使用

清除指令(clear),MA TLAB指令窗不关闭,这些变

量将一直保存在基本工作空间中。

M函数文件的特点:

从形式上看,与脚本文件不同,函数文件的笫一行

总是以“function”引导的“函数申明行”。从运行上

看,与脚本文件运行不同,M函数中使用的变量都

是局部变量,即在该函数返回之后,这些变量会自动

在MA TLAB的工作空间中清除掉。

3 试举例查看matlab工作空间内变量的命名

6 数组运算与矩阵运算的运算符有什么区别?

在加、减运算时数组运算与矩阵运算的运算符相同,

乘、除和乘方运算时,在矩阵运算的运算符前加一个

点即为数组运算,如a*b为矩阵乘,a.*b为数组乘

一维数组相当于向量,二维数组相当于矩阵.所以矩阵

是数组的子集

数组运算是指数组对应元素之间的运算,也称点运算.

矩阵的乘法、乘方和除法有特殊的数学含义,并不是

数组对应元素的运算,所以数组乘法、乘方和除法的

运算符前特别加了一个点。

矩阵是一个二维数组,所以矩阵的加、减、数乘等运

算与数组运算是一致的。但有两点要注意:

(1)对于乘法、乘方和除法等三种运算,矩阵运算与数

组运算的运算符及含义都不同:矩阵运算按线性变换

定义,使用通常符号;数组运算按对应元素运算定义,

使用点运算符;

(2)数与矩阵加减、矩阵除法在数学是没有意义的,在

MATLAB中为简便起见,定义了这两类运算

数组运算:

转置 A.' 非共轭转置,相当于(conj(A'))

数组加与减A+B与A-B 对应元素之间加减

数乘数组k.*A或A.*k k乘A的每个元素

数与数组加减k+A与k-A k加(减)A的每个元素

数组乘数组A.*B

数组乘方A.^k A的每个元素进行k次方运算

k.^A 以k底的,分别以A的元素为指数求幂值

数除以数组k./A和A.\k k分别被B的元素除

数组除法左除A.\B右除B./A

矩阵运算:

矩阵转置A' 共轭转置

加减A+B A-B

数乘矩阵k*A或A*k 上三项同数组运算

矩阵乘法A*B 按数学定义的矩阵乘法规则

矩阵乘方A^k k个矩阵A相乘

数与矩阵加减k+A与k-A 等价于k*ones(size(A))+-A 矩阵除法左除A\B,右除B/A 分别为AX=B和XA=B的解

plot(x,y,string)

其中string 是用单引号括起来的字符串,用来指定图形的属性(点、线的形状和颜色)

mesh(X,Y,Z,C)

绘制由矩阵X,Y,Z 所确定的曲面网格图,矩阵 C 用于确定网格颜色,省略时C=Z。

mesh(Z)

绘出矩阵Z 的三维消隐图。

mesh(x,y,Z)

x, y是向量时,length(x)=n,length(y)=m,[m,n]=size(Z) [X,Y]= meshgrid(x,y)

x, y 为给定的向量,X, Y 是网格划分后得到的网格矩阵

绘制等高线meshc

绘制边界面屏蔽meshz

surf(X,Y,Z) 绘制由矩阵X,Y,Z 所确定的曲面图,参数含义同mesh

sphere(n) 专用于绘制单位球面

colormap: 设置绘图的着色方案

二维曲线绘图:ezplot

ezplot(f,[a,b]):绘制f=f(x)在区间(a

ezplot(f(x,y)):f(x,y) =0在区间(-2pi

ezplot(f,g,[a,b]):x=f(t),y=g(t) 在区间(a

空间曲线绘图:ezplot3

MATLAB编程作业

《Matlab 编程训练》 作业 专 业 学生姓名 班级 学 号 指导教师 完成日期

实训一 MATLAB 语言介绍和数值计算 1.先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存变量。 12 2sin851z e =+ . 2. 已知 1234413134787,2033657327A B --???? ????==???? ????-???? ,求下列表达式的值: (1) A+6*B 和A-B+I (其中I 为单位矩阵) A+6*B:

A-B+I: (2)A*B和A.*B A*B程序: A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B 结果: A.*B程序: A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果:

(3)A^3和A.^3 A^3程序: A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果: A.^3程序: A=[12 34 -4;34 7 87;3 65 7] C=A.^3 (4)A/B及B\A A/B程序: A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] C=A/B 结果:

B\A程序: A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=B\A 结果: (5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并(3)保存变量(mat文件)程序: A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; C=B*inv(A); D=C(2:3,2:3) 结果:

西南交大交通运输系统分析离线作业

交通运输系统分析第一次作业 三、主观题(共12道小题) 41.开放系统 答:系统与外界环境之间有物质、能量、信息交换的系统称为开放系统. 42.请举出一个系统的实例,并分析其目的、要素集、关系、结构和环境. 答:(主观题,自行举例即可) 43.系统的目标集中若出现相互矛盾的情况,应当怎样处理?请举例说明. 答:分目标之间可能是矛盾的,因此采用某种形式的折衷是必要的,即在矛盾的分目标之间寻求平衡.具体的做法是通过计算每个分目标对总目标的贡献来确定最佳的妥协. (自行举例既可) 44.切克兰德软系统方法论的步骤是什么? 答: (1)系统现状说明 通过调查分析,对现存的不良结构系统的现状进行说明. (2)弄清关联因素 初步弄清与现状有关的各种因素及其相互关系. (3)建立概念模型 在不能建立数学模型的情况下,用结构模型或语言模型来描述系统的现状.(4)改善概念模型 随着分析的不断深入和“学习”的加深,进一步用更合适的模型或方法改进上述概念模型. (5)比较 将概念模型与现状进行比较,找出符合决策者意图而且可行的改革途径或方案.(6)实施 实施所提出的改革方案. 45.交通运输系统的作业特征是什么? 答: (1)交通运输系统是一个连续过程系统; (2)交通运输系统生产的多环节、多功能、超区域的特点; (3)交通运输系统生产具有网络性特点; (4)交通运输系统是一个动态系统. 46.怎样理解交通运输系统的网络性? 答:交通运输生产不仅仅是列车、车辆在轨道、道路上移动,而且在交通运输网(包括铁路网、公路网、水运网、航空网等)上运动.良好的交通运输系统首先要有合理的布局与结构,要建设成与内部、外部协调的交通运输网.在具有科学

西南交通大学网络教育工程经济学作业.

2015—2016第2学期离线作业 科目: 姓名: 学号: 专业: 西南交通大学远程与继续教育学院 直属学习中心

工程经济学 第一次作业 主观题: 12.什么是沉没成本?什么是机会成本? 答:沉没成本是指由于过去的决策已经发生了的,而不能由现在或将来的任何决策改变的、已经计入过去投资费用回收计划的成本费用。 机会成本是指资源用于某用途而放弃其他用途可带来的最大收益。 13.什么是现值、终值、年值? 答:现值:资金发生在某一特定时间序列起点的价值,一般用P表示。 终值:资金发生在某一特定时间序列终点的价值,即期末的本利和,一般用F表示。 年值:各年(期)收入或支付的等额金额,一般用A表示。 14.常用的还本付息方式有哪几种?哪种方式最好? 答:还本付息方式有等额本息还款、等额本金还款、一次还本付息、每期付息到期一次还本4种方式。 因为资金具有时间价值,所以不能把发生在不同时点的利息或还款额直接相加来判断哪种还款方式更好。事实上,以上四种还款方式的资金时间价值完全一样。 15.何谓内部收益率?其评价标准如何?有何优缺点? 答:内部收益率是指能使项目计算期内净现金流量现值累计等于零的折现率

评价准则:若(或i s),则认为项目在经济上可以接受; 若(或i s),则项目在经济上不可接受,应当予以拒绝。 优点:在于考虑了资金的时间价值,而且考虑了项目在整个计算期内的经济状况;内部收益率是由项目自身的现金流量特征决定的,是内生的,所以在计算时不需要事先确定一个基准收益率。 缺点:一是内部收益率指标计算烦琐,而且需要大量的与投资项目有关的数据;二是对非常规项目,其内部收益率往往不唯一,在某些情况下也可能不存在;另外,内部收益率指标也不适用于只有现金流入或现金流出的项目。 16.什么是差额投资内部收益率?如何运用差额投资内部收益率进行互斥方案比选? 答:差额投资内部收益率(D IRR)是指两个互斥方案的差额投资净现值等于零时的折现率或两个方案净现值相等时的折现率。 差额投资内部收益率用备选方案差额现金流按下式计算: 差额投资内部收益率法的计算步骤如下: (1)将方案按投资额由小到大排序; (2)计算各方案的IRR(或NPV或NAV),淘汰IRR i c,则说明初始投资

MATLAB程序设计作业

Matlab程序设计 班级 姓名 学号

《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线,给出MATLAB代码和运行结果。 代码如下: x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下:

1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0,3π]区间,绘制y=sin(x)曲线(要求消去负半波,即(π,2π)区间内的函数值置零),求出曲线y 的平均值,以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 代码如下: clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1);

大学物理(西南交大)作业参考答案5

NO.5 电势、导体与※电介质中的静电场 (参考答案) 班级: 学号: 姓名: 成绩: 一 选择题 1.真空中一半径为R 的球面均匀带电Q ,在球心O 处有一带电量为q 的点电荷,如图所示,设无穷远处为电势零点,则在球内离球心O 距离为r 的P 点处的电势为: (A )r q 04πε; (B )(041 R Q r q +πε; (C )r Q q 04πε+; (D ))(0 41 R q Q r q -+ πε; 参考:电势叠加原理。 [ B ] 2.在带电量为-Q 的点电荷A 的静电场中,将另一 带电量为q 的点电荷B 从a 点移动到b ,a 、b 两点距离点电荷A 的距离分别为r 1和r 2,如图,则移动过程中电场力做功为: (A )(2 101 1Q --; (B )(2 101 14r r qQ -πε; (C ) )(2 1 114r r qQ --πε; (D ) ) (4120r r qQ --πε。 参考:电场力做功=势能的减小量。A=W a -W b =q(U a -U b ) 。 [ C ] 3.某电场的电力线分布情况如图所示,一负电荷从M 点移到N 点,有人根据这个图做出以下几点结论,其中哪点是正确的? (A )电场强度E M <E N ; (B )电势U M <U N ; (C )电势能W M <W N ; (D )电场力的功A >0。 [ C ] 4.一个未带电的空腔导体球壳内半径为R ,在腔内离球心距离为d (d <R )处,固定一电量为+q 的点电荷,用导线把球壳接地后,再把地线撤去,选无穷远处为电势零点,则球心O 处的点势为: (A )0; (B )d q 4πε; (C )-R q 04πε; (D ))(1 1 40 R d q - πε。 参考:如图,先用高斯定理可知导体内表面电荷为-q ,外表面无电荷(可分析)。虽然内表面电荷分布不均,但到O 点的距离相同,故由电势叠加原理可得。 [ D ] ※5.在半径为R 的球的介质球心处有电荷+Q ,在球面上均匀分布电荷-Q ,则在球内外处的电势分别为: (A )内r Q πε4+,外r Q 04πε-; (B )内r Q πε4+,0; 参考:电势叠加原理。注:原题中ε为ε0 (C )R Q r Q πεπε44-+内 ,0; (D )0,0 。 [ C ] r 2 (-Q)A b r 1 B a (q )

实例matlab-非线性规划-作业

实例matlab-非线性规划-作业

现代设计方法-工程优化理论、方法与设计 姓名 学号 班级 研 问题 : 某厂向用户提供发动机,合同规定,第一、二、三季度末分别交货40台、60台、80台。每季度的生产费用为 (元),其中x 是该季生产的台数。若交货后有剩余,可用于下季度交货,但需支付存储费,每台每季度c 元。已知工厂每季度最大生产能力为100台,第一季度开始时无存货,设a=50、b=0.2、c=4,问工厂应如何安排生产计划,才能既满足合同又使总费用最低。讨论a 、b 、c 变化对计划的影响,并作出合理的解释。 问题的分析和假设: 问题分析:本题是一个有约束条件的二次规划问题。决策变量是工厂每季度生产的台数,目标函数是总费用(包括生产费用和存储费)。约束条件是生产合同,生产能力的限制。在这些条件下需要如何安排生产计划,才能既满足合同又使总费用最低。 问题假设: 1、工厂最大生产能力不会发生变化; 2、合同不会发生变更; 3、第一季度开始时工厂无存货; 4、生产总量达到180台时,不在进行生产; 5、工厂生产处的发动机质量有保证,不考虑退货等因素; 6、不考虑产品运输费用是否有厂家承担等和生产无关的因素。 符号规定: x1——第一季度生产的台数; x2——第二季度生产的台数; 180-x1-x2——第三季度生产的台数; y1——第一季度总费用; y2——第二季度总费用; y3——第三季度总费用; y ——总费用(包括生产费用和存储费)。 ()2bx ax x f +=

建模: 1、第一、二、三季度末分别交货40台、60台、80台; 2、每季度的生产费用为 (元); 3、每季度生产数量满足40 ≤x1≤100,0≤x2≤100,100≤x1+x2 ≤180; 4、要求总费用最低,这是一个目标规划模型。 目标函数: y1 2111x b x a Z ?+?= y2()4012222-?+?+?=x c x b x a Z y3()()()10018018021221213 -+?+--?+--?=x x c x x b x x a Z y x x x x x x Z Z Z Z 68644.04.04.0149201 212221321--+++=++= 40≤x1≤100 0≤x2≤100 100≤x1+x2≤180 ()2 bx ax x f +=

西南交通大学2016 大物作业No.01

?西南交大物理系_2016_02 《大学物理AI 》作业 No.01运动的描述 班级 ________ 学号 ________ 姓名 _________ 成绩 _______ 一、判断题:(用“T ”和“F ”表示) 【 F 】1.做竖直上抛运动的小球,在最高点处,其速度和加速度都为0。 解:对于竖直上抛运动,在最高点,速度为0,加速度在整个运动过程中始终不变,为重力加速度。 【 F 】2.在直线运动中,质点的位移大小和路程是相等的。 解:如果在运动过程中质点反向运动,必然导致位移的大小和路程不相等。 【 F 】3.质点做匀加速运动,其轨迹一定是直线。 反例:抛体运动。 【 T 】4.质点在两个相对作匀速直线运动的参考系中的加速度是相同的。 解:两个相对作匀速直线运动的参考系的相对加速度为0,根据相对运动公式知上述说法正确。 【 F 】5.在圆周运动中,加速度的方向一定指向圆心。 反例:变速率的圆周运动。 二、选择题: 1.一小球沿斜面向上运动,其运动方程为2125t t S -+=(SI ),则小球运动到最高点的时刻应是 (A )s 4=t (B )s 2=t (C )s 12=t (D )s 6=t [ D ] 解:小球运动速度大小 t t s v 212d d -== 。 当小球运动到最高点时v =0,即 0212=-t ,t =6(s )。 故选 D 2. 一列火车沿着一条长直轨道运行,如图所示,曲线图显示了火车的位置时间关系。这个曲线图说明这列火车 [ B ] (A) 始终在加速 (B) 始终在减速 (C) 以恒定速度运行 (D) 部分时间在加速,部分时间在减速 解:位置时间曲线的上某点的切线的斜率就表示该时刻质点运动速度。由图可知,该火车一直在减速。

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

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 here theta=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

房屋建筑学(西南交大)离线作业

21.墙体设计要求有哪些? 22.提高墙体的热阻的措施有哪些 23.墙体的隔热措施有哪些 24.实体砖墙的组砌方式主要有哪几种 25.普通粘土砖墙的砖模尺寸与建筑模数是否一致?如何协调二者关系 26.勒脚的作用是什么?常用的做法有哪几种 27.墙身防潮层的作用是什么?水平防潮层的位置如何确定?什么请款下须设垂直发防潮层 28.墙身防潮层的作用是什么?水平防潮层的位置如何确定?什么请款下须设垂直发防潮层

29.窗洞口上部过梁的常用做法有哪几种,各自的适用范围如何 30.砖砌平拱过梁的构造要点是什么 31.圈梁的作用有哪些?设置原则主要有哪些 32.构造柱的作用及加设原则是什么? 33.构造柱的构造要点有哪些 ? 34.简述墙体三种变形缝的异同 35.简述 1/2 砖隔墙构造要点。 36.简述加气混凝土砌块隔墙构造要点 37.墙面装修的作用是什么 ? 38.什么是抹灰类墙面装修?有哪些构造层次?简单介绍各构造层次的作用与做法

39.什么是贴面类装修?常见贴面类装修有哪些? 40.什么是涂料类墙面装修?涂料施涂方法有哪些? 17.楼梯的作用及设计要求有哪些 18.楼梯主要由哪些部分组成?各部分的作用和要求是什么? 19.楼梯坡度的表达方式有那些? 20.当楼梯底层中间平台下做通道而平台净高不满足要求时,常采取哪些办法解决? 21.现浇钢筋混凝土楼梯有哪几种结构形式?各有何特点 22.预制踏步有哪几种断面形式和支承方式?

23.栏杆扶手在平行楼梯的平台转弯处如何处理? 24.室外台阶的构造要求是什么?通常有哪些做法。 25.电梯井道的构造要求有哪些? 26.简述楼梯的设计步骤。 27.< 影响基础埋置深度的因素有哪些? 28.基础按构造形式不同分为哪几种?各自的适用范围如何? 29.确定地下室防潮或防水的依据是什么? 30.地下室卷材外防水的层数是如何确定的? 21. 现浇钢筋混凝土楼板的特点和适用范围是什么

第一学期西南交大理论力学C第1次作业答案

本次作业是本门课程本学期的第1次作业,注释如下: 一、单项选择题(只有一个选项正确,共26道小题) 1. 考虑力对物体作用的运动效应和变形效应,力是。 (A) 滑动矢量 (B) 自由矢量 (C) 定位矢量 正确答案:C 解答参考: 2. 考虑力对物体作用的运动效应,力是。 (A) 滑动矢量 (B) 自由矢量 (C) 定位矢量 正确答案:A 解答参考: 3. 图示中的两个力,则刚体处于。 (A) 平衡 (B) 不平衡 (C) 不能确定 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:B 解答参考: 4.

作用力的大小等于100N,则其反作用力的大小为。 (A) (B) (C) 不能确定 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:B 解答参考: 5. 力的可传性原理只适用于。 (A) 刚体 (B) 变形体 (C) 刚体和变形体 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:A 解答参考: 6. 图示结构,各杆自重不计,则杆BC是。

(A) 二力杆 (B) 不能确定 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:A 解答参考: 7. 图示作用于三角架的杆AB中点处的铅垂力如果沿其作用线移动到杆BC的中点,那么A、C处支座的约束力的方向。 (A) 不改变 (B) 改变 (C) 不能确定 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:B 解答参考: 8.

图示构架ABC中,力作用在销钉C上,则销钉C对杆AC的作用力与销钉C对杆B C的作用力。 (A) 等值、反向、共线 (B) 分别沿AC和BC (C) 不能确定 你选择的答案:[前面作业中已经做正确] [正确] 正确答案:B 解答参考: 9. 如图所示,物体处于平衡,,自重不计,接触处是光滑的,图中所画受力图。 (A) 正确 (B) 不正确

BP神经网络matlab实例(简单而经典)

p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化 net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络 net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); %对BP网络进行仿真 anew=postmnmx(anewn,mint,maxt); %还原数据 y=anew'; 1、BP网络构建 (1)生成BP网络 = net newff PR S S SNl TF TF TFNl BTF BLF PF (,[1 2...],{ 1 2...},,,) R?维矩阵。 PR:由R维的输入样本最小最大值构成的2 S S SNl:各层的神经元个数。 [ 1 2...] { 1 2...} TF TF TFNl:各层的神经元传递函数。 BTF:训练用函数的名称。 (2)网络训练 [,,,,,] (,,,,,,) = net tr Y E Pf Af train net P T Pi Ai VV TV (3)网络仿真 = [,,,,] (,,,,) Y Pf Af E perf sim net P Pi Ai T {'tansig','purelin'},'trainrp'

西南交通大学大物A作业解析

?西南交大物理系_2013_02 《大学物理AI 》作业 No.03角动量 角动量守恒定律 班级 ________ 学号 ________ 姓名 _________ 成绩 _______ 一、判断题:(用“T ”和“F ”表示) [ F ] 1.如果一个刚体所受合外力为零,其合力矩一定为零。 [ F ] 2.一个系统的动量守恒,角动量一定守恒。 [ T ] 3.一个质点的角动量与参考点的选择有关。 [ F ] 4.刚体的转动惯量反映了刚体转动的惯性大小,对确定的刚体,其转动惯量是一定值。 [ F ] 5.如果作用于质点的合力矩垂直于质点的角动量,则质点的角动量将不发生变化。 二、选择题: 1.有两个半径相同、质量相等的细圆环A 和B 。A 环的质量分布均匀,B 环的质量分布不均匀。它们对通过环心并与环面垂直的轴的转动惯量分别为A J 和B J [ C ] (A) A J >B J (B) A J

实验二--MATLAB程序的设计(含实验报告)

实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台(带有MATLAB7.0以上的软件环境)。 M 文件的编写: 启动MATLAB 后,点击File|New|M-File ,启动MATLAB 的程序编辑及调试器(Editor/Debugger ),编辑以下程序,点击File|Save 保存程序,注意文件名最好用英文字符。点击Debug|Run 运行程序,在命令窗口查看运行结果,程序如有错误则改正 三、 实验容 1.编写求解方程02=++c bx ax 的根的函数(这个方程不一定为一元二次方程,因c b a 、、的不同取值而定),这里应根据c b a 、、的不同取值分别处理,有输入参数提示,当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 (提示:提示输入使用input 函数) 2.输入一个百分制成绩,要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+,90分~99分为A ,80分~89分为B ,70分~79分为C ,60分~69分为D ,60分以下为E 。 要求:(1)用switch 语句实现。 (2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。 (提示:注意单元矩阵的用法) 3.数论中一个有趣的题目:任意一个正整数,若为偶数,则用2除之,若为奇数,则与3相乘再加上1。重复此过程,最终得到的结果为1。如: 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1

西南交通大学 道路工程 大作业

道路工程 姓名:学号: 1.基本资料 (1)设计任务书要求 我国某公路部分路段设计等级为高速公路,设计年限为15年,拟采用沥青路面结构,需进行路面结构设计。 (2)气象资料 该公路地处II2区,年降水量为620mm/年,最高气温35℃,最低气温-31℃,多年平均冻结指数为882℃·d极大冻结指数为1225℃·d。 (3)地质资料与筑路材料 沿线土质为中液限粘性土,填方路基高度1.8m,地下水位距路床底2.4m,一般路基处于中湿状态。公路沿线有丰富的砂砾,附近有小型采石场和石灰厂,筑路材料丰富。路面所用水泥和沥青均需外购。 (4)交通资料 据工程可行性报告得知近期交通组成与交通量如表所示;交通量年增长率取7%;不同车型的交通参数如表所示。

2.计算书和设计图纸如下所示: (1)初步设计两种沥青路面结构 本题是一条设计年限为15年的高速公路,并且要求道路路面面层结构使用沥青混凝土。由于该条道路施工沿线有丰富的砂砾、小型采石场和石灰厂。考虑道路施工石灰材料丰富,所以道路基层设计时使用石灰粉煤灰碎石材料,减小筑路材料的运输成本。根据沥青道路的一般设计方案以及本人的工作经验,现初步设计两种沥青路面结构:方案一的道路面层使用细粒式沥青混凝土,设计厚度为4cm;中粒式沥青混凝土,设计厚度为6cm;粗粒式沥青混凝土,设计厚度为7cm;道路基层使用18cm 水泥混凝土碎石和二灰土。二灰土的厚度根据题目给定的参数计算确定。方案一的道路结构图纸如下所示:

方案二的道路面层使用细粒式沥青混凝土,设计厚度为4cm;中粒式沥青混凝土,设计厚度为6cm;大粒径沥青碎石(ATB),设计厚度为10cm;道路基层使用二灰稳定砂粒和20cm的级配碎石。二灰稳定砂粒的厚度根据题目给定的参数计算确定。方案二的道路结构图纸如下所示: (2)交通轴载换算 根据题目参数表格中给出的交通参数进行荷载换算,如下所示:

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4.设???? ??????------=81272956313841A ,??????????-----=793183262345B ,求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5.设 ?? ????++=)1(sin 35.0cos 2x x x y ,把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6.产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 (mean var ) a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)

西南交通大学机械制造基础作业及答案

& 机械制造基础线下作业 第一次作业 一、主观题(共8道小题) 1.车外圆时,已知工件转速n=320r/min,车刀移动速度v f=64mm/min,试求切削速度v c,进给量f, 背吃刀量a p。已知待加工表面直径100mm, 已加工表面直径94mm. 答:V=s f=s ap=3mm 2.已知下列车刀的主要角度,试画出它们切削部分的示意图。 % 1)端面车刀:γ0=15°α0=10°κr=45°κr’=45° 2)切断刀:γ0=10°α0=6°κr=90°κr’=20° 3)外圆车刀:γ0=15°α0=8°κr=75°κr’=15° 答:(1)(2) (3) 3.实现切削加工的必要条件是什么 : 答:切削运动,刀具材料,刀具角度 加工方法机床主运动进给运动 ; 车 卧式车床立式车床工件的旋转运动刀具的纵向或横向进给 钻台式钻床立式钻床 摇臂钻床刀具的旋转运动¥ 刀具的轴向进给 镗卧式镗床(坐标镗床)刀具的旋转运动工作台纵向进给或横向进给 刀具径向或轴向进给 铣立式铣床、卧式铣床》 刀具的旋转运动 工作台纵向、横向和垂直进给 刨牛头刨床 龙门刨床刀具的往复移动 工作台的往复移动 工作台横向进给 !

5. 刀具常用三种刀具材料(高速钢、硬质合金及碳素工具钢)的特点及适用条件是什么。 ] 6. 切削三要素包括那些各自的定义是什么 答: 切削速度, 进给量, 切削深度, 切削速度:在单位时间内,刀具(或工件)沿主运动方向的相对位移. 进给量:主运动在每行程(或每转)内,工件与刀具之间沿进给运动方向的相对位移. 切削深度:工件已加工表面与待加工表面间的垂直距离,即在垂直于进给运动的方向上测量的主切削刃切入工件的深度,单位 mm 。 % 7. 刀具角度的参考坐标系有哪三个平面组成如何定义 答: 基面, 切削平面, 主剖面 基面: 通过切削刃某选定点, 与主运动假定方向相垂直的平面。 切削平面: 通过刀刃某选定点, 包含主切削刃且与基面相垂直的平面。 主剖面(正交平面) 通过刀刃某选定点, 同时垂直于基面与切削平面的平面 : 8. 刀具的基本角度有哪些它们是如何定义的 答: 前角, 后角, 主偏角, 副偏角, 刃倾角 前角: 主剖面内, 前刀面与基面之间的夹角。 后角: 主剖面内, 后刀面与切削平面之间的夹角。 主偏角: 基面内, 主切削刃与进给运动方向间的夹角。 副偏角: 基面内, 副切削刃与进给运动反方向间的夹角。刃倾角: 切削平面内, 主切削刃与基面间的夹角

matlab程序设计作业

Matlab程序设计作业 姓名: 学号: 专业:

? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线,给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10

7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合

x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在[0, 3n区间,绘制y二Sin(x)曲线(要求消去负半波,即(n 2n)区间内的函数值置零),求出曲线y 的平均值,以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >>

matlab函数计算的一些简单例子1

MATLAB作业一1、试求出如下极限。 (1) 23 25 (2)(3) lim (5) x x x x x x x ++ + →∞ ++ + ,(2) 23 3 1 2 lim () x y x y xy x y →- → + + ,(3) 22 22 22 1cos() lim ()x y x y x y x y e+ → → -+ + 解:(1)syms x; f=((x+2)^(x+2))*((x+3)^(x+3))/((x+5)^(2*x+5)) limit(f,x,inf) =exp(-5) (2)syms x y; f=(x^2*y+x*y^3)/(x+y)^3; limit(limit(f,x,-1),y,2) =-6; (3)syms x y; f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2); limit(limit(f,x,0),y,0) =0 2、试求出下面函数的导数。 (1 )() y x=, (2)22 atan ln() y x y x =+ 解; (1)syms x; f=sqrt(x*sin(x)*sqrt(1-exp(x))); g= diff(f,x); g== (sin(x)*(1 - exp(x))^(1/2) + x*cos(x)*(1 - exp(x))^(1/2) - (x*exp(x)*sin(x))/(2*(1 - exp(x))^(1/2)))/(2*(x*sin(x)*(1 - exp(x))^(1/2))^(1/2)) pretty(g)= (2)syms x y; f=atan(y/x)-log(x^2+y^2) pretty(-simple(diff(f,x)/diff(f,y)))= 2 x + y =------- x - 2 y (3) 假设1 cos u- =,试验证 22 u u x y y x ?? = ???? 。 解:syms x y; u=1/cos(sqrt(x/y)); diff(diff(u,x),y)-diff(diff(u,y),x)=0; 所以: 22 u u x y y x ?? = ????

大学英语IV 西南交大 离线作业

大学英语IV第一次作业: (一)41. Directions: For this part, you are allowed thirty minutes to write a composition on the topic “My Ideal Job”. You should write no less than 100 words. My Ideal Job Every person has his own ideal job. Some people wish to be doctors, while others want to be teachers and so on. As far as I am concerned, I have been dreaming of being a guide for a long time. First, my major is the manage of tourism, so I want to get a job about tourism. Secondly, good guides are knowledgeable and amiable, meanwhile, they have wide vision. Finally, guide can constantly expand his knowledge and understand lots of customs about different cultures and countries. To achieve my dream of being a qualified guide, I have to make sufficient preparations. Fro one thing, I must learn my specialty well and travel at my spare time, just as the saying goes: “Seeing is believing.” For another, I should enhance my eloquence so that I can express the allusion about the view clearly. What’s more important ,I should control my mood freely and treat every traveler kindly. Being a good guide is difficult, but I believe if I endeavor to do this things in all hands. I’ll be successful sooner or later. 大学英语IV第二次作业: (二)41. Directions: For this part, you are allowed thirty minutes to write a composition on the topic “My Favorite Means of Getting Information”. You should write no less than 100 words.

matlab仿真实例

matlab 仿真实例 实验五MATLAB 及仿真实验一、控制系统的时域分析 (一)稳定性 1、系统传递函数为G(s),试判断其稳定性。 程序: >> nu m=[3,2,5,4,6]; >> den=[1,3,4,2,7,2]; >> sys=tf( nu m,de n); >> figure(1); >> pzmap(sys); >> title(' 零极点图') 由图可知:在S 右半平面有极点,因此可知系统是不稳定的。 2、用MATLA 求 出 G(s)=(s A 2+2*s+2)/(s A 4+7*s A 3+5*s+2) 的极点。 程序及结果: >> sys=tf([1,2,2],[1,7,3,5,2]); >> p=pole(sys) 矿'. 赳 _ ■ —

-6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i -0.4100 (二)阶跃响应 1、二阶系统G(s)=10/s A2+2*s+10 1)键入程序,观察并记录单位阶跃响应曲线: 程序: >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线') 2)计算系统闭环跟、阻尼比、无阻尼振荡频率,并记录程序及结果: >> sys=tf(10,[1,2,10]); >> p=pole(sys)

p = -1.0000 + 3.0000i -1.0000 - 3.0000i >> [wn,z]=damp(sys) wn = 3.1623 3.1623 z = 0.3162 0.3162 3)记录实际测取的峰值大小,峰值时间和过渡过程时间,并填表实际值理论值峰值Cmax 1.35s 峰值时间tp 1.05s 过渡时间+5% 3.54s ts +2% 3.18s 程序: >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线')

西南交通大学信号处理期末作业

1、考虑两个谐波信号()x t 和()y t ,其中()cos()c x t A w t φ=+,()cos()c y t B w t =式中A 和c w 为正的常数,φ为均匀分布的随机变量,其概率密度函数为 1 ,02()20,f φπ φπ ?≤≤?=??? 其他, 而B 是一个具有零均值和单位方差的标准高斯随机变量,即其分布函数为 (1)求()x t 的均值()x u t 、方差2()x t σ、自相关函数()x R τ和自协方差函数()x c τ。 (2)若φ与B 为相互统计独立的随机变量,求()x t 和()x t 的互相关函数()xy R τ与互协方差函数()xy c τ。 解: (1) ()x t 的均值()x u t 为: 方差2()x t σ为: 自相关函数()x R τ为: 自协方差函数()x c τ为: (2)()y t 的均值为: ()(())(cos())()cos()0y B B c c u t E y t E B w t E B w t ====,所以()=0E B 由互相关函数的定义可知: 由题意知道φ与B 为相互统计独立的随机变量,所以有 互协方差函数()xy c τ 2.接收信号由下式给出:cos(),1,2,...,i c i y A i i N ωθω=++=,式中~(0,1)i N ω即i ω是零均值和单位方差的高斯噪声,c ω为载波角频率,而θ是未知的相位。假设12,,...N ωωω相互独立,求未知相位的最大似然估计^ ML θ。 解:由于12,,...N ωωω相互独立,所以1,..N y y 也相互独立并且服从高斯分布,可以得到1,..N y y 与θ的联合概率密度函数分布 由此,可以得到似然函数 该似然函数对θ求偏导,并令该偏导函数为零,即可得到如下公式: 因此,最大似然估计^ ML θ 为上述函数的零点值。 则 该函数为非线性方程,不容易求解,若忽略双倍频率2c ω ,则可简化到如下式子: 根据三角公式分解得到如下式子: 由此,可以得到如下公式 所以相位的最大似然估计如下: 3.离散时间的二阶AR 过程由差分方程12()(1)(2)()x n a x n a x n w n =-+-+ 描述,式中()w n 是 一零均值、方差为2 w σ 的白噪声。证明()x n 的功率谱为 证明: 由AR 过程的功率谱公式知 其中

相关文档 最新文档