东南大学matlab课程作业

Matlab Worksheet 1

Part A

1. Get into Matlab: Use the diary command to record your activity in to a file: diary mydiary01.doc before you start your work. (And diary off to switch off your diary when you finish your work.)

At the Command Window assign a v alue x=10, then use the Up Key ↑ to repeat the expression, editing it to show the effect of using a semicolon after the 10, namely x=10;

Answers:

>> x=10

x =

>> x=10;

2. Confirm whether the following names of variables are acceptable: a) Velocity Yes No

b) Velocity1 Yes No

c) Velocity.1 Yes No

d) Velocity_1 Yes No

e) Velocity-01 Yes No

f) velocityONE Yes No

g) 1velocity Yes No

3. Assign two scalar variables x and y values x=1.5, y=2.3, then construct Matlab expressions for the following: a) y x y x z +=35 b) ()3/27y x z = c) ()

y e x x z 25106

/111log -???????+-= d) )2cos()2sin(x e y x z x ππ+-=

Answers:

>> x=1.5;y=2.3

z1=5*x^3*y/(x+y)

z2=(x^7*y^0.5)^(2/3)

z3=(x^(1/6)/(log10(x^5-1))+1)*exp(-2*y)

z4=sin(2*pi*x-y)+exp(x).*cos(2*pi*x)

z1 =

10.2138

z2 =

8.7566

z3 =

0.0232

z4 =

-3.7360

4. Assign two variables with complex values u=2+3j and v=4+10j and then construct expression for:

a) v u b) j uv 2+ c) j

u 2 d) u j ve π2- Answers:

>> u=2+3j;v=4+10j;

z1=u/v

z2=u*v+2j

z3=u/2j

z4=v*exp(-j*2*pi*u)

z1 =

0.3276 - 0.0690i

z2 =

-22.0000 +34.0000i

z3 =

1.5000 - 1.0000i

z4 =

6.1421e+08 + 1.5355e+09i

https://www.360docs.net/doc/ee14925940.html,e the colon operator : to assign numerical values between 0 and 1 to vector array

variable a in steps of 0.1.

Answer:

>> V=0:0.1:1

V =

Columns 1 through 10

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000

0.9000

Column 11

1.0000

https://www.360docs.net/doc/ee14925940.html,e linspace function to assign 20 values to vector variable y between 20 and 30.

Answer:

>> V=linspace(20,30,20)

V =

Columns 1 through 10

20.0000 20.5263 21.0526 21.5789 22.1053 22.6316 23.1579 23.6842

24.2105 24.7368

Columns 11 through 20

25.2632 25.7895 26.3158 26.8421 27.3684 27.8947 28.4211 28.9474

29.4737 30.0000

7.Assign 20 values to a variable h increasing logarithmically between 10 and 1000.

Next, use the colon operator to assign the first 10 elements of h to a variable p.

Answers:

>> v=logspace(1,3,20)

v =

1.0e+03 *

Columns 1 through 10

0.0100 0.0127 0.0162 0.0207 0.0264 0.0336 0.0428 0.0546

0.0695 0.0886

Columns 11 through 20

0.1129 0.1438 0.1833 0.2336 0.2976 0.3793 0.4833 0.6158

0.7848 1.0000

>> p=v(1:10)

p =

10.0000 12.7427 16.2378 20.6914 26.3665 33.5982 42.8133 54.5559 69

.5193

88.5867

8.Create 6 element row vector z with values 1.0 1.2 1.6 -1.7 1.8 1.9, then construct

an expression for the sum of the 2nd 4th and 6th elements of z.

Answers:

>> z=[1.0 1.2 1.6 -1.7 1.8 1.9]

x=z(2)+z(4)+z(6)

z =

1.0000 1.2000 1.6000 -1.7000 1.8000 1.9000

x =

1.4000

https://www.360docs.net/doc/ee14925940.html,e the colon operator to create a vector array x between 10 and -10in steps of -1,

and second, an array vector y between 20 and -20 in steps -2 .

a)Add x and y, b) subtract 10x from 5y.

Answers:

>> x=10:-1:-10

y=20:-2:-20

x =

Columns 1 through 17

10 9 8 7 6 5 4 3 2 1 0 -1 -2 -3 -4 -5 -6

Columns 18 through 21

-7 -8 -9 -10

y =

Columns 1 through 17

20 18 16 14 12 10 8 6 4 2 0 -2 -4 -6 -8 -10 -12

Columns 18 through 21

-14 -16 -18 -20

>> x+y

ans =

Columns 1 through 17

30 27 24 21 18 15 12 9 6 3 0 -3 -6 -9 -12 -15 -18

Columns 18 through 21

-21 -24 -27 -30

>> 5*y-10*x

ans =

Columns 1 through 17

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Columns 18 through 21

0 0 0 0

https://www.360docs.net/doc/ee14925940.html,e the size,length, who and whos commands to establish the size and length of x

and y from Question 9, and use transpose operator ’ to convert vector x from

Question 9.

Answers:

>> size(x)

length(x)

ans =

1 21

ans =

>> who

Your variables are:

ans x y

>> whos

Name Size Bytes Class Attributes

ans 1x1 8 double

x 1x21 168 double

y 1x21 168 double

>> z=x'

z =

-1

-2

-3

-4

-5

-6

-7

-8

-9

-10

