4.1

for

t=-9:0.5:9

if(t>=0)

y=-(3*t^2)+5;

fprintf('y=%.2ft=%.2f\n',y,t);

else

y=(3*t^2)+5;

fprintf('y=%.2ft=%.2f\n',y,t);

end

end

编译结果：

y=248.00t=-9.00

y=221.75t=-8.50

y=197.00t=-8.00

y=173.75t=-7.50

y=152.00t=-7.00

y=131.75t=-6.50

y=113.00t=-6.00

y=95.75t=-5.50

y=80.00t=-5.00

y=65.75t=-4.50

y=53.00t=-4.00

y=41.75t=-3.50

y=32.00t=-3.00

y=23.75t=-2.50

y=17.00t=-2.00

y=11.75t=-1.50

y=8.00t=-1.00

y=5.75t=-0.50

y=5.00t=0.00

y=4.25t=0.50

y=2.00t=1.00

y=-1.75t=1.50

y=-7.00t=2.00

y=-13.75t=2.50

y=-22.00t=3.00

y=-31.75t=3.50

y=-43.00t=4.00

y=-55.75t=4.50

y=-70.00t=5.00

y=-85.75t=5.50

y=-103.00t=6.00

y=-121.75t=6.50

y=-142.00t=7.00

y=-163.75t=7.50

y=-187.00t=8.00

y=-211.75t=8.50

y=-238.00t=9.00

4.3

for

i=1:50

ii=i^2;

fprintf('i=%dii=%d\n',i,ii);

add=add+i+ii;

end

fprintf('和为%d',add);

编译结果：

i=2ii=4

i=3ii=9

i=4ii=16

i=5ii=25

i=6ii=36

i=7ii=49

i=8ii=64

i=9ii=81

i=10ii=100

i=11ii=121

i=12ii=144

i=13ii=169

i=14ii=196

i=15ii=225

i=16ii=256

i=17ii=289

i=18ii=324

i=19ii=361

i=20ii=400

i=21ii=441

i=22ii=484

i=23ii=529

i=24ii=576

i=25ii=625

i=26ii=676

i=27ii=729

i=28ii=784

i=29ii=841

i=30ii=900

i=31ii=961

i=32ii=1024

i=33ii=1089

i=34ii=1156

i=35ii=1225

i=36ii=1296

i=37ii=1369

i=38ii=1444

i=39ii=1521

i=40ii=1600

i=41ii=1681

i=42ii=1764

i=43ii=1849

i=44ii=1936

i=45ii=2025

i=46ii=2116

i=47ii=2209

i=48ii=2304

i=49ii=2401

i=50ii=2500

和为46750>>

4.5

n=input('请输入N的值');

if(n>0)

x=n*(n-1);

fprintf('N!=%d',x);

elseif(n==0);

disp('N！=1');

else

disp('error');

end

编译结果：

请输入N的值0

N！=1

请输入N的值-8

error

请输入N的值3

N!=6

4.7

a.ires=21循环21次；

b.ires=28循环4次；

c.ires=28循环4次；

d.ires=24循环11次；

4.9

a.

arr1

=

1-23-4

5-67-8

9-1011-12

b.

arr1

=

0000

0364964

81100121144

4.11

假设一个小球的初始位置为(x0，y0)为(0，0)，初速度为20m/s，水平速度和合速度之间

的夹角为θ

度，编写一个程序，画出这个小球的轨迹，并计算小球再次落地与初始位置之间的距离。这程序应当能画出多个抛物线，θ 的取值从5 到85

度，每隔10 度取一次，计算水平位移，θ 的取值从0 到90 度，每隔1 度取一次。最终应当确定那一个θ

值使得水平位移最大。还有打印不同的抛物线时要用不同的颜色。

conv =

pi / 180;

g =

-9.81;

vo =

20;

range

= zeros(1,91);

for ii

= 1:91

theta

= ii - 1;

vxo =

vo*

cos(theta*conv);

vyo =

vo * sin(theta*conv);

traj_time

= -2 * vyo / g;

range(ii)

= vxo * traj_time;

end

fprintf

