matlab第三章答案

当n等于10000时：1.6448
程序如下：
clear all
n=1:100;
b=1./(4.*n-3)-1./(4.*n-1);
c=sum(b);
disp([当n等于100时：,num2str(c)])
n=1:1000;
b=1./(4.*n-3)-1./(4.*n-1);
c=sum(b);
disp([当n等于1000时：,num2str(c)])
n=1:1000;
b=((2*n).*(2*n))./((2*n-1).*(2*n+1));
c=prod(b);
disp([当n等于1000时：,num2str(c)])
n=1:10000;
b=((2*n).*(2*n))./((2*n-1).*(2*n+1));
c=prod(b);
disp([当n等于10000时：,num2str(c)])
0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937-0.6976 0.7103 d =
-0.016600 0 1.48010 00 2.5365
plot(double(s.x),double(s.y),*);
截图：
p366第4题绘制极坐标曲线影响。，并分析对曲线形状的
function [ output_args ] = untitled2( input_args ) %untitled2 summary of this function goes here %detailed explanation goes here
for i=2:20
if ma(1,i)
m=a(1,i)
end
end
min函数
function n=min(a[])
m=a(1,20);
for i=2:20
if na(1,i)
n=a(1,i)
end
end
程序如下：
clear all
n=1:100;
b=1./(n.*n);
c=sum(b);
disp([当n等于100时：,num2str(c)])
-6.0000 26.6667 27.3333
4.求下列矩阵的秩、迹、条件数。
??1?12
3?（1）
a??
51?42???
3052?? ?111509??
（2）b???0.43432?
??8.9421??答：(1)
rank(a)；秩ans =
4
cond(a)；条件数ans =
11.1739
trace(a)；迹ans =
matlab第三章答案
【篇一：matlab第三章课后部分答案】
盘输入一个三位整数，将它反向输出。如输入639，输出为936
程序如下：
m=input(请输入一个三位整数：);
m1=fix(m/100);%求m的百位整数
m2=rem(fix(m/10),10);%求m的十位数字
m3=rem(m,10);%求m的个位数字
b=1./(4.^n);
c=sum(b);
disp([当n等于100时：,num2str(c)])
n=1:1000;
b=1./(4.^n);
c=sum(b);
disp([当n等于1000时：,num2str(c)])
n=1:10000;
b=1./(4.^n);
c=sum(b);
disp([当n等于10000时：,num2str(c)])
结果：当n等于100时：1.5669
当n等于1000时：1.5704
当n等于10000时：1.5708
3-9若两个连续自然数的乘积减1是素数，则称这两个连续自然数是亲密数对，该素数是亲密素数。例如，2*3-1=5，由于5是素数，所以2和3是亲密数对，5是亲密素数。求[2,50]区间内：
（1）亲密数对的对数。
case num2cell(60:69)
disp(d成绩合理);
case num2cell(0:59)
disp(e成绩合理);
otherwise
disp(成绩错误);
end
3-4建立5*6矩阵，要求输出矩阵第n行元素。当n值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。
程序如下：
clear all
y1=2*x1-0.5;
x2=sin(3*t).*cos(t);
y2=sin(3*t).*sin(t);
plot(x1,y1,r-,x2,y2,g-)
axis([-1,2源自文库-1.5,1])
hold on
s=solve(y=2*x-0.5,x=sin(3*t)*cos(t),y=sin(3*t)*sin(t));
结果：当n等于100时：0.33333
当n等于1000时：0.33333
当n等于10000时：0.33333
程序如下：
clear all
n=1:100;
b=((2*n).*(2*n))./((2*n-1).*(2*n+1));
c=prod(b);
disp([当n等于100时：,num2str(c)])
a=fix(rand(5,6)*89)+10;
n=input(请输入行数n（输出结果是5*6矩阵的元素）：);
if n5|n0
b=a(5:5,:);
disp([输入行数错误，程序输出最后一行：,num2str(b)])
else
c=a(n:n,:);
disp([输出结果为：,num2str(c)])
end
（2）与上述亲密数对对应的所有亲密素数之和。
程序如下：
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),)])
n=1:10000;
b=1./(4.*n-3)-1./(4.*n-1);
c=sum(b);
disp([当n等于10000时：,num2str(c)])
结果：当n等于100时：0.78415
当n等于1000时：0.78527
当n等于10000时：0.78539
程序如下：
clear all
n=1:100;
if sum3(sum-1/(2*n-1))3
disp(n-1)
disp(sum-1/(2*n-1))
end
end
方法2
clear all
n=1;
y=0;
while(y3)
x=1/(2*n-1);
n=1+n;
y=y+x;
end
disp([y3时n的最大值是：,num2str(n-2)])
disp([相应的y值是：,num2str(y-x)])
m=m3*100+m2*10+m1%反向输出m
3-3输入一个百分制成绩，要求输出成绩等级a,b,c,d,e。其中90~100分为a，80~89分为b，70~79分为c，60~69分为d，60分以下为e。
要求：（1）分别用if语句和switch语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。程序如下：
p=p+1;
h=[h,c];
end
l=0;
end
sum=sum(h);
disp([亲密数对的个数是：,num2str(p)])
disp([亲密数对的和是：,num2str(sum)])
（1）y3时的最大n值。
（2）与（1）的n值对应的y值。
程序如下：sum=0;
for n=1:100
sum=sum+1/(2*n-1);
theta=0:0.01:2*pi;
a=input(请输入a的值：);
b=input(请输入b的值：);
n=input(请输入n的值：);
rho=a*sin(b+n*theta);
polar(theta,rho,k);
end
下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。
对a的值进行改变：对比发现a只影响半径值的整倍变化
对b的值进行改变：对比发现b的值使这个圆转换了一定的角度
对n的值进行改变：对比发现当n=2时有如下规律
1、当n为整数时，图形变为2n个花瓣状的图形
2、当n为奇数时，图形变为n个花瓣状的图形分别让n为2、3、4、
5
【篇三：matlab答案第三章】
下列操作的命令。（1）建立3阶单位矩阵a。
（2）建立5x6随机矩阵a，其元素为[100，200]范围内的随机整数。（3）产生均值为1，方差为0.2的500个正态分布的随机数。
n=1:1000;
b=1./(n.*n);
c=sum(b);
disp([当n等于1000时：,num2str(c)])
n=1:10000;
b=1./(n.*n);
c=sum(b);
disp([当n等于10000时：,num2str(c)])
结果：当n等于100时：1.635
当n等于1000时：1.6439
3-5产生20个两位随机整数，输出其中小于平均值的偶数。
程序如下：
x=fix(rand(1,20)*89)+10;
x1=fix(sum(x)/20);
disp([平均数是：,num2str(x1)])
m=(rem(x,2)==0xx1);
n=find(m);
disp([小于平均数的数是：,num2str(x(n))]);
的行列式的值，并验证a与a-1
是互逆的。答：
b=inv(a);a的逆矩阵c=det(a);a的行列式的值d=a*b e=b*a
d=e=同阶单位矩阵因此a与a-1是互逆的。
3.求下面线性方程组的解。
??
4x1?2x2?x3?2?3x1?x2?2x?3?10 ?12x1
?3x2?8
答：
a=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(a)*b x =
3-6输入20个数，求其中最大数和最小数。要求分别用循环结构和调用matlab的max函数、min函数来实现。
程序如下：
主函数
c=input(请输入20个数：)
disp(最大值是：)
m=max(c)
disp(最小值是：)
n=min(c)
max函数
function m=max(a[])
m=a(1,20);
3-12编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在命令文件中调用该函数。程序如下：
mul函数%用于求两个矩阵的乘积和点乘
function [x,y]=mul(a,b)
x=a*b;
y=a.*b;
【篇二：matlab程序设计第三章课后习题答案】
注两曲线交叉点。
t=0:0.01:pi;
x1=t;
16（2）略
5.求矩阵a的特征值和相应的特征向量。
?1
10.5?a???110.25? ??
?0.50.252??
答：
a=[1 1 0.5;1 1 0.25;0.5 0.25 2] a =
1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [v,d]=eig(a) v =
（4）产生和a同样大小的幺矩阵。（5）将矩阵a对角线的元素加30。答：
(1)a=eye(3)
(2) a=100+100*rand(5,6)
(3)a=1+sqrt(0.2)*randn(10,50) (4)b=ones(size(a)) (5)a+30*eye(size(a))
2.建立一个方阵a，求a的逆矩阵和a
else
disp(成绩错误);
end
（2）switch语句
c=input(请输入成绩：);
switch fix(c)
case num2cell(90:100)
disp(a成绩合理);
case num2cell(80:89)
disp(b成绩合理);
case num2cell(70:79)
disp(c成绩合理);
（1）if语句
c=input(请输入成绩：);
if c=90c=100
disp(a成绩合理);
elseif c=80c=89
disp(b成绩合理);
elseif c=70c=79
disp(c成绩合理);
elseif c=60c=69
disp(d成绩合理);
elseif c60
disp(e成绩合理);