11.Show that if w=[ 2i 3i 3+i] the .’ operator creates the transpose. What effect does the

operator ’ applied to w have on its own?

Answer:

>> w=[ 2i 3i 3+i]

w =

0 + 2.0000i 0 + 3.0000i 3.0000 + 1.0000i

>> x=w'

x =

0 - 2.0000i

0 - 3.0000i

3.0000 - 1.0000i

>> y=w.'

y =

0 + 2.0000i

0 + 3.0000i

3.0000 + 1.0000i

https://www.360docs.net/doc/ee14925940.html,e the ones function to create a 4 by 6 array of 1’s. Considering just the shape of

the resulting array, what do the expression ones(3), ones(5) and ones(7) all have in common?

Answer:

>> v=ones(4,6)

v =

1 1 1 1 1 1

>> ones(3)

ans =

1 1 1

>> ones(5)

ans =

1 1 1 1 1

Ones(n)是n阶全为1的方阵

13.Create using the rand function a 5 by 4 random matrix and assign it to matrix array

variable A and observe carefully what A(:,3) A(1:2) and A(3,[2 4]) mean.

Answer:

>> A=rand(5,4)

A =

0.8147 0.0975 0.1576 0.1419

0.9058 0.2785 0.9706 0.4218

0.1270 0.5469 0.9572 0.9157

0.9134 0.9575 0.4854 0.7922

0.6324 0.9649 0.8003 0.9595

>> A(:,3)

ans =

0.1576

0.9706

0.9572

0.4854

0.8003

>> A(1:2)

ans =

0.8147 0.9058

>> A(3,[2 4])

ans =

0.5469 0.9157

https://www.360docs.net/doc/ee14925940.html,ing array subscripts, create an expression for the sum of the element in the top

right-hand-corner of A and the bottom left-hand-corner of A. Also assign the 2nd column of A to a column vector b, and assign the 3rd row of A to a row vector d.

Answer:

>> A=rand(5,4)

A =

0.7513 0.9593 0.8407 0.3500

0.2551 0.5472 0.2543 0.1966

0.5060 0.1386 0.8143 0.2511

0.6991 0.1493 0.2435 0.6160

0.8909 0.2575 0.9293 0.4733

>> x=A(1,4)+A(4,1)

x =

1.0491

>> b=A(:,2)

b =

0.9593

0.5472

0.1386

0.1493

0.2575

>> d=A(3,:)

d =

0.5060 0.1386 0.8143 0.2511

>>diary off to switch off your diary now.

15. Using the colon operator, assign a row vector array t, values between 0 and 10 in steps of 0.01. Use the ; operator to prevent displaying the information. Obtain the term-by-term values of functions:

