p=prime(p) %再对对换后的素数矩阵找出所有的素数 运算结果：
p = 11 31 71 13 73 17 37 97 79 4. 设f(x)?11，编写一个MATLAB函数文件fx.m，使得调用f(x)?24(x?2)?0.1(x?3)?0.01时，x可用矩阵代入，得出的f(x)为同阶矩阵。
解： 函数fx.m文件： function f= fx(x) %fx fx求算x矩阵下的f(x)的函数值 A=0.1+(x-2).^2; B=0.01+(x-3).^4; f=1./A+1./B; 命令文件： clc; x=input('输入矩阵x='); f=fx(x) 运算结果：
>> x=input('输入矩阵x='); f=fx(x) 输入矩阵x=[7 2;12 5] f = 0.0437 10.9901 0.0101 0.1724 5. 已知y?f(40)
f(30)?f(20)(1) 当f(n)=n+10ln(n2+5)时，求y的值。
(2) 当f(n)=1×2+2×3+3×4+...+n×(n+1)时，求y的值。 解：(1) 函数f.m文件: function f=f(x) f=x+10*log(x^2+5); 命令文件： clc; n1=input('n1='); n2=input('n2='); n3=input('n3='); y1=f(n1); y2=f(n2); y3=f(n3); y=y1/(y2+y3) 运算结果如下： n1=40 n2=30 n3=20 y = 0.6390 (2). 函数g.m文件 function s= g(n) for i=1:n g(i)=i*(i+1); end s=sum(g); 命令文件： clc; n1=input('n1='); n2=input('n2='); n3=input('n3='); y1=g(n1); y2=g(n2); y3=g(n3); y=y1/(y2+y3)
运算结果如下：
n1=40 n2=30 n3=20 y = 1.7662
实验六 高层绘图操作
一、实验目的
1. 掌握绘制二维图形的常用函数。 2. 掌握绘制三维图形的常用函数。 3. 掌握绘制图形的辅助操作。 二、实验内容
1. 设y??0.5???3sinx?cosx，在x=0~2π区间取101点，绘制函数的曲线。 2?1?x?解：M文件如下：
clc; x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)); plot(x,y) 运行结果有：
2. 已知y1=x，y2=cos(2x)，y3=y1×y2，完成下列操作： (1) 在同一坐标系下用不同的颜色和线型绘制三条曲线。 (2) 以子图形式绘制三条曲线。
(3) 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 解：（1） M文件：
2
clc; x=-pi:pi/100:pi; y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'k--') 运行结果：
（2）M文件：
clc; x=-pi:pi/100:pi; y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1); plot(x,y1,'b-'); title('y1=x^2'); subplot(1,3,2); plot(x,y2,'r:'); title('y2=cos(2x)'); subplot(1,3,3); plot(x,y3,'k--'); title('y3=y1*y2'); .运行结果：
（3）M文件：