1,下面为m文件代码
i=-0.1;
j=1;
if i==0
i=i+0.005;%此步为了避免fzero(f,0)解不出
end
while(i<=0.1)
f=@(x)x*sin(1/x);%创建函数句柄
x(j)=fzero(f,i);%用fzero函数求解i附近的零点
j=j+1;
i=i+0.005;
end
x%此时x为一向量
2,
>> length=@(x)5./cos(x)+10./sin(x)
length =
@(x)5./cos(x)+10./sin(x)
>> [x,fval]= fminbnd(length,0,pi/2)
x =
0.8999
fval =
20.8097
>>length2=@(x)sqrt((5./cos(x)+10./sin(x))^2+10^2)
>> [x,fval]= fminbnd(length2,0,pi/2)
x =
0.8999
fval =
23.0877