a) t e z 05.0-= b) )sin(05.0t e z t π-= c) )2cos()sin(05.0t t e z t ππ-=

Answers:

0.5874 0.6186 0.6473 0.6731 0.6958 0.7153 0.7314 0.7439 0.7528

0.7581

Columns 551 through 560

0.7596 0.7573 0.7513 0.7417 0.7284 0.7117 0.6916 0.6684 0.6421

0.6131

Columns 561 through 570

0.5815 0.5476 0.5117 0.4741 0.4350 0.3948 0.3538 0.3123 0.2706

0.2291

Columns 571 through 580

0.1880 0.1477 0.1085 0.0706 0.0344 0.0000 -0.0322 -0.0621

-0.0895 -0.1141

Columns 581 through 590

-0.1359 -0.1548 -0.1705 -0.1832 -0.1928 -0.1992 -0.2025 -0.2027 - 0.2000 -0.1944

Columns 591 through 600

-0.1861 -0.1753 -0.1621 -0.1467 -0.1295 -0.1105 -0.0901 -0.0686 -

0.0462 -0.0232

Columns 601 through 610

0.0000 0.0232 0.0461 0.0684 0.0898 0.1099 0.1287 0.1457 0.1608 0.1737

Columns 611 through 620

0.1843 0.1923 0.1976 0.2001 0.1997 0.1962 0.1897 0.1801 0.1675 0.1518

Columns 621 through 630

0.1332 0.1117 0.0875 0.0607 0.0315 -0.0000 -0.0335 -0.0687 -

0.1055 -0.1435

Columns 631 through 640

-0.1824 -0.2221 -0.2621 -0.3022 -0.3420 -0.3813 -0.4197 -0.4569 -

0.4927 -0.5267

Columns 641 through 650

-0.5587 -0.5885 -0.6157 -0.6403 -0.6619 -0.6804 -0.6957 -0.7076 -

0.7161 -0.7211

Columns 651 through 660

-0.7225 -0.7204 -0.7147 -0.7055 -0.6929 -0.6770 -0.6579 -0.6358 -

0.6108 -0.5832

Columns 661 through 670

-0.5532 -0.5209 -0.4868 -0.4510 -0.4138 -0.3756 -0.3366 -0.2971 -

0.2574 -0.2179

Columns 671 through 680

-0.1788 -0.1405 -0.1032 -0.0672 -0.0327 -0.0000 0.0307 0.0591

0.0851 0.1085

Columns 681 through 690

0.1293 0.1472 0.1622 0.1743 0.1834 0.1895 0.1926 0.1928 0.1902 0.1849

Columns 691 through 700

0.1771 0.1667 0.1542 0.1396 0.1231 0.1051 0.0857 0.0652 0.0439 0.0221

Columns 701 through 710

0.0000 -0.0221 -0.0439 -0.0650 -0.0854 -0.1046 -0.1224 -0.1386 - 0.1530 -0.1653

Columns 711 through 720

-0.1753 -0.1829 -0.1880 -0.1903 -0.1899 -0.1866 -0.1805 -0.1713 - 0.1593 -0.1444

Columns 721 through 730

-0.1267 -0.1063 -0.0832 -0.0577 -0.0299 0.0000 0.0318 0.0654

0.1003 0.1365

Columns 731 through 740

0.1735 0.2113 0.2493 0.2874 0.3253 0.3627 0.3992 0.4346 0.4686 0.5010

Columns 741 through 750

0.5315 0.5598 0.5857 0.6090 0.6296 0.6472 0.6618 0.6731 0.6812

0.6859

Columns 751 through 760

0.6873 0.6853 0.6798 0.6711 0.6591 0.6440 0.6258 0.6048 0.5810 0.5548

Columns 761 through 770

0.5262 0.4955 0.4630 0.4290 0.3936 0.3573 0.3201 0.2826 0.2449 0.2073

Columns 771 through 780

0.1701 0.1336 0.0981 0.0639 0.0311 0.0000

-0.0292 -0.0562 -

0.0809 -0.1033

Columns 781 through 790

-0.1230 -0.1400 -0.1543 -0.1658 -0.1744 -0.1802 -0.1832 -0.1834 -

0.1810 -0.1759

Columns 791 through 800

-0.1684 -0.1586 -0.1467 -0.1328 -0.1171 -0.1000 -0.0815 -0.0621 -

0.0418 -0.0210

Columns 801 through 810

0.0000 0.0210 0.0417 0.0619 0.0812 0.0995 0.1164 0.1318 0.1455 0.1572

Columns 811 through 820

0.1667 0.1740 0.1788 0.1811 0.1807 0.1775 0.1717 0.1630 0.1516 0.1374

Columns 821 through 830

0.1205 0.1011 0.0792 0.0549 0.0285 0.0000 -0.0303 -0.0622 -

0.0954 -0.1298

Columns 831 through 840

-0.1651 -0.2010 -0.2372 -0.2734 -0.3094 -0.3450 -0.3797 -0.4134 - 0.4458 -0.4766

Columns 841 through 850

-0.5055 -0.5325 -0.5571 -0.5793 -0.5989 -0.6157 -0.6295 -0.6403 - 0.6480 -0.6525

Columns 851 through 860

-0.6538 -0.6518 -0.6467 -0.6384 -0.6270 -0.6126 -0.5953 -0.5753 - 0.5527 -0.5277

Columns 861 through 870

-0.5005 -0.4714 -0.4405 -0.4081 -0.3744 -0.3398 -0.3045 -0.2688 - 0.2329 -0.1972

Columns 871 through 880

-0.1618 -0.1271 -0.0934 -0.0608 -0.0296 -0.0000 0.0277 0.0535

0.0770 0.0982

Columns 881 through 890

0.1170 0.1332 0.1468 0.1577 0.1659 0.1714 0.1743 0.1745 0.1721 0.1673

Columns 891 through 900

0.1602 0.1509 0.1395 0.1263 0.1114 0.0951 0.0776 0.0590 0.0398 0.0200

Columns 901 through 910

0.0000 -0.0200 -0.0397 -0.0589 -0.0773 -0.0946 -0.1108 -0.1254 - 0.1384 -0.1495

Columns 911 through 920

-0.1586 -0.1655 -0.1701 -0.1722 -0.1718 -0.1689 -0.1633 -0.1550 - 0.1442 -0.1307

Columns 921 through 930

-0.1147 -0.0962 -0.0753 -0.0522 -0.0271 -0.0000 0.0288 0.0591

0.0908 0.1235

Columns 931 through 940

0.1570 0.1912 0.2256 0.2601 0.2944 0.3281 0.3612 0.3932 0.4240 0.4533

Columns 941 through 950

0.4809 0.5065 0.5300 0.5511 0.5697 0.5856 0.5988 0.6091 0.6164 0.6207

Columns 951 through 960

0.6219 0.6200 0.6151 0.6072 0.5964 0.5827 0.5663 0.5472 0.5257 0.5020

Columns 961 through 970

0.4761 0.4484 0.4190 0.3882 0.3562 0.3233 0.2897 0.2557 0.2216 0.1876

Columns 971 through 980

0.1539 0.1209 0.0888 0.0578 0.0281 0.0000

-0.0264 -0.0509 -

0.0732 -0.0934

Columns 981 through 990

-0.1113 -0.1267 -0.1396 -0.1500 -0.1578 -0.1631 -0.1658 -0.1660 -

0.1637 -0.1592

Columns 991 through 1000

-0.1524 -0.1435 -0.1327 -0.1201 -0.1060 -0.0905 -0.0738 -0.0562 -

0.0378 -0.0190

Column 1001

-0.0000

16.For Question 15 c), use the plot function to plot the value z against t. And use

xlabel and ylabel to label the t axis ‘time’ and z axis ‘Response’ and use title to give your figure a title. Copy the figure into your mydiary01.doc (from the figure edit menu, select copy, then open mydiary01.doc and paste in the figure as an object).

