| 查看: 2627 | 回复: 8 | ||
[求助]
继续求助Matlab绘图！
|
|
在帖子http://muchong.com/bbs/viewthread.php?tid=4346449中部分实现了绘图目标。感谢sogooda朋友！ 现在情况是：得到了520*3的矩阵，前两列是点的XY坐标，第三列是与圆心成的角度。 现在想要实现： 一、colorbar按矩阵行数从黑色[0 0 0]过渡到绿色[0 1 0]， 二、在图上划出形状，首先是质心与圆心连线和Y轴角度的值是第三列数据，然后画出图上所示的形状并填充采样时刻（所在的第几行）对应的颜色。 希望得到大家的建议与帮助！效果图和mat文件如下。  效果图[ Last edited by plxh on 2012-4-8 at 09:43 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab.mat
2012-04-08 09:42:48, 9.33 K
» 本帖@通知
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab画图求助? (有人说, 标题一定要长------------------------->)
已经有5人回复
- 求助matlab三维画图中：Inner matrix dimensions must agree.
已经有4人回复
- MATLAB 画图小问题
已经有11人回复
- 求助matlab画图问题
已经有3人回复
- 求助：谁会用matlab画图？
已经有13人回复
- 【求助】利用MATLAB求解的外围轮廓包络线为什么不与图形外围轮廓重合？怎样才可以精确
已经有8人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab 遗传算法工具箱绘图命令是什么？
已经有6人回复
2楼2012-04-08 12:02:00
3楼2012-04-08 12:04:33
|
大哥，如果要一张图上画出来原理完全一样但半径不一样的三个呢？ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 3个.bmp
- 附件 2 : matlab3.mat
2012-04-08 13:19:06, 2.52 M
2012-04-08 13:19:06, 30.89 K
4楼2012-04-08 13:20:16
xiegangmai: , 帖子中代码部分可以用“插入代码”功能 2012-04-08 22:52:58
|
clear;clc;close all load matlab3 for i=1:3; if i==1 %必须是用两个等号。 eta12=eta1(1:5:end,  ;x=eta12(:,2); y=eta12(:,1); % theta=eta1(:,3); % ? theta=pi/2-eta12(:,3); figure1=figure('color','w');plot(x,y);hold on;axis equal myr=linspace(0,0,length(x)); myg=linspace(0,1,length(x)); myb=linspace(0,0,length(x)); for jj=1:length(x) plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]); end colormap([myr', myg', myb']) colorbar hold on; elseif i==2 eta22=eta2(1:5:end, ;x=eta22(:,2); y=eta22(:,1); % theta=eta1(:,3); % ? theta=pi/2-eta22(:,3); figure1=figure('color','w');plot(x,y);hold on;axis equal myr=linspace(0,0,length(x)); myg=linspace(0,1,length(x)); myb=linspace(0,0,length(x)); for jj=1:length(x) plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]); end colormap([myr', myg', myb']) colorbar hold on; elseif i==3 eta32=eta1(1:5:end, ;x=eta32(:,2); y=eta32(:,1); % theta=eta1(:,3); % ? theta=pi/2-eta32(:,3); figure1=figure('color','w');plot(x,y);hold on;axis equal myr=linspace(0,0,length(x)); myg=linspace(0,1,length(x)); myb=linspace(0,0,length(x)); for jj=1:length(x) plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]); end colormap([myr', myg', myb']) colorbar hold on; end end function plotSquare2(x0,y0,theta,a,mycolor) xs=[-2*a,-2*a,2*a,3*a,2*a,-2*a]; ys=[a,-a,-a,0,a,a]; xy=[cos(theta),-sin(theta); sin(theta),cos(theta)]*[xs;ys] ; fill(xy(1, +x0,xy(2,+y0,mycolor)为什么会报警呢？ 报警： ??? Error: File: Untitled4.m Line: 70 Column: 1 Function definitions are not permitted at the prompt or in scripts. |
5楼2012-04-08 13:44:10
【答案】应助回帖
★ ★ ★ ★ ★ ★
xiegangmai: 金币+1, 谢谢参与！ 2012-04-08 22:53:07
plxh: 金币+5, ★★★很有帮助, 大哥你说得对，我就是赶鸭子上架，好高骛远了。昨天又学习了下function文件和M文件的不同。现在已基本实现了把三个1/4圆画到一张图上了。现在的问题就是我想把产生数据的m文件的内容和画图的函数文件放到一起。可又总是出错。 2012-04-09 08:02:52
xiegangmai: 金币+1, 谢谢参与！ 2012-04-08 22:53:07
plxh: 金币+5, ★★★很有帮助, 大哥你说得对，我就是赶鸭子上架，好高骛远了。昨天又学习了下function文件和M文件的不同。现在已基本实现了把三个1/4圆画到一张图上了。现在的问题就是我想把产生数据的m文件的内容和画图的函数文件放到一起。可又总是出错。 2012-04-09 08:02:52
|
针对你这个问题的解决方法是把 function plotSquare2(x0,y0,theta,a,mycolor) 以下的部分剪切出来，另存在一个m文件里就好了。 btw: 关于matlab的基础还是要学的。 |
6楼2012-04-08 21:22:32
7楼2012-04-09 08:05:27
8楼2012-04-09 09:14:52
匿名
用户注销 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1963.9
- 散金: 3479
- 红花: 4
- 帖子: 820
- 在线: 2782.4小时
- 虫号: 0
- 注册: 2013-01-04
- 性别: GG
- 专业: 理论和计算化学
9楼2014-12-02 15:55:44
