Matlab练习题
0.1窗口操作
建立更大的魔方矩阵，对矩阵的各行和各列进行求和等操作。
源代码：
format compact
A=magic(12)
sum(A)
sum(A,2)
trace(A)
trace(rot90(A))
0.2削顶全波整流
正弦波为sin
y x
，画一削顶全波整流曲线，顶部的值为3/4。
源代码：
clear
x=0:5:720;
xx=x*pi/180;
y=sin(xx);
figure
y(y<0)=-y(y<0);
y(y>=3/4)=3/4;
plot(xx,y,xx,sin(xx),'--')
title('削顶的全波整流曲线','FontSize',16)
grid on
0.3无穷级数的渐近线
下面无穷级数右边需要取多少项才与左边值的误差小于310-？画出级数的渐近线。
2222
111...6123π=+++clear
n=1;
x=0;
xx=[];
while 1
y=1/n^2;
x=x+y;
xx=[xx,x];
if n>2000,return,end
if abs(x-pi^2/6)<1e-3,break,end
n=n+1;
end
n %显示n
figure
plot(xx,'x-')
text(length(xx),x,num2str(x),'FontSize',16)
title('级数的渐近线','FontSize',16)
grid on
hold on
plot([0,1000],[pi^2/6,pi^2/6],'r')
0.4方格中的魔方
根据魔方的阶数画方格，将魔方填入方格中。
clear
n=input('请输入魔方的阶数：');
a=0:90:360;
aa=a*pi/180;
aa=aa+pi/4;
x=sqrt(2)/2*cos(aa);y=sqrt(2)/2*sin(aa);
M=magic(n);
figure
title('方格中的魔方','FontSize',16)
axis equal
hold on
for i=1:n
for j=1:n
plot(x+j,y-i)
text(j,-i,num2str(M(i,j)));
end
end
问题：转成.jpg格式后，最上面的直线没有了。
0.5正多边形的画法
从键盘输入正多边形的边数，画正多边形，再加外接圆。
clear
n=input('请输入正多边形的边数:');
a=linspace(0,360,n+1);
b=0:5:360;
aa=a*pi/180;
bb=b*pi/180;
x=cos(aa);y=sin(aa);
figure
plot(x,y)
axis equal
grid on
hold on
plot(cos(bb),sin(bb))
title(['正',num2str(n),'边形及其外接圆'],'FontSize',16)