Answers:

>> plot(t,z3,'m','LineWidth',2)

xlabel('time')

ylabel('Response')

title('z=exp(-0.05*t)*sin(πt)*cos(2πt))')

Part B

1.Create a new Script m-file, called myfirst.m, to assign values to three vector arrays

a, b, and c; add each component together and assign to array d, where

A=[ 1, 2, 3], b=[4 5 6], c=[7 8 9]. Run myfirst.m. Examine the variable values in the WorkSpace . Edit myfirst.m to change variable b to b=[10 20 30]. Re-run myfirst.m

Answers:

myfirst.m

A=[ 1, 2, 3];

B=[4 5 6];

C=[7 8 9];

D=A+B+C

result:

>> myfirst

D =

12 15 18

myfirst.m

A=[ 1, 2, 3];

B=[10 20 30];

C=[7 8 9];

D=A+B+C

(2-2)result:

>> myfirst

D =

18 30 42

2.Create a function m-file, called myfun.m, that accepts two arrays x and y as

arguments, and produces the sum of x and y, and gives as the function output, the maximum of x+y.

Answer:

(1)myfun.m

function

[M] = myfun(x,y)

%MAXIMUM

sum=x+y;

t=length(x)-1;

for

i=1:t;

sum(1+i)>sum(i)

max=sum(1+i);

sum(1+i)=sum(i);

end

M=max

end

(2)result

>> x=[1 2 3]

y=[4 5 6]

x =

1 2 3

y =

4 5 6

>> myfun(x,y)

M =

ans =

3.Modify your script m-file myfirst.m by calling in it myfun.m twice, namely with a

and b as parameters in the argument list, and second with parameters b and c. Run myfirst.mand afterwards. Examine the WorkSpace variable values.

Answers:

matlab 作业

实验一 1、熟悉MATLAB的窗口结构（命令窗口、历史命令窗口、工作区窗口、当前目录窗口） 2、掌握命令窗口中基本命令的使用 3、在命令窗口中，给定圆的半径r，求得圆的周长c和面积s，并查看工作区窗口的变化 4、将r,c,s变量保存到磁盘文件abc.mat中，并删除内存变量r,c,s，查看工作区窗口的变化 5、将abc.mat文件中变量装入内存，查看工作区窗口的变化 6、将历史命令窗口中的命令再装入命令窗口中使用 7、改变当前目录，查看当前目录窗口的变化 8、掌握命令窗口中 cd,quit,help,date,dir,ls,what,who,clocl,fix(clock),format,save,loa d,clc,clear等命令的使用 9、注意各种MATLAB版本的差别实验二 1、在命令窗口中，输入长方形的长和宽，求长方形的周长和面积 2、输入三角形的三条边（要满足构成三角形的条件），求三角形的周长和面积 3、掌握MATLAB中各标准函数的使用（sin,cos,sind,fix,mod,…） 4、用fprintf输出各种类型的数据（如fprintf('a=%d\n',123) a=123 >> fprintf('b=%f\n',123.456) b=123.456000 >> fprintf('c=%c\n','A') c=A……) 实验三 1、在编辑窗口中：输入学生成绩，输出该成绩的等级。等级规定如下：[90， 100]为A等，[80，90)为B等，[70，80)为C等，[60，70)为D等，[0，60)为E等。要求用if和 switch两种方法实现。 2、商场购物，100件以下，不优惠，100~199件95折，200~399件90折，400~799 件85折，800~1499件80折，1500件以上，75折。输入所购货物的单价、件数，求实际付款数目。要求用if和 switch两种方法实现（在编辑窗口中实现）。实验四 1、求两个正整数的最大公约数和最小公倍数（在编辑窗口中实现，命令窗口中调用）。 2、求100~300内所有素数（在编辑窗口中实现，命令窗口中调用）。

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 所示。系统模型设计模型分析模型执行模型计算机图1 从模型设计到模型分析经历的过程，即对实物进行仿真可概括为以下几个方面： 1) 问题的描述; 2) 建立概念模型; 3) 建立仿真模型; 4) 收集数据; 5) 编写程序; 6) 在计算机进行模型试验; 7) 模型和数据的验证; 8) 仿真结果显示; 9) 仿真结果分析和评估。仿真是建立模型在计算机上运行,但这属于数学仿真,随着技术的发展,许多应用领域建立仿真系统时除了模型外还要求将实物和人员包含在仿真回路中。由此可见，基于仿真设计与传统设计的方法和流程两者有很大区别(图2) ,基于仿真设计可以在计算机上建立虚拟样机,对产品的外形、结构、强度、动力

数学应用软件作业2Matlab作图

注意：上机作业文件夹以自己的姓名学号命名，文件夹中包括如下上机报告和Matlab程序。

