函数极值的matlab编程实现
%不讨论区间的情况。
clear
clc
syms x
f=input('请输入所要求极值的函数：');
s=diff(f); %一阶导数
s2=diff(f,2); %二阶导数
h=double(solve(s)); %一阶导数为零的点可能就是极值点，注意是可能，详情请见高数课本
for i=1:length(h)
if subs(s2,x,h(i))<0
disp(['函数在' num2str(h(i)) '处取得极大值，极大值为' num2str(subs(f,x,h(i)))])
elseif subs(s2,x,h(i))>0
disp(['函数在' num2str(h(i)) '处取得极小值，极小值为' num2str(subs(f,x,h(i)))])
else
disp(['函数在' num2str(h(i)) '处二阶导数也为0，故在该点处函数可能有极大值、极小值或无极值']) %详情见高数课本
end
end
%讨论区间求函数极值的情况。
clear;
clc;
syms x
f=input('请输入所要求极值的函数：');
A=input('请输入所要求极值的函数的区间：');
b=[];
s=diff(f); %一阶导数
s2=diff(f,2); %二阶导数
h=double(solve(s)); %一阶导数为零的点---可能极值点
for i=1:length(h)%判断可能极值点在不在区间A中
if A(1)
0
disp(['函数在区间[' num2str(A(1)) ',' num2str(A(2)) ']上' num2str(b(i)) '处取得极小值，极小值为' num2str(subs(f,x,b(i)))])
else
disp(['函数在区间[' num2str(A(1)) ',' num2str(A(2)) ']上' num2str(b(i)) '处二阶导数也为0，故在该点处函数可能有极大值、极小值或无极值'])
end
end
end
关于
级数，无理数
，欧拉常数
的matlab编程实现
一、
，
其中
是整数
clc
clear;
m=input('请输入的求和个数')
syms n;
b=[];
for k=2:m;
b(k-1)=(pi^k)/eval(symsum(1/n^k,1,Inf)); %求级数
的收敛值
end
digits(100);
bb=vpa(b)'
二、计算无理数
n=input('输入所要计算的阶乘数n:');
a=1;b=1;c=1;
digits(30)
for k=1:n %计算阶乘
a=a/k;
b=b+a;
end
bb=vpa(b)
三、计算欧拉常数
clc;
clear;
syms k
q=[];
p=[];
m=input('请输入要求和的项数：')
for n=10:m
q(n-9)=eval(symsum(1/k,1,n))-log(n+1);
p(n-9)=eval(symsum(1/k,1,n))-log(n);
end
a=q';
b=p';
c=[a,b]
r=(a>b)
浙公网安备 33021202002483