24小时热门版块排行榜    

查看: 1989  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

plxh

木虫 (正式写手)

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

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

plxh

木虫 (正式写手)

引用回帖:
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);
...

兄弟,我看了你的新浪博客,你真是牛人!
3楼2012-04-08 12:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

sogooda

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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)



2楼2012-04-08 12:02:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plxh

木虫 (正式写手)

引用回帖:
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);
...

大哥,如果要一张图上画出来原理完全一样但半径不一样的三个呢?

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 3个.bmp
  • 2012-04-08 13:19:06, 2.52 M
  • 附件 2 : matlab3.mat
  • 2012-04-08 13:19:06, 30.89 K
4楼2012-04-08 13:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plxh

木虫 (正式写手)

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +9 yangfz 2026-03-17 9/450 2026-03-18 12:38 by 尽舜尧1
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 307求调剂 +3 冷笙123 2026-03-17 3/150 2026-03-18 09:55 by macy2011
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见