4.15 修改4.7

中的程序lsqfit，使它能夠從input1.dat

文件中讀取它的輸入值。文件中的數據是以行組織的，每一行都有一對(x，y)，如下所示:

1.1 2.2

2.2 3.3

...

用例4.7 中的數據檢測你的程序。(提示:我們用load

命令從input1 數組讀數據。然后把input1 的第一列賦值于數組x，把input1

的第二列賦值于數組y)。

解答：

disp('This

program performs a leastsquares fit of an ');

disp('input

data set to a straight line.');

n_points =

load('sample_file.txt')

% Read the

input data

for ii =

1:n_points

temp =

input('Enter [x y] pair: ');

x(ii) =

temp(1);

y(ii) =

temp(2);

end

% Accumulate

statistics

sum_x =

0;

sum_y =

0;

sum_x2 =

0;

sum_xy =

0;

for ii =

1:n_points

sum_x = sum_x

+ x(ii);

sum_y = sum_y

+ y(ii);

sum_x2 =

sum_x2 + x(ii)^2;

sum_xy =

sum_xy + x(ii) * y(ii);

end

% Now

calculate the slope and intercept.

x_bar = sum_x

/ n_points;

y_bar = sum_y

/ n_points;

slope =

(sum_xy - sum_x * y_bar) / ( sum_x2 - sum_x * x_bar);

y_int = y_bar

- slope * x_bar;

% Tell

user.

disp('Regression coefficients

for the leastsquares line:');

fprintf('

Slope (m) = %8.3f\n', slope);

fprintf('

Intercept (b) = %8.3f\n', y_int);

fprintf(' No

of points = \n', n_points);

% Plot the

data points as blue circles with no

% connecting

lines.

plot(x,y,'bo');

hold

on;

xmin =

min(x);

xmax =

max(x);

ymin = slope *

xmin + y_int;

ymax = slope *

xmax + y_int;

plot([xmin

xmax],[ymin ymax],'r','LineWidth',2);

hold

off;

title

('\bfLeastSquaresFit');

xlabel('\bf\itx');

ylabel('\bf\ity');

legend('Input

data','Fitted line');

grid

on

4.17 在例4.3

中，已知年月日，計算相應的thedayofyear。在這個程序中，并沒有檢測是否輸入了正確的年月日，它能接收無效的月和日，并產生無意義的結果。修改你的程序使之只能輸入有效的年月日。如果輸入的值無效，則提示用戶出錯，并且跳出執行。我們要求年應當大于0月只能是1

到12 之間的整數。日只能1 到那一月的最大數之間的整數。用switch 結構檢查日是否正確。

解答：

disp('This

program calculates the day of year given the ');

disp('current

date.');

month = input('Enter

current month (1-12):');

day = input('Enter

current day(1-31):');

year = input('Enter

current year(yyyy): ');

% Check for leap

year, and add extra day if necessary

if mod(year,400) ==

0

leap_day = 1; %

Years divisible by 400 are leap years

elseif mod(year,100)

== 0

leap_day = 0; %

Other centuries are not leap years

elseif mod(year,4)

== 0

leap_day = 1; %

Otherwise every 4th year is a leap year

else

leap_day = 0; %

Other years are not leap years

end

% Calculate day of

year by adding current day to the

% days in

previous months.

day_of_year =

day;

for ii =

1:month - 1

% Add days in

months from January to last month

switch

(ii)

case

{1,3,5,7,8,10,12},

day_of_year =

day_of_year + 31;

case

{4,6,9,11},

day_of_year =

day_of_year + 30;

case

2,

day_of_year =

day_of_year + 28 + leap_day;

end

end

% Tell

user

fprintf('The

date -/-/M is day of year %d.\n', ...

month, day,

year, day_of_year);

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 循環進行計算。

解答：

n=input('請輸入一個數n');

f(1)=1;

f(2)=2;

if(n>2)

for

i=3:n

f(i)=f(i-1)+f(i-2);

end

end

fprintf('f(%d)的值為%d',n,f(n));

4.21

輕繩上的拉力。一重200

英磅的物體被固定在一水平桿的末端。如圖4.5 所示這一水平桿由一輕繩固定。繩子上的拉力為

T 代表繩子的拉力，W 代表物體的重量，lp 代表桿的長度，lc

為繩長，d 代表繩與桿的結點到墻面的距離。編寫一個程序，以確定d 為多大時，繩的拉力最小。為達此目的，d

應從1 英尺到7 英尺，每隔1

英尺取一次值，并找出使拉力最小的d。

4.23 分貝

工程師們經常用分貝或dB 來描述兩功率之比。1dB

的定義如下，

P2 是已測量的功率，P1 代表參考功率。假設參考功率P1 是1

瓦。P2 從1 到20 瓦每隔0.5

瓦取一次值，編寫程序，計算相應的dB 值，并畫出dB-P2

解答：

P1=1;

P2=1:0.5:20;

dB=10*log10(P2);

plot(P2,dB);

4.25

均方根平均數(rmsaverage)。均方根平均數是另一種計算數據平均數的方法。它的定義如下

編寫一個程序，它能接受任意個數的正輸入值，并計算它們的算太平均數和幾何平均數。用while

循環讀取輸入值，當輸入一個負數中止輸入數據。計算數列10，5，2，5 的均方根平均數，用以檢測程序。

解答：

x=input('請輸入一個數：(以負數結束輸入)');

sum=0;

n=0;

m=1;

ad=0;

while(x>=0)

sum=sum+x;

m=m*x;

ad=ad+1/x;

x=input('請輸入下一個數：');

n=n+1;

end

hm=n/ad;

ss=sum/n;

mm=m^(1/n);

fprintf('算術平均數為%f\n幾何平均數為%f\n調和平均數為',ss,mm,hm);

4.27

編寫一個程序，能夠計算一系列正數的算術平均數，幾何平均數，均方根平均數，調和均數。可使用任意算法讀取輸入值。用下列數測試你的程序。

總結

以上是生活随笔為你收集整理的matlab作业参考4,matlab第四章作业的全部內容，希望文章能夠幫你解決所遇到的問題。

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