Evaluating a function with three variables

3 次查看（过去 30 天）
Alin Brad
Alin Brad 2018-9-26
评论： Walter Roberson 2018-9-28
Dear all, Since I am new in MatLab, I want to evaluate the following equation y=(1./((x.^2)+1))-(p./((x+w).^2+1))-(p./((x-w).^2+1)); I want to sweep the variables x,p and w as follows x from -5 to 5 with 100 steps ,p from 0 to 1 with 100 steps ,w from 0 to 3 with 100 steps and then I want to normalized the functiony and then I want to extract the specific values of x (lets name new variable r) which give me specific value of y (lets say y=0.5) for each combination of variables x,w,p. After that I want to make a surface graph of ( r,w and p)
thanks

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

采纳的回答

KSSV
KSSV 2018-9-26
编辑：KSSV 2018-9-26
N = 100 ;
x = linspace(-5,5,N) ;
p = linspace(0,1,N) ;
w = linspace(0,3,N) ;
[x,p,w] = meshgrid(x,p,w) ;
y=(1./((x.^2)+1))-(p./((x+w).^2+1))-(p./((x-w).^2+1));
figure
hold on
for i = 1:N
surf(x(:,:,i),p(:,:,i),w(:,:,i),y(:,:,i))
end
% GEt manually y = 0.5
idx = abs(y-0.5)<=10^-3 ;
xr = x(idx) ;
pr = p(idx) ;
wr = w(idx) ;
figure
scatter3(xr,pr,wr,10,wr)
% use isosurface
p = patch(isosurface(x,p,w,y,0.5));
p.FaceColor = 'red';
p.EdgeColor = 'none';
daspect([1 1 1])
view(3);
axis tight
camlight
lighting gouraud
  6 个评论
显示 4更早的评论
Alin Brad
Alin Brad 2018-9-26
I have made the modification and change the value of y to idx = (max(y(:))-(abs(min((y(:)))))/2) however, I got one value of idx and now values of xr and got the following message #Array indices must be positive integers or logical values#.
Walter Roberson
Walter Roberson 2018-9-26
You would not use val as idx. You would use
idx = abs(y-val)<=10^-3 ;

请先登录，再进行评论。

更多回答（1 个）

Alin Brad
Alin Brad 2018-9-26
I have made this modification to make surface and make xr as Z dimension , but it doesnot work properly
N = 100 ;
x = linspace(-5,5,N) ;
p = linspace(0,1,N) ;
w = linspace(0,3,N) ;
[x,p,w] = meshgrid(x,p,w) ;
y=(1./((x.^2)+1))-(p./((x+w).^2+1))-(p./((x-w).^2+1));
figure
hold on
for i = 1:N
surf(x(:,:,i),p(:,:,i),w(:,:,i),y(:,:,i))
end
% GEt manually y
val=(max(y(:))-(abs(min((y(:)))))/2)
idx = abs(y-val)>=10^-3 ;
xr = x(idx) ;
pr = p(idx) ;
wr = w(idx) ;
figure
scatter3(xr,pr,wr,10,wr)
% use isosurface
p = patch(isosurface(xr,pr,wr,y,val));
p.FaceColor = 'red';
p.EdgeColor = 'none';
daspect([1 1 1])
view(3);
axis tight
camlight
lighting gouraud
  4 个评论
显示 2更早的评论
Alin Brad
Alin Brad 2018-9-26
Thank you The bellow program returns a message The size of X must match the size of V or the number of columns of V. and show a plot of p as z axis which should be xr as z axis
N = 100 ;
x = linspace(-5,5,N) ;
p = linspace(0,1,N) ;
w = linspace(0,3,N) ;
[x,p,w] = meshgrid(x,p,w) ;
y=(1./((x.^2)+1))-(p./((x+w).^2+1))-(p./((x-w).^2+1));
figure
hold on
for i = 1:N
surf(x(:,:,i),p(:,:,i),w(:,:,i),y(:,:,i))
end
% GEt y
val=(max(y(:))-(abs(min((y(:)))))/2)
idx = abs(y-val)>=10^-3 ;
yt = y;
yt(~idx) = nan;
xr = x(idx) ;
pr = p(idx) ;
wr = w(idx) ;
figure
scatter3(xr,pr,wr,10,wr)
% use isosurface
p= patch(isosurface(xr, pr, wr, yt, val));
p.FaceColor = 'red';
p.EdgeColor = 'none';
daspect([1 1 1])
view(3);
axis tight
camlight
lighting gouraud
Walter Roberson
Walter Roberson 2018-9-28
Notice in my call I did not select subsets of x and so on for the isosurface call: I left them the original size but assigned nan to locations in the grid that were not of interest.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Volume Visualization 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by