数值分析课程设计 沈阳航空航天大学 课程设计用纸 附 录
附 录
源程序: （1.1.1）
x=-2:0.02:2;%从-2，增量为0.02 y=x; %从-2，增量为0.02
[X,Y]=meshgrid(x,y);%对X.Y取x,y的网格坐标矩阵 Z=zeros(201);%取Z为201*1个的零矩阵 for i=1:201 %利用for循环从1到201 for j=1:201%利用for循环从1到201 if X(i,j)+Y(i,j)>1
Z(i,j)=0.5*exp(-0.75*Y(i,j).^2-3.75*X(i,j).^2-1.5*X(i,j)); %如果在X+Y>1范围内将X(i,j), Y(i,j)利用相对应的表达式带入Z零矩阵相对应的点中
else if X(i,j)+Y(i,j)<=-1
Z(i,j)=0.5*exp(-0.75*Y(i,j).^2-3.75*X(i,j).^2+1.5*X(i,j)); %否则，如果在X+Y<=1范围内将X(i,j), Y(i,j) 利用相对应的表达式带入Z零矩阵相对应的点中
else Z(i,j)=0.8*exp(-Y(i,j).^2-6*X(i,j).^2);
%否则将X(i,j), Y(i,j) 利用相对应的表达式带入Z零矩阵其余的点中
end end
end %利用条件语句写出题中给的分段函数
mesh(X,Y,Z);%作出三维曲面图 end
（1.1.2.1）
x=-50:1:50; %从-50到50，增量为1 y=atan(x); %将带入表达式中 plot(x,y); %作出图像
（1.1.2.2）
function y=myfun(n) %建立变量为n的函数y x=-50:1:50 %从-50到50，增量为1 for k=1:n %用for循环从1到n p=0;
p=p+(-1)^(k-1)/(2*k-1).*x.^(2*k-1); %带入表达式
第 11 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 附 录
end
plot(x,p); %作出图像
（1.2）
function s=l(n) %建立变量为n的函数s s=0;
for k=1:n %对k从1一直取n到 s=s+1./k^2; %带入表达式 end s （1.3）
syms t %定义符号变量t
y=t^3+t^2; %带入t对y的表达式 x=t-log(1+t); %带入t对x的表达式 f1=diff(y)/diff(x); %求y对x的一阶导数. f2=diff(f1)/diff(x); %求y对x的二阶导数. f2
(1.4)
syms x y %定义符号变量x y f=log(x)+exp(-y/x); %写出f(x)表达式 diff(f,'x'); %求f对x的偏导数. diff(f,'y'); %求f对y的偏导数. z=(-1)* diff(f,'x')/diff(f,'y') %求y对x的导数
（1.5）：
syms x ; %定义符号变量x. f=4*(x.^3)-6*(x.^2)+x-2; df=diff(f,x);
x0=solve(df,x); %解函数df的零点. y1=4*(x0(1).^3)-6*(x0(1).^2)+x0(1)-2;
y2=4*(x0(2).^3)-6*(x0(2).^2)+x0(2)-2;%将零点分别带入表达式其结果为y1,y2 x=0:0.002:1;
b=4*(x.^3)-6*(x.^2)+x-2; plot(x,b)
hold on %作出原函数图像
plot(x,y1,'r') %将x0(1),y1带入切线表达式中 hold on %作出一条切线图像
plot(x,y2,'r') %将x0(2),y2带入切线表达式中 hold on %作出另一条切线图像 二
第 12 页
数值分析课程设计 沈阳航空航天大学 课程设计用纸 附 录
function ans=Untitled5(n,a,b,g)
%建立变量为n,a,b,g的函数ans; n为被区间被分个数；a,b为被积区间；g为被积函数 f=@(x)g(x) %定义句柄g为函数f(x) h=((b-a)/n); %计算h
m=0; %将0赋值给m
for i=1:n-1 %用for循环从1到n-1 k1=a+i*h; %将a+i*h赋值给k1
m=m+feval(f,k1); %调用函数f(x)在k1点的值，利用for循环计算m值 end
t=0; %将0赋值给t
for j=0:n-1 %用for循环从0到n-1
k2=a+(1/2+j)*h; %将a+(1/2+j)*h赋值给k2
t=t+feval(f,k2); %调用函数f(x)在k2点的值，利用for循环计算t值 end
p=feval(f,a)+feval(f,b); %调用函数计算f(a),f(b)的值 ans=h/6*(p+4*t+2*m); %带入公式计算积分值
第 13 页