2-D graph help
显示 更早的评论
Hi everyone. I would like to produce a graph like the one shown in the picture on MatLAB. Any help with writing the script will be greatly appreciated. The radius of each band increases from the central point 0,0.
采纳的回答
Star Strider
2016-2-10
One possibility:
w = linspace(0, 2*pi);
c = [cos(w); sin(w)];
r = [1 2 3];
figure(1)
plot(r(1)*c(1,:), r(1)*c(2,:))
hold on
plot(r(2)*c(1,:), r(2)*c(2,:))
plot(r(3)*c(1,:), r(3)*c(2,:))
hold off
axis equal
5 个评论
Ogen
2016-2-10
Star Strider
2016-2-10
My pleasure.
To change the radii of the circles, change the values in my ‘r’ vector. For example, to put the second and third closer together:
r = [1 2 2.3];
You can make the values in ‘r’ whatever you want. To change the other characteristics of the circles (line widths, colours, and other properties), see the documentation for the plot function.
Ogen
2016-2-10
Ogen
2016-2-10
Star Strider
2016-2-10
My pleasure.
The mushroom cloud problem would seem to be a computational fluid dynamics problem far beyond the scope of my knowledge. However, if you want to plot a cone plotting the radius at different heights, and perhaps the shape of the cloud itself, I would investigate the cylinder function. It, possibly in combination with the sphere function, could give you the result you want. This code uses both:
[Xc,Yc,Zc] = cylinder(0.5, 20);
[Xs,Ys,Zs] = sphere(20);
Xsm = Xs * 2;
Ysm = Ys * 2;
Zsm = Zs + 2;
figure(1)
surf(Xc,Yc,Zc)
hold on
surf(Xsm,Ysm,Zsm)
hold off
grid on
axis equal
view([-35 20])
to produce:
更多回答（2 个）
Ogen
2016-2-10
0 个投票
5 个评论
Star Strider
2016-2-10
It is complicated, with rapid convection and all the rest of the atmospheric physics involved. I’m not certain what you want, so I’ll leave the web search to you. The same sort of physics might apply to the rapid development of a thunderstorm, so there are likely articles on the Internet that you could implement in MATLAB that would give you the behaviour you want to model. I imagine it would involve solving partial differential equations in time and distance, but there could be published solutions for you to use.
My plot is simply a way to create the image. You can easily experiment with my code to make it look the way you want, and also to animate it. See the documentation on Vector Fields and Surfaces, Volumes, and Polygons, particularly the section on Volume Visualization for more plotting options.
When you get your project coded and working, write it up and submit it as a File Exchange contribution. It seems to me to be something that atmospheric physicists and other fluid dynamicists would be interested in.
Cheers — and good luck! We’ll help you here as we can!
Ogen
2016-3-6
Star Strider
2016-3-6
I’m glad you got it working! Consider contributing it to the File Exchange when you get it to the point you’re happy with it, and have the time to document it and make it into a funciton.
Ogen
2016-3-6
Star Strider
2016-3-6
My pleasure.
The drawnow function is useful inside a loop, so that you calculate and draw each frame and animate it in real time. This works best if the calculations creating the image are reasonably fast, and the resolution relatively low.
The getframe function and its friends allow you to create a movie out of your MATLAB images. This is best if the calculations for each frame take a while, are relatively high resolution, and so are not suited to real time display. If you have detailed images that take some time to render and you don’t need to change the parameters of the simulation in real time, this is likely the better option.
It all depends on what you want to depict. I last did this sort of thing many years ago when MATLAB, and PCs in general, were slower, so getframe was the best option for me. I've used drawnow recently with good results, but with simpler code.
Ogen
2016-2-10
0 个投票
4 个评论
Star Strider
2016-2-10
I believe this is what you want:
[Xc,Yc,Zc] = cylinder(0.5, 20);
Xcm = Xc;
Ycm = Yc;
Zcm = 2 * Zc; % <— Multiplied By ‘2’ To Double Cylinder Height
[Xs,Ys,Zs] = sphere(20);
Xsm = Xs * 2;
Ysm = Ys * 2;
Zsm = Zs + 3; % <— Added ‘1’ To Move Spheroid Up
figure(1)
surf(Xcm,Ycm,Zcm)
hold on
surf(Xsm,Ysm,Zsm)
hold off
grid on
axis equal
view([-35 20])
If it isn’t, let me know.
Note the changes I initially made to the original sphere matrices to distort it. You can do similar things with the cylinder radius. You can make it look like a bottle with the appropriate radius vector. It has to retain its essential cylinder properties, but it doesn’t have to just be a cylinder if you don’t want it to.
Ogen
2016-2-10
Star Strider
2016-2-11
As always, my pleasure.
Star Strider
2016-2-11
UPDATE — I did some Web searching, and discovered Dr. Jos Stam’s website at the University of Toronto http://www.dgp.toronto.edu/people/stam/reality/Research/pub.html. See specifically ‘Real-Time Fluid Dynamics for Games’. There’s also a CDROM, but I doubt the code’s written in MATLAB, so you’d have to adapt it. I found it through Mushroom Cloud Physics on ‘gamedev.net’. (I’m a gamer, so I can’t explain this not occurring to me earlier, other than that I do primarily medieval-themed games where mushroom clouds are rare and dragons are plentiful.)
If I find anything else, I’ll post it here.
类别
在 帮助中心 和 File Exchange 中查找有关 Lighting, Transparency, and Shading 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