5、用surf ，mesh 绘制曲面22 2z x y =+，]3,3[],3,3[-∈-∈y x 。 6、用polar 绘制阿基米德螺线r a θ=和三叶玫瑰线cos3r a θ=。（a=100） 7、在同一平面的两个窗口中分别画出心形线和马鞍面。三. 上机方法与步骤给出相应的问题分析及求解方法，并写出Matlab 程序。并有上机程序显示。第1题：要在同一坐标系中分别画出四个函数的图形，可以利用plot 来将这四条曲线画在一起。 Matlab 程序： x=linspace(15,200,50); y1=(1+1./x).^x; y2=(1+1./x).^(x+1); y3=[1+1./(x+1)].^x; y4=2.7183; plot(x,y1,'r',x,y2,'g',x,y3,'b',x,y4,'co') 第2题：直接用ezplot 命令绘制函数的图形。 Matlab 程序： ezplot('（exp(x*y)-sin(x+y)）',[-3,3])

第3题：取a=1，直接用ezplot命令绘出这两条曲线。 Matlab程序： ezplot('[t-sin(t)]','[1-cos(t)]',[0,2*pi]) 第4题：利用subplot命令将这四个函数的图像画在两行两列的同一个图形上。Matlab程序： x=[0:0.1:5]; subplot(2,2,1); plot(x,sin(5*x)); grid on; title('plot-y=sin(5*x)'); subplot(2,2,2); plot(x,cos(3*x));

科技写作结课作业(时域有限差分法的Matlab仿真开题报告)

开题报告论文题目：基于matlab的时域有限差分法的电磁仿真研究（10分）学院：电气工程及其自动化学院学号：1103000105姓名：__杨志刚___ 一、论文选题的目的和意义（300字以内；15分）时域有限差分法，因具有多种优点被运用到电磁场理论研究的各个方面，而且其使用成效和应用领域还在迅速扩大和提高，在现代电磁场理论研究中具有很大的重要性和很强的可操作性。但是同时这种方法也存在一定的缺陷，主要表现在对无边界问题需要吸收边界条件处理，有色散误差，消耗内存大等方面。本课题在利用时域有限差分法对一些实际的算例进行实验仿真和验证，同时对这种方法在解决实际问题的缺陷进行一定程度的研究和分析。 Matlab作为一种工程仿真工具得到了广泛应用。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。二、国内外关于该论题的研究现状和发展趋势（500字以内；15分）时域有限差分方法作为一种典型的全波时域分析方法，因其原理直观、编程简便、实用性强在目前的计算电磁学领域内被人们广泛深入地研究，并取得巨大应用成功的方法。时域数值技术的一个突出优点是可以给出关于问题空间的丰富的时域信息，而且经过简单的时频变换，即可得到宽带范围的频域信息，相对频域方法显著地节约了计算量。最近几十年，是电磁场数值计算时域技术蓬勃发展的时期，各具优势和特色的新颖时域算法层出不穷。但是到目前为止国内关于时域有限差分法中的PML 算法文献较少，其中绝大多数文献集中在综述和应用方面。而在国际的学报和杂志上对于这方面的文献非常多。时域有限差分法经过了三十年多年的高速发展之后，仍然还是计算电磁学制高点的研究热潮，而且其应用的范围和成效还在迅速的扩大和提高。本课题正是利用时域有限差分法的基础理论，利用matlab对一些实际的电磁场问题进行仿真研究。三、论文的主攻方向、主要内容、研究方法及技术路线（1000字左右；40分）通过对时域有限差分法理解基础之上，利用matlab仿真软件按照这种方法编程，实现对三种情况下的电磁场情况的仿真研究。

matlab教程课后作业

【例1.3-5】图示复数i z i z 21,3421+=+=的和。 z1=4+3*i;z2=1+2i; z12=z1+z2 clf,hold on plot([0,z1,z12],'-b','LineWidth',3) plot([0,z12],'-r','LineWidth',3) plot([z1,z12],'ob','MarkerSize',8) hold off,grid on axis equal axis ([0,6,0,6]) text(3.5,2.3,'z1') text(5,4.5,'z2') text(2.5,3.5,'z12') xlabel('real') ylabel('image') shg z12 = a=-8; r_a=a^(1/3) p=[1,0,0,-a]; R=roots(p) MR=abs(R(1)); t=0:pi/20:2*pi; x=MR*sin(t); y=MR*cos(t); plot(x,y,'b:'),grid on hold on plot(R(2),'.','MarkerSize',30,'Color','r') plot(R([1,3]),'o','MarkerSize',15,'Color','b') axis([-3,3,-3,3]),axis square

hold off r_a = 1.0000 + 1.7321i R = -2.0000 1.0000 + 1.7321i 【例1.3-10】画出衰减振荡曲线t e y t 3sin 3-=，t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y=exp(-t/3).*sin(3*t); plot(t,y,'r','LineWidth',2) axis([0,4*pi,-1,1]) xlabel('t'),ylabel('y')

春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)。

东南大学信号与系统试题含答案

东南大学考试卷（A 、B 卷）（答案附后）课程名称信号与线性系统考试学期 03-04-3 得分适用专业四系，十一系考试形式闭卷考试时间长度 120分钟一、简单计算题（每题8分）： 1、已知某连续信号()f t 的傅里叶变换为 21 ()23F j j ωωω= -+，按照取样间隔1T =对其进行取样得到离散时间序列()f k ，序列()f k 的Z 变换。 2、求序列{} 10()1,2,1 k f k ==和2()1cos ()2f k k k πε????=+ ???????的卷积和。 3、已知某双边序列的Z 变换为 21 ()1092F z z z = ++，求该序列的时域表达式()f k 。

