matlab第四次作业

实验四 3：

a=input('请输入a的值'); b=input('请输入b的值'); X0=1.0; X1=a/(b+X0); n=1;

while (abs(X1-X0)<=(10^-5))&&(n<=500) X1=a/(b+X0); X0=X1; n=n+1; end disp(X0);

r1=(-b+sqrt(b*b+4*a))/2; r2=(-b-sqrt(b*b+4*a))/2; disp(r1); disp(r2); 结果： 请输入a的值1 请输入b的值1 1

0.6180

-1.6180

请输入a的值8 请输入b的值3 1

1.7016

-4.7016

请输入a的值10 请输入b的值0.1 1

3.1127

-3.2127

%第五题：

n=0; s=0; fori=2:49 a=i*(i+1)-1; m=fix(sqrt(a)); for j=2:m if rem(a,j)==0 break end end if j==m s=s+a; n=n+1; end end

结果： n n = 28 >>s s =

21066

实验五

%第二题程序一

function [a,b,N,M]=shiyanwu2(m,n,t)

A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1]; B=[0,9.8*m,0,9.8*n]; C=inv(A)*B'; a=C(1); b=C(2); N=C(3); M=C(4);

%在命令窗口调用该函数文件: m1=input('m1='); m2=input('m2=');

theta=input('theta=');

[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)

%第二题程序二

function X=qiujie1(m1,m2,t) g=9.8;

A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0;m1*sin(t*pi/180),0,cos(t

*pi/180),0;0,m2,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1]; B=[0;m1*g;0;m2*g]; X=inv(A)*B;

结果：

>> a=qiujie1(2,3,30) a =

7.0000 2.4249 14.5492 42.0000

%第五题 %(1)

function f1=mat5(n) f1=n+10*log(n*n+5);

结果：

>> y=mat5(40)/(mat5(30)+mat5(20)) y =

0.6390

%(2)方法二

function f2=mat52(i) n=1:i;

m=n.*(n+1); f2=sum(m); end

结果：

>> b=mat52(3) b =

20

本文来源：https://www.bwwdw.com/article/z1vf.html