matlab上机报告绘制函数,Matlab程序设计上机报告

y=(exp(0.3*a(i))-exp(-0.3*a(i)))/2*sin(a(i)+0.3)+log((0.3+a(i))/2) End

5.当n 分别取100、1000、10000时，求下列格式的值：要求分别用顺序结构和向量运算(使用sum 函数)来实现。

(1))6

...(1 (3121112)

2222π=+++++n

>> A=[0,0,0];

n=[100,1000,10000]; for i=1:3 for j=1:n(i)

A(i)=1/(j^2)+A(i); end end A A =

1.6350 1.6439 1.6448 >> A=[0,0,0];

n=[100,1000,10000]; for i=1:3 j=1:1:n(i);

A(i)=sum(1./j.^2); end A A =

1.6350 1.6439 1.6448 (2)

???

?

?=???? ??+-??? ??????? ??????? ????2...)12)(12()2)(2(...756653443122πn n n n >> A=[1,1,1];

n=[100,1000,10000]; for i=1:3 for j=1:n(i)

A(i)=A(i)*((2*j)^2)/((2*j-1)*(2*j+1)); end end >> A

A =

1.5669 1.5704 1.5708

6.建立5×6矩阵，要求输出矩阵第n 行元素。当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。 >> A=rand(5,6)

n=input('输入第几行？'); if n>5 n=5;

disp('超出范围，输出最后一行的元素'); A(n,:) else A(n,:) end

第4章 1.绘制下列曲线

(1)!

33

x x y -=

>> x=0:pi/100:2*pi; >> y=x-x.^3/factorial(3);

>> plot(x,y)

(2)64222=+y x

>> ezplot('x.^2+2*y.^2-64',[-8,8,-6,6])