[未答复] plot3 ,for语句,图像与函数不符

[复制链接]
fanfanfanfxx 发表于 2018-1-30 01:49:37

想用plot3来画函数图像,出来的图跟函数不太相符,不知道是for语句用错了,还是函数输入有问题,检查了几遍,自己没发现函数输入的问题,不知道问题出在哪里,请帮忙看看,谢谢啦。

  1. clear
  2. clc
  3. CI=200;
  4. CP=10;
  5. d=1000;
  6. ci=2;
  7. cp=10;
  8. hu=4;
  9. hs=3;
  10. eta=0.2;
  11. delta=0.0002;
  12. y=0.4;
  13. np=1;
  14. M=zeros(1,400);
  15. Q=zeros(1,400);
  16. TC=zeros(400,400);



  17. for i=1:400
  18.     M(1,i)=0.6+0.0002*i;
  19. end
  20. for j=1:400
  21.     Q(1,j)=j;
  22. end
  23. for i=1:400
  24.     for j=1:400
  25.    
  26.       
  27.        TC(1,i)=-2*eta*(-1+M(1,i)+y)/(delta*(-1+M(1,i)-y))+ d*(CI+np*CP)/ Q(1,j)-(hs*Q(1,j)*y^2/(2*np))*((-3+3*M(1,i)+y)/(-1+M(1,i)+3*y))-cp*d*(M(1,i)-1)+ci*d*M(1,i)+hu*Q(1,j)*M(1,i)/2+hs*Q(1,j)*M(1,i)^2/2;
  28.       
  29.     end
  30. end

  31. set(gcf,'Position',[100 100 260 220]);
  32. set(gca,'Position',[.13 .17 .80 .74]);
  33. set(findobj(get(gca,'Children'),'LineWidth',0.5),'LineWidth',1);
  34. set(gca, 'Fontname', 'Times newman', 'Fontsize', 10.5);


  35. plot(M,TC,'g--');

  36. xlabel(' The value of M');
  37. ylabel('The  total cost');

  38. box off;

  39. xl=xlim;

  40. yl=ylim;

  41. line([xl(1),xl(2)],[yl(2),yl(2)],'color',[0 0 0]);   %画上边框,线条的颜色设置为黑色

  42. line([xl(2),xl(2)],[yl(1),yl(2)],'color',[0 0 0]);    %画右边框 ,线条的颜色设置为黑色

  43. surf(TC(1:400,1:400),'DisplayName','TC(1:400,1:400)');figure(gcf);hold on;

  44. A=min(min(TC));
  45. [x,y]=find(TC==min(min(TC)));
  46. plot3(y,x,A,'ro');

复制代码


QQ截图20180130014757.png
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表