查看: 244|回复: 1|关注: 0

[未答复] 有积分的隐函数问题

[复制链接]

新手

13 麦片

财富积分


050


6

主题

15

帖子

0

最佳答案
发表于 2020-8-10 10:19:48 | 显示全部楼层 |阅读模式
麻烦各位大神帮忙看一下
我的目标方程是比较复杂的一个含有积分的三维隐函数,目的是通过画出三维图像来得到z的一系列值
曾经尝试着用了for循环来求解方程值,但是由于没有matlab基础,经过高人指点之后效果还是不如意

也尝试了不用for循环  但又一直报错

实在是不知道该怎么解决了  特此请求大家的帮忙  (有什么表达不清楚的地方还请见谅)
  1. clc;
  2. clear;
  3. close all;
  4. a = linspace(-1.5,0,10);
  5. b = linspace(-3.5,0,10);
  6. c = linspace(0,5,10);
  7. [x,y,z]=meshgrid(a,b,c);
  8. v= -7.93.*10.^14.*2./sqrt(pi).*integral(@(r)((((r-1.23)./0.025852).^0.5)./(1+exp((r-z)./0.025852))),1.23,15)...
  9.     +1.57.*10.^17.*exp(-z./0.025852)...
  10.     +5.8917.*10^21.*2.*exp((0.9314+y-x)./0.025852).*exp((2.0354-z)./0.025852)./(2+exp((2.0354-z)./0.025852))...
  11.     +2.*2.9458.*10^21.*6.*exp((x-1.8411)./0.025852).*exp((3.506-2.*z)./0.025852)./(6+exp((3.506-2.*z)./0.025852))...
  12.         +2.*2.3567.*10^22.*6.*exp((0.25.*x+0.5.*y-0.5006)./0.025852).*exp((3.4079-2.*z)./0.025852)./(6+exp((3.4079-2.*z)./0.025852));      
  13. p=patch(isosurface(x,y,z,v,0));
  14. p.FaceColor='none';
  15. p.EdgeColor='red';
  16. daspect([1,1,1])
  17. view(3);axis tight
  18. grid on;
  19. xlabel('x'),ylabel('y'),zlabel('z');
  20. view(45,30)
  21. F = p.Faces;
  22. V = p.Vertices;
复制代码


回复主题 已获打赏: 0 积分

举报

新手

13 麦片

财富积分


050


6

主题

15

帖子

0

最佳答案
 楼主| 发表于 2020-8-11 16:42:47 | 显示全部楼层
:(有大神帮忙看一下的么?
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表