matlab第二次

Matlab素质拓展课第二次作业

第二次课知识点if 选择语句for while 循环语句子函数（多个输入参数输出参数）

递归二分法秦九韶算法等

3.3 写一个名为nexthour的函数，接收一个整数参数，该参数是一天内的某个小时，然后返回下一个小时。假定是12小时制，所以例子中12的下一个小时是1.

% 3.3 写一个名为nexthour的函数，

% 接收一个整数参数，该参数是一天内的某个小时，然后返回下一个小时。

% 假定是12小时制，所以例子中12的下一个小时是1.

function y=nexthour(n)

if n>=1&&n<=11

y=n+1;

elseif n==12

y=1;

else

y=-1;

disp('输入错误！');

end

end

3.6 写一个脚本，提示用户输入分数的分子和分母。如果分母是0，则打印错误信息，指出0不能做分母。如果分母不是0，则打印分数的结果。

% 3.6 写一个脚本，提示用户输入分数的分子和分母。

% 如果分母是0，则打印错误信息，指出0不能做分母。

% 如果分母不是0，则打印分数的结果。

clc;

clear;

fenzi=input('请输入分子:');

fenmu=input('\n请输入分母:');

if fenmu==0

disp('打印错误！分母不能为0！')

else

fprintf('分数=%d/%d\n',fenzi,fenmu);

end

3.11 写一个函数createMToN，它能够创建并返回一个从m到n的整数组成的向量（m是第一个输入参数，n是第二个），不管m是小于n还是大于n。如果m等于n，该向量正好是1*1的，或者说是标量。

%dy_chapter3_11.m cerateMToN

function vec=dy_chapter3_11(m,n)

if m>n

vec=m:-1:n;

elseif m

vec=m:n;

else

vec=m;

end

3.17 写一个脚本，它能生产一个随机整数，并且打印该随机整数是一个偶数还是一个奇数。%dy_chapter3_17.m

clc;

clear;

A=randint(1,1,[0,10000]);

if mod(A,2)==0

fprintf('%d 是偶数\n',A);

else

fprintf('%d 是奇数\n',A);

end

3.18.写一个名为isdivby4的函数，接收整数作为输入参数，如果输入参数能够被4整除，则返回逻辑1，如果不能整除，则返回逻辑假。

function R=dy_chapter3_18(n)

%写一个名称为mydiv4的函数，接收整数作为输入参数，

%如果输入整数能被4整除，则返回1，如果不能整除，则返回0

if mod(n,4)==0

R=1;

else

R=0;

end

end

3.20 一个毕达哥拉斯（勾股）三角形是一组正整数(a,b,c)，它满足a2=b2+c2.写一个名为ispythag的函数接收三个正整数(按照a,b,c的顺序)，如果这三个数能形成该三角形则返回1，不能则返回0.

function y=ispythag(a,b,c)

if a^2==b^2+c^2

y=1;

else

y=0;

end

end

3.26 云层一般分为高、中和低三层。云的高度是决定性因素，但是划分范围还取决于温度。例如，在热带地区的分类可以取决于下面的高度范围（单位英尺）

低：0—6500

中6500—20000

高>20000

%dy_chapter3_26.m

clc;

clear;

A=input('请输入以英尺为单位的云层的高度:');

if A>0&A<6500

fprintf('低\n');

elseif A>6500&A<20000

fprintf('中\n');

else

fprintf('高\n');

end

4.1 写一个for循环，打印出从1.5到3.1的实数列，步长为0.2.

clc;

clear;

for i=1.5:0.2:3.1

fprintf('%f ',i);

end

4.2 写一个函数sumsteps2,计算并返回1到n(每次递增2)的和，n作为参数传递给函数。例如，如果把11传递给了函数，就返回1+3+5+7+9+11的值。

%R=dy_chapter4_2.m

function R=dy_chapter4_2(n)

% 写一个函数dysum2，计算并返回从1到n(每次递增2)的和，n作为参数传递给函数.

% 例如，如果把9传递给了函数，就返回1+3+5+7+9的值.

R=0;

for i=1:2:n

R=R+i;

end

end

4.3 写一个函数prodby2,接收一个正整数n，计算并返回从1到n的奇数的乘积(如果n是偶数，即从1到n-1).

function y=prodby2(n)

y=1:

for i=1:2:n

y=y*i;

end

4.4 提示用户输入一个整数n，并打印出”I love this stuff!”n次。

clc;

clear;