4、已知某连续系统的特征多项式为： 269111063)(234567+++++++=s s s s s s s s D 试判断该系统的稳定情况，并指出系统含有负实部、零实部和正实部的根各有几个？ 5、已知某连续时间系统的系统函数为： 323 2642 ()21s s s H s s s s +++=+++。试给出该系统的状态方程。 6、求出下面框图所示离散时间系统的系统函数。 ) (k

二、（12分）已知系统框图如图（a ），输入信号e(t)的时域波形如图（b ），子系统h(t)的冲激响应波形如图(c)所示，信号()f t 的频谱为 ()jn n F j e πω ω+∞ =-∞ = ∑ 。图(a) y(t) ) (t f e(t)图(b) h(t)图(c) 试：1）分别画出)(t f 的频谱图和时域波形； 2）求输出响应y(t)并画出时域波形。 3）子系统h(t)是否是物理可实现的？为什么？请叙述理由；

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

MATLAB结课作业

4.10 上机操作步骤 1在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on n=1:1000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') 观察数列的散点图22,当n 趋于无穷大时,数列趋于 0 subplot(1,2,2) hold on grid on n=500:10000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') fplot('0.001',[500,10000]) fplot('-0.001',[500,10000]) axis([500,10000,-0.005,0.005]) 观察图23,当001.0=ε时,可以取N= 1000 ,当n>N 时有επε<< -2 co s n 1n . 图22 图23 2 在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on fplot('x.*x',[1,3])

观察函数图24, 当2x →时,2x y =的极限是 4 subplot(1,2,2) hold on grid on fplot('x.*x',[1.9,2.1]) fplot('4.001',[ 1.9,2.1]) fplot('3.999',[ 1.9,2.1]) axis([1.9997,2.0005,3.9989,4.0011]) % 调整显示图形的范围是该实验的重点观察图25,当001.0=ε时, δ取 0.003 δ<-<2 0x 时,001.04<-y ? 图24 图25 3 在MatLab 的命令窗口输入: syms x limit((2.^x-log(2.^x)-1)./(1-cos(x)),x,0) 运行结果为 ans = log(2)^2 理论上用洛必达法则计算该极限: x x x cos 112ln 2lim 0x ---→= 1 4 在MatLab 的命令窗口输入如下命令序列: （1）syms x y=sqrt(x+2)*(3-x)^4/(x+1)^5 diff(y,x) ％求一阶导数运行结果＝y'1/2/(x+2)^(1/2)*(3-x)^4/(x+1)^5-4*(x+2)^(1/2)*(3-x)^3/(x+1)^5-5*(x+2)^(1/2)*(3-x )^4/(x+1)^6 x=1; eval(y) ％求导数在x ＝1处的值运行结果 1'=x y = 0.8660

有限元大作业matlab课程设计例子

有限元大作业程序设计学校：天津大学院系：建筑工程与力学学院专业：01级工程力学姓名：刘秀学号：\\\\\\\\\\\ 指导老师：连续体平面问题的有限元程序分析 [题目]：如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]：由于连续平板的对称性，只需要取其在第一象限的四分之一部分

参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。 [ 用和单元信息文件DATA.OUT。位移模式：用用线性位移模式载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷材料性质：弹性体由单一的均匀材料组成约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束方程求解：针对半带宽刚度方程的Gauss消元法输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

（1） ID ： ID=2时为平面应变问题（平面问题），LJK_ELE(I,1)，LJK_ELE(I,2)， X(I),Y(I)分别存放节点I 的x ，y 表示第I 个作用有节点载荷的节点x,y 方向的节点载荷数值存放节点载荷向量，解方程后该矩（2 READ_IN ：读入数据 BAND_K ：形成半带宽的整体刚度矩阵 FORM_KE ：计算单元刚度矩阵 FORM_P ：计算节点载荷 CAL_AREA ：计算单元面积 DO_BC ：处理边界条件 CLA_DD ：计算单元弹性矩阵 SOLVE ：计算节点位移 CLA_BB ：计算单元位移……应变关系矩阵 CAL_STS ：计算单元和节点应力（3）文件管理：源程序文件： chengxu.for 程序需读入的数据文件：

利用Matlab实现Romberg数值积分算法----系统建模与仿真结课作业

利用Matlab 实现Romberg 数值积分算法一、内容摘要针对于某些多项式积分，利用Newton —Leibniz 积分公式求解时有困难，可以采用数值积分的方法，求解指定精度的近似解，本文利用Matlab 中的.m 文件编写了复化梯形公式与Romberg 的数值积分算法的程序，求解多项式的数值积分，比较两者的收敛速度。二、数值积分公式 1.复化梯形公式求解数值积分的基础是将区间一等分时的Newton —Cotes 求积公式: I =(x)[f(a)f(b)]2 b a b a f dx -≈ +? 其几何意义是，利用区间端点的函数值、与端点构成的梯形面积来近似(x)f 在区间[a,b]上的积分值，截断误差为： 3" (b a)()12 f η-- (a,b)η∈ 具有一次的代数精度，很明显，这样的近似求解精度很难满足计算的要求，因而，可以采用将积分区间不停地对分，当区间足够小的时候，利用梯形公式求解每一个小区间的积分近似值，然后将所有的区间加起来，作为被求函数的积分，可以根据计算精度的要求，划分对分的区间个数，得到复化梯形公式： I =1 1 (b a)(b a) (x)dx [f(a)f(b)2(a )]2n b a k k f f n n -=--≈+++∑? 其截断误差为：

2" (b a)h ()12 R f η--= (a,b)η∈ 2.Romberg 数值积分算法使用复化的梯形公式计算的数值积分，其收敛速度比减慢，为此，采用Romberg 数值积分。其思想主要是，根据I 的近似值2n T 加上I 与2n T 的近似误差，作为新的I 的近视，反复迭代，求出满足计算精度的近似解。用2n T 近似I 所产生的误差可用下式进行估算： 12221 ()3 n n n I T T T -?=-=- 新的I 的近似值： 122 n n j T T -=?+ j =（0 1 2 ….） Romberg 数值积分算法计算顺序 i=0 (1) 002T i=1 (2) 102T (3) 012T i=2 (4) 202T (5) 112T (6) 022T i=3 (7) 302T (8) 212T (9) 122T (10) 032T i=4 (11) 402T (12) 312T (13) 222T (14) 132T … … … … 其中，第一列是二阶收敛的，第二列是四阶收敛的，第三列是六阶收敛的，第四列是八阶收敛的，即Romberg 序列。

《科学计算与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)]] %输出矩阵

