matlab考试.docVIP

2012年Matlab程序设计考核题目 一、函数文件的编写 1．已知阶勒让德多项式的表达式为，其中表示不超过的最大整数，试编写双输入变量单输出变量函数my_legendre，当给定和的值时，由该函数可求出勒让德多项式的值，其中为整数，. 注意：此函数还应该支持关于的数组运算，即当输入变量为数组时，返回变量也应是与同维的数组。此外利用此函数在同一坐标系内绘制的图像。 function P=my_legendre(x,l) if nargin~=2 error(要求双输入变量); end if nargout1 error(输出变量个数过多); end if abs(l-floor(l))eps||l0 error(要求为整数); end if x1||x0 error(x要求在0-1之间); end M=l/2; m=fix(M); P=0; for k=1:m; P=P+((-1)^k).*(factorial(2*l-2*k))./(2^l.*(factorial(k)).*(factorial(l-k)).*(factorial(l-2*k))).*x^(l-2*k); end function y=my_legendre(l,x) %求勒让德多项式 x=0:0.01:1;%x的取值范围及步长 y0=legendre(0,x); y1=legendre(1,x); y2=legendre(2,x); y3=legendre(3,x); y4=legendre(4,x); y5=legendre(5,x); y6=legendre(6,x); y7=legendre(7,x);%分别列出0到7的表达式 plot(x,y0(1,:),x,y1(1,:),x,y2(1,:),x,y3(1,:),x,y4(1,:),x,y5(1,:),x,y6(1,:),x,y7(1,:)) %作图 title(Legendre) %图像命名 2．在量子力学中，粒子的自旋量子数只能有两种可能的取值，一种是正整数与零，如，另一种称为半整数，如。对于一个确定的，其磁量子数有个取值，即，相应自旋算符为维方阵，其矩阵元为 试编写单输入变量三输出变量函数spinam(j)，用来求自旋为的粒子的角动量算符的矩阵形式。 function A=spinam(j); if nargin~=1 error(输入变量个数错误，只能有一个输入变量) end if nargout~=3 error(输出变量个数错误，只能有三个输出变量) end if abs(2.*j-floor(2.*j))eps | j0 error(j只能取0，正整数，半整数) end for m=-j:j for n=-j:j if m==n F(m,n)=1; else F(m,n)=0; end A(m,n)=[m.*F(m,n);1/2.*sqrt(j(j+1)-m(m+1)).*F(m+1,n)+1/2.*sqrt(j(j+1)-m(m-1)).*F(m-1,n);1/(2*i).*sqrt((j(j+1)-m(m+1)).*F(m+1,n)-1/(2*i).*sqrt(j(j+1)-m(m-1)).*F(m-1,n))] end end 3．玻色-爱因斯坦积分的定义为： ， 其中，根据以上信息编写函数双输入单输出函数BEint(n,z)用来求函数的值，并利用此函数绘制的图像。 function plotBEint() % plotBEint函数调用BEint函数求值并画出函数的图形 clear z=0:0.1:1; %自变量x f=0*z; m=length(z); for n=1:m %循环计算gnz的值 g(n)=BEint(3/2,z(n)); end plot(z,double(g)) %画出函数图形 function gnz=BEint(n,z) % BEint函数功能是求玻色-爱因斯坦积分gnz的值 if nargin~=2 error(输入参数个数不正确，应为两个输入参数); end if nargout~=1 error(输出参数个数不正确，应为1个输出参数); end [mn,nn]=size(n); if(mn~=1 | nn~=1) error(输入数据n应为标量); end [mz,nz]=size(z); if(mz~=1 | nz~=1) error(输入数据z应为标量); end syms x %定义符号