n=input('请输入一个整数：');

for i=1:n

fprintf('I love this stuff!\n');

end

4.14写一个函数，输入参数为矩阵，计算并返回矩阵中所有数据的平均值。

function y=dy_chapter4_14(A)

[m,n]=size(A);

y=sum(sum(A))/(m*n);

end

4.16 创建一个3*5的矩阵，使用内置函数和循环(如果必要的话可以使用if语句)两种方法来完成下面的操作：

A找出每列中的最大值；B 找出每行中的最大值；C找出整个矩阵的最大值。

%dy_chapter4_16.m

%dy_chapter4_16.m

clc;

clear;

A=rand(3,5);

%使用内置函数

colA=max(A) %每列中的最大值

rowA=max(A') %每行中的最大值

maxA=max(max(A)) %整个矩阵的最大值

%使用循环

colA2=A(:,1);

rowA2=A(1,:);

maxA2=A(1,1);

for j=1:5

for i=1:3

if A(i,j)>colA2(i)

colA2(i)=A(i,j); %每列中的最大值

end

end

end

for i=1:3

for j=1:5

if A(i,j)>rowA2(j)

colA2(j)=A(i,j); %每行中的最大值

end

end

end

for i=1:3

for j=1:5

if A(i,j)>maxA2

maxA2=A(i,j); %矩阵中的最大值

end

end

end

4.19 写一个beautyofmath脚本，产生下面的输出。脚本应该从第1行到第9行重复按左对齐方式生成如下表达式，执行指定的操作来获得等式右边显示的结果，严格按下面显示格式打印。

4.20 写一个脚本，打印出下面乘法表：

1

2 4

3 6 9

4 8 12 16

5 10 15 20 25

%dy_chapter4_20.m

clc;

clear;

for i=1:5

for j=1:i

fprintf('%2d ',i*j);

end

fprintf('\n');

end

4.29 写一个脚本，连续提示用户输入正数，并把它们存储到一个向量变量中，直到用户输入一个负数为止。

0.17 给出20个数，1,3,7,13,…,其规律是：第一个数是1，第二个数比第一个数大2，第三个数比第二个数大4….，请写出相应的matlab源码。

0.18 中国网通规定：拨打市内电话时，如果不超过3分钟，则收取话费0.22元；如果通话时间超过3分钟，则超出部分按每分钟0.1元收取通话费，不足一分钟按一分钟计算。设通话时间为t(分钟)，通话费用y(元)，设计一个程序，计算通话的费用.请写出相应的matlab 源码。

0.19 写出已知函数，输入x的值，求y的值。请写出相应的matlab源码。

0.20. 编写一个程序，求1！+2！+…+10!的值.

0.21. 求满足1+2+3+…+n>500的最小的自然数n.

0.22. 求的值。

0.23. 已知函数，请编写程序，输入自变量x的值，输出其相应的函数值。

0.24. 已知S=5+9+13+…+101,分别用for和while写出相应的源码。

0.25. 某市固定电话(市话)的收费标准为：3分钟之内(包括3分钟)收取0.2元;超过3分钟，每分钟(不足1分钟，按1分钟计算)按0.1元收费。设计一个算法，根据通话时间计算话费。

16. 用秦九韶算法求多项式f(x)=2x4+3x3+x2+5x-4，当x=2是的函数值。

17.已知编写一个程序，对每输入一个x值，都得到相应的函数值。

18. 编写一个程序，计算n个数的和：。

19. 利用，编写程序求π的近似值。

20. 求10~999中的水仙花数，所谓水仙花数是一个三位数，它的各位数字的立方和等于该数，例如153是一个水仙花数，因为153=13+53+33.试编写程序，找出所有的水仙花数。clc; %清空command window中的内容

clear;%清空内存中的临时变量

%水仙花数是指一个n 位数( n≥3 )，

%它的每个位上的数字的n 次幂之和等于它本身。

%（例如：1^3 + 5^3 + 3^3 = 153）

%三位水仙花共有4个，153,370,371,407

%方法一，用字符串有利于位数较多的情况

% % n = 3;

% % p = (0:9) .^ n; %预先算好0到9的n次方% % tempStr = '0' * ones(1, n);

% %

% % for x = (10^(n-1) : (10^n - 1))

% % digits = num2str(x, '%d') - tempStr;

% % s = sum(p(digits + 1));

% % if s == x

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

% % end

% % end

%方法二，用子函数

clc;

clear;

%调用my42(i)子函数

for i=100:999

y=my42(i);

if y~=0

disp(y)

end

end

%判断某个三位数是否为水仙花数

function y=my42(x)

y=0;

x1=fix(x/100);% 百位数字

x2=rem(fix(x/10),10); %十位数字

x3=rem(x,10); %个位数字

if x==x1^3+x2^3+x3^3

y=x;

end

end

21.设计算法，交换两个变量a,b的值，并输出交换前后的值。

22.输入三个正实数x,y,z，以它们为三边是否可以构成直角三角形？编写程序输出结果。

23. 写出按从小到大的顺序重新排列x,y,z三个数值的算法。

24. 判断直线ax+by+c=0与圆(x-x0)2+(y-y0)2=r2的位置关系。

25. 输入一个x的值，输出f(x)=2x2+3的一个函数值。

26. 有一函数，写一程序,输入自变量的值,输出函数值。

27. 产生20个两位随机整数，输出其中小于平均值的偶数。

28. 输入20个数，求其中最大数和最小数。要求分别用循环结构和调用MATLAB的max 函数、min函数来实现。

29.编写一个函数文件，求小于任意自然数n的Fibnacci数列各项。Fibnacci数列定义如下：

30.找出m～n之间的全体素数。

function u=fq(m,n)

％找出m到n的全体素数fq(m,n)

％若仅有一个参数则m= 2，给出参数作为n

％若没有参数则m = 2，n = 100

％1999.3

if nargin == 0,

m = 2; n = 100;

elseif nargin == 1,

n = m; m = 2;

end

i=1;

for t = m:n,

if t ~= 2;

for r = 2:t–1,

if fix(t/r)＊r == t/r＊r,

break

end

end

if r==t–1,u(i)=t;i=i+1;end

else

u(i)=t;i=i+1;

end

end

31.输出全部三位的水仙花数。

clc; %清空command window中的内容

clear;%清空内存中的临时变量

%水仙花数是指一个n 位数( n≥3 )，

%它的每个位上的数字的n 次幂之和等于它本身。%（例如：1^3 + 5^3 + 3^3 = 153）

%三位水仙花共有4个，153,370,371,407

%方法一，用字符串有利于位数较多的情况

% % n = 3;

% % p = (0:9) .^ n; %预先算好0到9的n次方

% % tempStr = '0' * ones(1, n);

% %

% % for x = (10^(n-1) : (10^n - 1))

% % digits = num2str(x, '%d') - tempStr;

% % s = sum(p(digits + 1));

% % if s == x

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

% % end

% % end

%方法二，用子函数

clc;

clear;

%调用my42(i)子函数

for i=100:999

y=my42(i);

if y~=0

disp(y)

end

end

%判断某个三位数是否为水仙花数

function y=my42(x)

y=0;

x1=fix(x/100);% 百位数字

x2=rem(fix(x/10),10); %十位数字

x3=rem(x,10); %个位数字

if x==x1^3+x2^3+x3^3

y=x;

end

end

32.编程求出Fibonacci数列中第一个大于1000的元素及序号。%主函数

clc;

clear;

[Fib index]=my43()

%子函数

function [Fib index]=my43()

a=1;

b=1;

Fib=a+b;

index=3;

while Fib<=1000

a=b;

b=Fib;

Fib=a+b;

index=index+1;

end

end

33.求[100,1000]范围内的全部素数之和。%主函数

clc;

clear;

x=150;

result=[];

for i=100:1000

y=my44(i);

if y==1

result=[result i];

end

end

result

%子函数

function result=my44(x)

result=1;

for i=2:x/2

if mod(x,i)==0

result=0;

break;

end

end

end

34.求1—100的平方、平方根和立方根之和。%主函数

clc;

clear;

for x=1:100

[my1 my2 my3]=my45(x)

end

%子函数

function [my1 my2 my3]=my45(x)

my1=x*x; %平方

my2=sqrt(x); %平方根

my3=x^(1/3); %立方根

end

35.求n!, n=1,2, (10)

%主函数

clc;

clear;

for i=1:10

y=my46(i)

end

%子函数，递归

function y=my46(n)

if n==1

y=1;

else

y=n*my46(n-1);

end

end

36.建一个命令文件将变量a,b(包括矩阵)的值互换。%主函数

clc;

clear;

A=[1 2];

B=[3 4;5 6];

[M N]=my47(A,B)

%子函数,交换A，B

function [M N]=my47(A,B)

M=B;

N=A;

end

37.求一元二次方程的根。

%主函数

clc;

clear;

a=11;b=2;c=1;

result=my48(a,b,c)

%子函数

function result=my48(a,b,c)

%a,b,c为系数

if b*b-4*a*c<0

result=NaN;

disp('没有实数解');

elseif b*b-4*a*c==0

result=-b/(2*a);

else

delta=sqrt(b*b-4*a*c);

x1=(-b+delta)/(2*a);

x2=(-b-delta)/(2*a);

result=[x1 x2];

end

38.输入矩形的两边，求面积。

%主函数

clc;

clear;

M=input('请输入两个数,例如[1 2]:')

result=my49(M)

%子函数,求面积

function result=my49(M)

result=M(1)*M(2);

end

39.输入10个数，求其中最大数和最小数。要求分别用子函数和调用max命令。%主函数

clc;

clear;

M=input('请输入十个数，例[1 2 ... 10]:') [mymax mymin]=my50(M)

mymax2=max(M)

mymin2=min(M)

%子函数

function [mymax mymin]=my50(M) [m,n]=size(M);

mymax=M(1);

mymin=M(1);

for i=2:n

if M(i)>mymax

mymax=M(i);

end

if M(i)

mymin=M(i);

end

end

end

40.多项式相加，p1=[1 5 4], p2=[ 1 6 11 6]; %主函数

clc;

clear;

p1=[1 2 1 5 4]

p2=[ 1 6 11 6]

p=my51(p1,p2)

%子函数

function p=my51(p1,p2)

[m1,n1]=size(p1);

[m2,n2]=size(p2);

if n1-n2>0

temp=zeros(1,n1-n2);

p2=[temp p2];

elseif n1-n2<0

temp=zeros(1,n2-n1);

p1=[temp p1];

end

p=p1+p2;

end

41.鸡兔同笼，头有36，脚共100，求鸡兔各多少？%主函数

clc;

clear;

head=36;

foot=100;

[chicken rabbit]=my52(head,foot)

%子函数

function [chicken rabbit]=my52(head,foot)

chicken=-1;

rabbit=-1;

for i=0:head

for j=0:head

if i+j==head && 2*i+4*j==foot

chicken=i;

rabbit=j;

break;

end

end

end

end

42.分别用if语句和switch语句实现以下计算，其中、、的值从键盘输入。

43.产生20个两位随机整数，输出其中小于平均值的偶数。

44. 设，数列是否收敛？若收敛，其值为多少？编程计算其收敛的近似值，精确到6位有效数字。

45.请至少用三种方法编程求。

46.计算1---100之间所有偶数之和。

47.创建Hilbert矩阵，其中a(i,j)=1/(i+j-1)。

%方法一

clc; %清除command window中的内容

clear;%清除Matlab工作空间中保存的变量

M=3;

N=4;

for i=1:M

for j=1:N

a(i,j)=1/(i+j-1); %Hilbert矩阵近似解

end

end

a

%方法二

clc; %清除command window中的内容

clear;%清除Matlab工作空间中保存的变量

M=3;

N=4;

%format %恢复默认数据格式

format rat %用分数形式表示

for i=1:M

for j=1:N

a(i,j)=1/(i+j-1); %Hilbert矩阵

end

end

a

48.任给一组输入x，求函数y=|x|的值。

%

clc; %清除command window中的内容

clear;%清除Matlab工作空间中保存的变量

x=-1

y=abs(x) %对单个变量求绝对值

t=[-2,-1,0,1,2]

y=abs(t) %对一维矩阵求绝对值

m=[-1 -2;

1 2]

y=abs(m) %对多维矩阵求绝对值

49.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。%本程序实现从键盘输入若干数据，遇0结束，然后求其平均值!

clc; %清除command window中的内容

clear;%清除Matlab工作空间中保存的变量

m=input('请输入数据，输入0结束:');

M=[];

while m~=0

M=[M m];

m=input('请输入数据，输入0结束:');

end

M %取值结束后的矩阵B=sum(M) %求和

[row col]=size(M);

avg_B=B/col %求平均

50.求[108,288]之间第一个能被13整除的整数；

求第一个同时能被11、13整除的整数。

clc; %清除command window中的内容

clear;%清除Matlab工作空间中保存的变量

m=108;

n=288;

for i=m:n

if mod(i,11)==0 &&mod(i,13)==0

num=i;

break;

end

end

num