matlab结课论文

精选文库山西大同大学matlab课程结课作业MATLAB程序应用姓名：课程序号： 2 班级：学号： 2013年12月

1.实验内容：已知!123n n =?????L ，编写一个程序求满足100!10n ≤的最大的n 值以及此时!n 的值。 function n n=2;m=1; while m<=10^100 m=m.*n;n=n+1; end m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131 1(22Λ+--++--+=π，试根据公式编出计算pi 的Mat lab 主程序文件，pi 的精度为0.00001。程序： k=0;n=1;b=0;a=0; while abs((pi-a))>0.00001 a=2*sqrt(2)*k; k=( bcos( *pi/2)+sin(b*pi/2))/n+k; n=n+2; b=b+1; end a 输出a=3.141602572083633 ; a-pi= 9.918493839577991e-006 3.有两个矩阵A 和B 如下：????????????---=771175420132861-1A ，????????????------=0162310013125673B ，将A 中所有等于-1的元素改为-2，将B 中所有小于0的元素改为1，然后将B 中等于0的元素的值改为A 的相应位置元素的值。请用Matlab 函数文件实现上述运算。 clear;

clc; A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0]; C=A;A(A==-1)=-2;U=A; D=B;B(B<0)=1;V=B; A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D; A,B,W,U,V %用函数文件实现矩阵中元素的变换。 %A、B为输入变量。 %U、V、W分别存放A、B中间变换结果。 ; 4.用matlab主程序文件产生动画：呈现一小圆（半径为1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。 clea close;clc;r; axis([-6 6 -6 6],'equal','manual');hold on; ezplot('x^2+y^2-9'); h=ezplot('x^2+y^2-1'); x=get(h,'xdata'); y=get(h,'ydata'); for t=1:7200 set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t)); drawnow; end

东南大学信号与系统本科试卷答案

1、已知某连续信号()f t 的傅里叶变换为 21 ()23F j j ωωω= -+，按照取样间隔1T =对其进行取样得到离散时间序列()f k ，序列()f k 的Z 变换。解法一：f(t)的拉普拉斯变换为 21 11)2)(1(1321)(2+- +=++=++= s s s s s s s F ， 2111)(Re )(--===---=-=??? ???-=∑∑e z z e z z e z z K e z z s F s z F n i T s i s s n i sT i i 解法二：f(t)=L -1{F(jw)}=(e -t - e -2t )ε(t) f(k)= (e -k - e -2k )ε(k)=)())()((21k e e k k ε--- F(z)=Z[f(k)]= 21 -----e z z e z z 2、求序列{} 10()1,2,1 k f k ==和2()1cos ()2f k k k πε?? ??=+ ???????的卷积和。解：f 1(k)={1,2,1}=δ(k)+2δ(k -1)+ δ(k -2) f 1(k)* f 2(k)= f 2(k)+ 2f 2(k -1)+ f 2(k -2) 3、已知某双边序列的Z 变换为 21 ()1092F z z z = ++，求该序列的时域表达式()f k 。解： 5.01 4.01)(+- += z z z F ，两个单阶极点为-0.4、-0.5 当收敛域为|z|>0.5时，f(k)=(( -0.4)k -1-( -0.5)k -1)ε(k -1) 当收敛域为0.4<|z|<0.5时，f(k)= ( -0.4)k -1ε(k -1)+( -0.5)k -1ε( -k) 当收敛域为|z|<0.4时，f(k)= - ( -0.4)k -1ε(-k)+( -0.5)k -1ε( -k) 点评：此题应对收敛域分别讨论，很多学生只写出第一步答案，即只考虑单边序列。 4、已知某连续系统的特征多项式为： 269111063)(234567+++++++=s s s s s s s s D 试判断该系统的稳定情况，并指出系统含有负实部、零实部和正实部的根各有几个？解构作罗斯-霍维茨阵列 611617s 29 1036s 3168 385s 2314s 342 (00) 32s s s ++此时出现全零行，有辅助多项式 3 4646,4,6s s +求导可得以代替全零行系数。 21 32 2232s s s 由罗斯-霍维茨数列可见，元素符号并不改变，说明s 右半平面无极点。再由

MATLAB大作业

M A T L A B大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。第一类：绘制图形。（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三

角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如），请分析构图原理并用MATLAB 实现。问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。第二类：插值与拟合。（B级）问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未减速。在一组时间点上测得汽车的速度如表所示。（1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。（2）绘制插值图形并标注样本点。

(完整版)有限元大作业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 程序的原理如框图：

matlab结课论文

山西大同大学matlab课程结课作业MATLAB程序应用姓名：课程序号： 2 班级：学号： 2013年12月

1.实验内容：已知!123n n =????? ，编写一个程序求满足100!10n ≤的最大的n 值以及此时!n 的值。 function n n=2;m=1; while m<=10^100 m=m.*n;n=n+1; end m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131 1(22 +--++--+=π，试根据公式编出计算pi 的Mat lab 主程序文件，pi 的精度为0.00001。程序： k=0;n=1;b=0;a=0; while abs((pi-a))>0.00001 a=2*sqrt(2)*k; k=( bcos( *pi/2)+sin(b*pi/2))/n+k; n=n+2; b=b+1; end a 输出a=3.141602572083633 ; a-pi= 9.918493839577991e-006 3.有两个矩阵A 和B 如下：????????????---=771175420132861-1A ，????????????------=0162310013125673B ，将A 中所有等于-1的元素改为-2，将B 中所有小于0的元素改为1，然后将B 中等于0的元素的值改为A 的相应位置元素的值。请用Matlab 函数文件实现上述运算。

clear; clc; A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0]; C=A;A(A==-1)=-2;U=A; D=B;B(B<0)=1;V=B; A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D; A,B,W,U,V %用函数文件实现矩阵中元素的变换。 %A、B为输入变量。 %U、V、W分别存放A、B中间变换结果。 ; 4.用matlab主程序文件产生动画：呈现一小圆（半径为1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。 clea close;clc;r; axis([-6 6 -6 6],'equal','manual');hold on; ezplot('x^2+y^2-9'); h=ezplot('x^2+y^2-1'); x=get(h,'xdata'); y=get(h,'ydata'); for t=1:7200 set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t)); drawnow; end