('Range versus angle theta:\n');

for ii

= 1:91

theta

= ii - 1;

fprintf('

- %8.4f\n',theta, range(ii));

end

[maxrange

index] = max(range);

maxangle

= index - 1;

fprintf

('\nMax range is %8.4f at - degrees.\n',maxrange,

maxangle);

for ii

=[5 15 25 35 45 55 65 75 85];

theta

= ii;

vxo =

vo * cos(theta*conv);

vyo =

vo * sin(theta*conv);

traj_time

= -2 * vyo / g;

x =

zeros(1,21);

y =

zeros(1,21);

for jj

= 1:21

time =

(jj - 1) * traj_time/20;

x(jj)

= vxo * time;

y(jj)

= vyo * time + 0.5 * g * time^2;

end

plot(x,y,'b');

if ii

== 5

hold

on;

end

end

title

('\bfTrajectory of Ball vs Initial Angle \theta');

xlabel

('\bf\itx \rm\bf(meters)');

ylabel

('\bf\ity \rm\bf(meters)');

axis

([0 45 0 25]);

grid

on;

vxo =

vo * cos(maxangle*conv);

vyo =

vo * sin(maxangle*conv);

traj_time

= -2 * vyo / g;

x =

zeros(1,21);

y =

zeros(1,21);

for jj

= 1:21

time =

(jj - 1) * traj_time/20;

x(jj)

= vxo * time;

y(jj)

= vyo * time + 0.5 * g * time^2;

end

plot(x,y,'r','LineWidth',3.0);

hold

off

4.13

conv = pi / 180;

g = -9.81;

vo = input(':');

range = zeros(1,91);

forii

= 1:91

theta = ii - 1;

vxo = vo * cos(theta*conv);

vyo = vo * sin(theta*conv);

traj_time = -2 * vyo / g;

range(ii) = vxo * traj_time;

end

fprintf ('Range

versus angle theta:\n');

forii

= 1:91

theta = ii - 1;

fprintf('

- %8.4f\n',theta,

range(ii));

end;

[maxrange index] = max(range);

maxangle = index - 1;

fprintf ('\nMax

range is %8.4f at - degrees.\n',maxrange,

maxangle);

forii

= 5:10:85

theta = ii;

vxo = vo * cos(theta*conv);

vyo = vo * sin(theta*conv);

traj_time = -2 * vyo / g;

x = zeros(1,21);

y = zeros(1,21);

forjj

= 1:21

time =(jj - 1) * traj_time/20;

x(jj) = vxo * time;y(jj) = vyo * time + 0.5 * g *

time^2;

end

plot(x,y,'b');

ifii

== 5

holdon;

end

end

title ('\bfTrajectory

of Ball vs Initial Angle \theta');

xlabel ('\bf\itx

\rm\bf(meters)');

ylabel ('\bf\ity

\rm\bf(meters)');

axis ([0 11 0 6]);

gridon;

vxo = vo * cos(maxangle*conv);

vyo = vo * sin(maxangle*conv);

traj_time = -2 * vyo / g;

x = zeros(1,21);

y = zeros(1,21);

forjj

= 1:21

time = (jj - 1) * traj_time/20;

x(jj) = vxo * time;

y(jj) = vyo * time + 0.5 * g * time^2;

end

plot(x,y,'r','LineWidth',3.0);

holdoff

输入为10m/s时：

输入为20m/s时：

输入为30m/s时：

4.19

斐波那契数列。含有n个数的斐波那契数列的定义如下：

f(1)

=1

f(2)

=2

f(n) = f(n-1)

+ f(n-2)

所以f(3)=f(2)+f(1)=2+1=3，还有更多的数。在M文件中编写一程序，计算并写斐波那

契数列中第n(n>2)个数的值，n由用户输入。用for循环进行计算。

4.21

4.23

工程师们经常用分贝或dB来描述两功率之比。1dB的定义如下，

P2是已测量的功率，P1代表参考功率。假设参考功率P1是1瓦。P2从1到20瓦每隔

0.5瓦取一次值，编写程序，计算相应的dB值，并画出dB-P2

4.25

4.27