實驗五 函數(shù)文件
第一題
function?y=mat1(x)????%建立函數(shù)文件mat1.m
y=[exp(x),log(x),sin(x),cos(x)];
%在命令窗口調(diào)用上述函數(shù)文件：
y=mat1(1+i)
%第二題程序一
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);
%在命令窗口調(diào)用該函數(shù)文件:
m1=input('m1=');
m2=input('m2=');
theta=input('theta=');
[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)
%第二題程序二
function?X=mat2(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;
%在命令窗口調(diào)用該函數(shù)文件：
X=mat2(1,1,60)
%第三題
function?flag=mat3(x)
flag=1;
for?i=2:sqrt(x)
if?rem(x,i)==0
????flag=0;
????break;
end
end
%在命令窗口調(diào)用該函數(shù)文件：
for?i=10:99
j=10*rem(i,10)+fix(i/10);
if?mat3(i)&mat3(j)
???disp(i)
end
end
%第四題
function?y=fx(x)
y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);
%在命令窗口調(diào)用該函數(shù)文件：
y=fx(2)
a=[1,2;3,4];
y=fx(a)
%第五題
%(1)
function?f1=mat5(n)
f1=n+10*log(n*n+5);
%在命令窗口中調(diào)用該函數(shù)文件：
y=mat5(40)/(mat5(30)+mat5(20))
%(2)方法一
function?f2=mat6(n)
f2=0;
for?i=1:n
????f2=f2+i*(i+1);
end
%在命令窗口中調(diào)用該函數(shù)文件如：
y=mat6(40)/(mat6(30)+mat6(20))
%(2)方法二
function?f2=mat7(n)
i=1:n;
m=i.*(i+1);
f2=sum(m);
end
%在命令窗口中調(diào)用該函數(shù)文件如：
y=mat7(40)/(mat7(30)+mat7(20))
實驗六 高層繪圖操作
第一題：
x=linspace(0,2*pi,101);
y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);
plot(x,y)
%第二題：
%（1）
x=linspace(-2*pi,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
plot(x,y1,'b-',x,y2,'r:',x,y3,'y--');
text(4,16,'\leftarrow?y1=x^2');
text(6*pi/4,-1,'\downarrow?y2=cos(2*x)');
text(-1.5*pi,-2.25*pi*pi,'\uparrow?y3=y1*y2');
%（2）
x=linspace(-2*pi,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
subplot(1,3,1);%分區(qū)
plot(x,y1);
title('y1=x^2');%設(shè)置標題
subplot(1,3,2);
plot(x,y2);
title('y2=cos(2*x)');
subplot(1,3,3);
plot(x,y3);
title('y3=x^2*cos(2*x)');
%（3）
x=linspace(-2*pi,2*pi,20);
y1=x.^2;
subplot(2,2,1);%分區(qū)
bar(x,y1);
title('y1=x^2的條形圖');%設(shè)置標題
subplot(2,2,2);
stairs(x,y1);
title('y1=x^2的階梯圖');
subplot(2,2,3);
stem(x,y1);
title('y1=x^2的桿圖');
subplot(2,2,4);
fill(x,y1,'r');%如果少了'r'則會出錯
title('y1=x^2的填充圖');
%其他的函數(shù)照樣做。
%第三題
x=-5:0.01:5;
y=[];%起始設(shè)y為空向量
for?x0=x
if?x0<=0??%不能寫成x0=<0
y=[y,(x0+sqrt(pi))/exp(2)];????%將x對應的函數(shù)值放到y中
else
y=[y,0.5*log(x0+sqrt(1+x0^2))];
end
end
plot(x,y)
%第四題：
a=input('a=');
b=input('b=');
n=input('n=');
t=-2*pi:0.01:2*pi;
r=a*sin(b+n*t);
polar(t,r)
%第五題
x=linspace(-5,5,21);
y=linspace(0,10,31);
[x,y]=meshgrid(x,y);%在[-5,5]*[0,10]的范圍內(nèi)生成網(wǎng)格坐標
z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);
subplot(2,1,1);
surf(x,y,z);
subplot(2,1,2);
contour3(x,y,z,50);%其中50為高度
的等級數(shù)，越大越密
%第六題
ezsurf('cos(s)*cos(t)','cos(s)*sin(t)','sin(s)',[0,0.5*pi,0,1.5*pi]);?%利用ezsurf隱函數(shù)
shading?interp??%進行插值著色處理
