syms y
x=linspace(-2,2,200)
n=length(x)
g=[];
for i=1:n
a=x(i);
t=vpa(solve(y/(1+a^2+y^2)-sin(a+cos(y)),y),3);
g=[g t];
end
对于这个代码主要是想算出对应x的y,因为直接用get函数得到的x和y是按照顺序先后排列而不是一一对应。下面那个代码得出的结果是错误的。
syms x y
f=y/(x^2+y^2+1)-sin(x+cos(y))
ezplot(f,[-2 2],[-2 2])
h=get(gca,'Children');
x=get(h,'xdata');
y=get(h,'ydata');
c=x.^2+y^2;%解离原点距离 *这边如果y不转置根本运行不了*
n=find(c==min(c));%找最短
x(n)%结果
y(n)
但是可以用第二个 代码的find去做。
求min(x^2+y^2)
方程在实数区域没有根,可能在虚数区域是有根的
F = @(y)y/(1+a^2+y^2)-sin(a+cos(y));
syms yr yi
H = F(yr+yi*1i);
HH = matlabFunction([real(H),imag(H)],'vars',{[yr,yi]});
fsolve(HH,[-.2,1.2])
这个就将根的形式写成复数的形式。
fsolve函数主要通过迭代法求根的,如果你给的是复数形式,那么就有可能得到复数形式的根
clear all;clc
syms x y
c=[];
for y=linspace(-2,2,200)
c=[c solve(y/(1+x^2+y^2)-sin(x+cos(y)),x)];
end
y=linspace(-2,2,200);x=c;
d=[x' y' (y.^2+x.^2)'];
d=double(d);
[i j]=find(d==min(d(:,3)));
d(i,:)
gg=d(i,3);
cc=sqrt(gg);
h=ezplot('y/(1+x^2+y^2)-sin(x+cos(y))',1000);
hold on;
rectangle('Position',[-cc,-cc,2*cc,2*cc],'Curvature',[1,1]),axis equal
hold on;
aa=ones(size(x))*0;
plot(x,aa);
hold on;
bb=ones(size(y))*0;
plot(bb,y);