24小时热门版块排行榜    

查看: 2374  |  回复: 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的回帖

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);
...

兄弟,我看了你的新浪博客,你真是牛人!
3楼2012-04-08 12:04:33
已阅   回复此楼   关注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的回帖

sogooda

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
xiegangmai: 金币+1, 谢谢参与! 2012-04-08 22:53:07
plxh: 金币+5, ★★★很有帮助, 大哥你说得对,我就是赶鸭子上架,好高骛远了。昨天又学习了下function文件和M文件的不同。现在已基本实现了把三个1/4圆画到一张图上了。现在的问题就是我想把产生数据的m文件的内容和画图的函数文件放到一起。可又总是出错。 2012-04-09 08:02:52
引用回帖:
5楼: Originally posted by plxh at 2012-04-08 13:44:10:
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=eta ...

针对你这个问题的解决方法是把
function plotSquare2(x0,y0,theta,a,mycolor)
以下的部分剪切出来,另存在一个m文件里就好了。

btw: 关于matlab的基础还是要学的。
6楼2012-04-08 21:22:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plxh

木虫 (正式写手)

引用回帖:
6楼: Originally posted by sogooda at 2012-04-08 21:22:32:
针对你这个问题的解决方法是把
function plotSquare2(x0,y0,theta,a,mycolor)
以下的部分剪切出来,另存在一个m文件里就好了。

btw: 关于matlab的基础还是要学的。

看了你的回复,我受到了启发,是不是我可以在产生数据的m文件里面引用画图函数?
7楼2012-04-09 08:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plxh

木虫 (正式写手)

引用回帖:
6楼: Originally posted by sogooda at 2012-04-08 21:22:32:
针对你这个问题的解决方法是把
function plotSquare2(x0,y0,theta,a,mycolor)
以下的部分剪切出来,另存在一个m文件里就好了。

btw: 关于matlab的基础还是要学的。

谢谢,已经实现了一个m文件即取得数据又绘出图。

向你请教一下,我是否可以用什么方式可以看到小方块  一个一个地方的变动过程呢?

Matlab的GUI或者Real Time Workshop 或者S-Function?

8楼2012-04-09 09:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

本帖仅楼主可见
9楼2014-12-02 15:55:44
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 plxh 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 15/750 2026-08-20 17:28 by wuchongjun
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见