how to Shading the intersection area?

2 次查看（过去 30 天）

sadeem alqarni 2018-9-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/421511-how-to-shading-the-intersection-area

评论： sadeem alqarni 2018-10-1

x0 = -5;
 x1 = 5;
 Nx = 301;
 % Specify y range and number of points
 y0 = -5;
 y1 = 5;
 Ny = 301;
 % Construct mesh
 xv = linspace(x0,x1,Nx);
 yv = linspace(y0,y1,Ny);
 [x,y] = meshgrid(xv,yv);
 % Calculate z
 q = x + i*y;
 % 2nd order Runge-Kutta growth factor
   g1=((196.*q.^4 - 108.*q.^3 - 1143.*q.^2 + 3888.*q + 5184).^(1/2) - 27.*q + 14.*q.^2 + 72)./(12.*(q.^2 - 6.*q + 12))
   g2=-(27.*q + (196.*q.^4 - 108.*q.^3 - 1143.*q.^2 + 3888.*q + 5184).^(1/2) - 14.*q.^2 - 72)./(12.*(q.^2 - 6.*q + 12))
   % Calculate magnitude of g
   gmag1 = abs(g1);
   gama2=abs(g2)
   % Plot contours of gmag
   contourf(x,y,gmag1,[1 1],'k');
   hold on
   contourf(x,y,gama2,[1 1],'r');
   hold on
   axis([x0,x1,y0,y1]);
   axis('square');
   xlabel('Real \lambda\Delta t');
   ylabel('Imag \lambda\Delta t');
   grid on;
   figure 
   contour(x,y,[gama2,gmag1],[1 1],'r')

3 个评论
显示 1更早的评论隐藏 1更早的评论

KSSV 2018-10-1

You need to explain more....question is not clear.

sadeem alqarni 2018-10-1

   iwant to  Shading the intersection area for g1,g2 ?
  g1=((196.*q.^4 - 108.*q.^3 - 1143.*q.^2 + 3888.*q + 5184).^(1/2) - 27.*q + 14.*q.^2 + 72)./(12.*(q.^2 - 6.*q + 12))
   g2=-(27.*q + (196.*q.^4 - 108.*q.^3 - 1143.*q.^2 + 3888.*q + 5184).^(1/2) - 14.*q.^2 - 72)./(12.*(q.^2 - 6.*q + 12))

请先登录，再进行评论。

请先登录，再回答此问题。