没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB程序设计与应用 刘卫国主编
需积分: 14 11 下载量 175 浏览量
2010-06-19
19:09:25
上传
评论
收藏 111KB PDF 举报
温馨提示
试读
21页
MATLAB程序设计与应用 刘卫国主编
资源推荐
资源详情
资源评论
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验3.txt
SY301
clear all
m=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];
for x=m
if x<0&x~=-3
y1=x^2+x-6;
disp(['y= ',num2str(y1)])
elseif x>=0&x<5&x~=2&x~=3
y2=x^2-5*x+6;
disp(['y= ',num2str(y2)])
else
y3=x^2-x-1;
disp(['y= ',num2str(y3)])
end
end
SY302
clear all
x=input('请输入成绩:');
while (x>100|x<0)
x=input('输入错误,请重新输入成绩:');
end
switch fix(x/10)
case {9,10}
disp('A')
case {8}
disp('B')
case {7}
disp('C')
case {6}
disp('D')
case {0,1,2,3,4,5}
disp('E')
end
SY303
clear all
x=input('请输入工号:');
y=input('请输入工作时长:');
if y>120
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验3.txt(第 1/2 页)2009-5-10 14:32:04
课后答案网 www.khdaw.com
www.khdaw.com
课后答案网
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验3.txt
a=84*120+(y-120)*84*0.15;
disp(['您本月工资是: ',num2str(a)])
elseif y<60
b=84*y-700;
disp(['您本月工资是: ',num2str(b)])
else
c=84*y;
disp(['您本月工资是: ',num2str(c)])
end
SY304
clear all
a=fix(rand(1)*89)+10;
disp(['a=',num2str(a)])
b=fix(rand(1)*89)+10;
disp(['b=',num2str(b)])
x=input('请输入一个四则运算符号','s');
if abs(x)==43
c=a+b;
elseif abs(x)==45
c=a-b;
elseif abs(x)==42
c=a*b;
elseif abs(x)==47
c=a/b;
end
disp(['a于b的计算结果是:',num2str(c)])
SY305
clear all
a=fix(rand(5,6)*89)+10;
n=input('请输入一个数(输出结果是一个5*6的矩阵的该行元素):');
if n>5|n<0
b=a(5:5,:);
disp(['输入错误,程序输出最后一行的元素:',num2str(b)])
else
c=a(n:n,:);
disp(['输出结果是:',num2str(c)])
end
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验3.txt(第 2/2 页)2009-5-10 14:32:04
课后答案网 www.khdaw.com
www.khdaw.com
课后答案网
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt
SY401
clear all
s1=0;
for n=1:100
x=1/(n^2);
s1=s1+x;
end
disp(['当N等于100时: ',num2str(s1)])
s2=0;
for n=1:1000
x=1/(n^2);
s2=s2+x;
end
disp(['当N等于1000时: ',num2str(s2)])
s3=0;
for n=1:10000
x=1/(n^2);
s3=s3+x;
end
disp(['当N等于10000时: ',num2str(s3)])
SY402
clear all
n=1;
y=0;
while (y<3)
x=1/(2*n-1);
n=1+n;
y=y+x;
end
disp(['y<3时n的最大值是:',num2str(n-2)])
disp(['相应的的y值是: ',num2str(y-x)])
SY403
clear all
x=input('请输入x的初始值:');
a=input('请输入常数a的值:');
b=input('请输入常数b的值:');
g=a/(b+x);
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt(第 1/4 页)2009-5-10 14:32:05
课后答案网 www.khdaw.com
www.khdaw.com
课后答案网
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt
n=0;
while ((abs(x-g))>10^-5)
g=a/(b+x);
x=a/(b+g);
n=n+1;
if n>499
break
end
end
disp(['迭代结果是:',num2str(x)])
disp(['循环次数是',num2str(n),'次'])
r1=(-b+sqrt(b^2+4*a))/2;
r2=(-b-sqrt(b^2+4*a))/2;
disp(['如果迭代收敛,准确值是:',num2str(r1),' 或者是:',num2str(r2)])
SY404
clear all
for n=1:4
if n==1
f1=1;
elseif n==2
f2=0;
elseif n==3
f3=1;
else
a=f3-2*f2+f1;
b=a-2*f3+f2;
c=b-2*a+f3;
d=c-2*b+a;
H=[1,0,1,a,b,c,d];
for m=8:4:99
a=d-2*c+b;
b=a-2*d+c;
c=b-2*a+d;
d=c-2*b+a;
H=[H,a,b,c,d];
end
f100=d-2*c+b;
end
end
max=max(H);
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt(第 2/4 页)2009-5-10 14:32:05
课后答案网 www.khdaw.com
www.khdaw.com
课后答案网
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt
min=min(H);
sum=sum(H);
disp(['最大值是:',num2str(max)])
disp(['最小值是:',num2str(min)])
disp(['各数和是:',num2str(sum)])
k=0;
l=0;
p=0;
for e=H
if e>0
k=k+1;
elseif e<0
l=l+1;
else
p=p+1;
end
end
disp(['正数的个数是:',num2str(k)])
disp(['负数的个数是:',num2str(l)])
disp(['零的个数是: ',num2str(p)])
SY405
clear all
p=0;
l=0;
H=[];
for m=2:49
a=m;
b=m+1;
c=a*b-1;
for k=1:c
if rem(c,k)==0
l=l+1;
end
end
if l==2
disp(['亲密数对是','(',num2str(a),' ,',num2str(b),')'])
p=p+1;
H=[H,c];
end
l=0;
file:///F|/khdaw/下载答案/5月10日/matlab/matlab/matlab/实验4.txt(第 3/4 页)2009-5-10 14:32:05
课后答案网 www.khdaw.com
www.khdaw.com
课后答案网
剩余20页未读,继续阅读
资源评论
litao11046
- 粉丝: 7
- 资源: 3
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功