24小时热门版块排行榜

返回列表

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

plxh

木虫 (正式写手)

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

[求助] Matlab绘图再求助

在帖子http://muchong.com/bbs/viewthread.php?tid=4354741中部分实现了绘图目标。感谢sogooda朋友！

利用sogooda同学给出的代码实现了画出如下效果的图。
现在有新的想法，用什么方式可以看到小方块一个一个地方的变动过程呢？

用Matlab的GUI或者Real Time Workshop 或者S-Function？
本人新手，请各位同学不吝赐教！！
PS：mat文件里的数据是用把simulink里的数据用To Workspace导出到工作空间的。
效果图

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : matlab.mat

2012-04-09 09:27:44, 9.33 K

1楼2012-04-09 09:35:29

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
仿真EPI: 6
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
plxh: 金币+5, ★★★很有帮助, 这位兄弟谢谢你！你可能没看到。我要把三个数据回到一张图上，不过你这个很有帮助。 2012-04-10 10:17:04

CODE:

function aa

clear;clc;close all

load matlab3

for i=1:3;

    hold on;

    if i==1   %必须是用两个等号。

          figure1=figure('color','w')  

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

        x=eta12(:,2);

        y=eta12(:,1);

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

        theta=pi/2-eta12(:,3);

        myr=linspace(0,0,length(x));

        myg=linspace(0,1,length(x));

        myb=linspace(0,0,length(x));

        for jj=1:length(x)   

           plot(x,y);hold on;

            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);

            axis equal 

            pause(0.1)

    hold off;

        end

        colormap([myr', myg', myb'])

        colorbar

    elseif i==2

         figure1=figure('color','w')

        eta22=eta2(1:5:end,:);

        x=eta22(:,2);

        y=eta22(:,1);

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

        theta=pi/2-eta22(:,3);

        myr=linspace(0,0,length(x));

        myg=linspace(0,1,length(x));

        myb=linspace(0,0,length(x));

        for jj=1:length(x) 

            plot(x,y);hold on;

            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);

            axis equal

            pause(0.1)

    hold off;

        end

        colormap([myr', myg', myb'])

        colorbar

    elseif i==3        

          figure1=figure('color','w') 

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

        x=eta32(:,2);

        y=eta32(:,1);

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

        theta=pi/2-eta32(:,3);

        myr=linspace(0,0,length(x));

        myg=linspace(0,1,length(x));

        myb=linspace(0,0,length(x));

        for jj=1:length(x)   

            plot(x,y);hold on;

            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);

            axis equal

            pause(0.1)

    hold off;

        end

        colormap([myr', myg', myb'])

        colorbar

    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)