Matlab程序设计(2016大作业)

Matlab程序设计课程大作业题目名称：_________________________________ 班级：_________________________________ 姓名：_________________________________ 学号：_________________________________ 课程教师：温海骏学期：2015-2016学年第2学期完成时间： MATLAB优化应用 §1 线性规划模型一、线性规划问题：问题1：生产计划问题假设某厂计划生产甲、乙两种产品，现库存主要材料有A类3600公斤，B类2000公斤，C类3000公斤。每件甲产品需用材料A类9公斤，B类4公斤，C类3公斤。每件乙产品，需用材料A类4公斤，B类5公斤，C类10公斤。甲单位产品的利润70元，乙单位产品的利润120元。问如何安排生产，才能使该厂所获的利润最大。问题2：投资问题某公司有一批资金用于4个工程项目的投资，其投资各项目时所得的净收益(投入资金百分比)如下表：工程项目收益表工程项目 A B C D 收益(%) 15 10

12 由于某种原因，决定用于项目A的投资不大于其他各项投资之和而用于项目B和C的投资要大于项目D的投资。试确定该公司收益最大的投资分配方案。问题3：运输问题有A、B、C三个食品加工厂，负责供给甲、乙、丙、丁四个市场。三个厂每天生产食品箱数上限如下表：工厂 A B C 生产数 60 40 50 四个市场每天的需求量如下表：市场甲乙丙丁需求量 20 35 33 34 从各厂运到各市场的运输费(元/每箱)由下表给出：收点发点市场甲乙丙丁工厂 A 2 1 3 2 B

Matlab程序设计大作业(终审稿)

M a t l a b程序设计大作业公司内部档案编码：[OPPTR-OPPT28-OPPTL98-OPPNN08]

MATLAB优化应用 §1 线性规划模型一、线性规划问题：问题1：生产计划问题假设某厂计划生产甲、乙两种产品，现库存主要材料有A类3600公斤，B类2000公斤，C类3000公斤。每件甲产品需用材料A类9公斤，B类4公斤，C类3公斤。每件乙产品，需用材料A类4公斤，B类5公斤，C类10公斤。甲单位产品的利润70元，乙单位产品的利润120元。问如何安排生产，才能使该厂所获的利润最大。问题2：投资问题某公司有一批资金用于4个工程项目的投资，其投资各项目时所得的净收益(投入资金百分比)如下表：工程项目收益表由于某种原因，决定用于项目A的投资不大于其他各项投资之和而用于项目B和C的投资要大于项目D的投资。试确定该公司收益最大的投资分配方案。问题3：运输问题

有A 、B 、C 三个食品加工厂，负责供给甲、乙、丙、丁四个市场。三个厂每天生产食品箱数上限如下表：四个市场每天的需求量如下表：从各厂运到各市场的运输费(元/每箱)由下表给出：求在基本满足供需平衡的约束条件下使总运输费用最小。 §2 多目标规划模型多目标规划定义为在一组约束下，多个不同的目标函数进行优化设计。数学模型： 12min ()() ().()0,1,2, ,m j f x f x f x st g x j k ???? ≤= 其中x=(x 1 ,x 2 , … ,x n )为一个n 维向量；f i (x)为目标函数，i=1, 2, … ,m; g j (x)为系统约束, j=1, 2, … ,k 。