Matlab期末复习（一）

1
1.

clear
syms n
F = (1 + 1/n)^n;
limit(F,n,inf)

函数解释：

limit(F,x,a);计算x—>a时符号表达式F=F(x)的极限值

clear
A = [11:19;21:29;31:39;41:49;51:59;61:69;71:79;81:89;91:99;];
k = rank(A);

函数解释：

rank(A);求矩阵A的秩
rank(A,tol);tol为给定误差

clear
syms  b a
x = -2:0.01:1;
y = x.^3 .* sin(x) + x.^2/3 + x.*cos(x);
plot(x,y)
grid on;
h = line([-2,1],[0,0]);
set(h,'Color','k')

b = @(a)a^3*sin(a) + a^2/3 + a*cos(a);
k1 = fzero(b,[-1,-0.5])
k2 = fzero(b,0)

答案：

k1 =

-0.7889

k2 =

0

函数解释

plot(X,Y);基本平面图形绘制函数
grid on/off;给坐标加上网格线
X = fzero(fun,x0); 参数fun表示的是一元函数，x0表示求解的初始值
x0:也可以是区间

clear
x = -6:0.01:6;
y = (-x - pi) .*(x < -pi) + sin(x).*(x >= -pi & x <= pi) + ((x - pi)/2).*(x > pi);
plot(x,y)
grid on

答案：

函数解释：

分段函数求解，可以直接 .* 参数x的取值范围

clear
x=0:0.1:1;
y=1./(1+x.^2);
sum1=sum(y);
s1=0.1*sum1
s2=0.1*trapz(y)
s3=quad(@(x)1./(1+x.^2),0,1)
disp(['sum 求得积分',blanks(3),'trapz 求得积分',blanks(3),'Simpson 求得积分'])
disp([s1,s2,s3])
disp(['sum pi近似值',blanks(3),'trapz pi近似值',blanks(3),'Simpson pi近似值'])
disp([4*s1,4*s2,4*s3])

解：

clear
A = [6 2 1 -1;2 4 1 0;1 1 4 -1;-1 0 -1 3];
b = [6 1 5 -5]';
k = rank(A);
x = A\b;

答案：

k =

     4

x =

    0.7906
   -0.3613
    0.8639
   -1.1152

解：

clear
syms  y x a b
y = (2/3) * x^(3/2);
diff(y,x);
y1 = sqrt(1 + diff(y,x)^2);
h = int(y1,x,a,b)

答案：

h =
 
(2*(b + 1)^(3/2))/3 - (2*(a + 1)^(3/2))/3

函数解释：

diff(f,'a');求符号微分,对A变量求微分
sqrt();求根号的函数
int(f,v,a,b);符号函数积分.对变量v在[a,b]上求积分

n = input('输入一个四位正整数:n =') ;
flag = [0,1] ;%falg(0)计算第i次的n的值，flag(1)记录第i+1次的n的值
while((flag(1) - flag(2)) ~= 0)%当计算的前后两次的结果相同时便终止
flag(1) = n ;
strva = num2str(n) ;
N = length(strva) ;
num = [] ;
for i = 1:1:N
    num = [num,str2num(strva(i))] ;
end
 
%从小到大排序
for i =1:1:N - 1
    for j = 1:1:N-i
        if (num(j) > num(j+1))
            tmp = num(j) ;
            num(j) = num(j+1) ;
            num(j+1) = tmp ;
        end
    end
end
 
%计算最大值与最小值
max = num(N) ;
min = num(1) ;
for i =1:1:N
    if(i < N)
       max = max * 10 + num(N - i) ;
    end
    if(i > 1)
        min = min * 10 + num(i) ;
    end
end
 
n = max -min ;%计算最大值与最小值之差
flag(2) = n ;
end
n;