24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

plxh

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3257.5
散金: 6
帖子: 527
在线: 92.2小时
虫号: 1286569
注册: 2011-05-05
专业: 控制理论与方法

[求助] 继续求助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

1楼2012-04-08 09:35:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

plxh

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3257.5
散金: 6
帖子: 527
在线: 92.2小时
虫号: 1286569
注册: 2011-05-05
专业: 控制理论与方法

xiegangmai: , 帖子中代码部分可以用“插入代码”功能 2012-04-08 22:52:58

引用回帖:

2楼: Originally posted by sogooda at 2012-04-08 12:02:00:
兄弟你的数据太密了，画出来都看不出来效果。
下面的程序里是将数据1/5稀疏之后的效果图。

function Untitled2
clear;clc;close all
load matlab
eta1=eta1(1:5:end,

;

x=eta1(:,2);
y=eta1(:,1);
...

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

sogooda

金虫 (正式写手)

应助: 19 (小学生)
金币: 1497.4
红花: 2
帖子: 737
在线: 94.6小时
虫号: 438143
注册: 2007-09-16
专业: 防灾工程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
plxh: 金币+5, ★★★很有帮助, 先谢谢！我理解理解！ 2012-04-08 12:03:51
xiegangmai: 金币+2, 谢谢应助！ 2012-04-08 22:52:38

兄弟你的数据太密了，画出来都看不出来效果。
下面的程序里是将数据1/5稀疏之后的效果图。

CODE:

function Untitled2

clear;clc;close all

load matlab

eta1=eta1(1:5:end,:);

x=eta1(:,2);

y=eta1(:,1);

% theta=eta1(:,3); % ?

theta=pi/2-eta1(:,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

saveas(gcf,'untitled2.jpg')

function plotSquare(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)