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

創作挑戰賽新人創作獎勵來咯，堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的matlab4.0,matlab 4.0的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。