是否可以通过在matlab中调用c/c++代码来加速matlab绘图？

一般来说很容易调用mexMatlab 中的文件（用 c/c++ 编写）以加速某些计算。然而，根据我的经验，Matlab 的真正瓶颈是数据绘图。创建句柄非常昂贵，即使您只更新句柄数据（例如 XData、YData、ZData），这也可能需要很长时间。更糟糕的是，由于 Matlab 是单线程程序，因此不可能同时更新多个绘图。

因此我的问题: 是否可以编写一个 Matlab GUI 并调用 C++（或其他一些可并行代码）来处理绘图/可视化？我正在寻找一个可以在 Windows、Mac 和 Linux 上运行的跨平台解决方案，但是任何能让我在任一操作系统上开始使用的解决方案都非常感谢！

我找到了一个C++ 库 http://code.google.com/p/matplotpp/好像是用Matlab的plot()语法，但我不确定这是否会加快速度，因为我担心如果我绘制到 Matlab 中figure()窗口，事情可能会再次放缓。

我非常感谢以前处理过这种情况的人的任何评论和反馈！

EDIT:显然，我已经分析了我的代码，瓶颈是绘图（数十个带有大量数据的面板）。

EDIT2:为了让你获得赏金，我需要一个现实生活中的、最小的工作示例来说明如何做到这一点 - 暗示性答案对我没有帮助。

EDIT3:关于要绘制的数据：在最简单的情况下，考虑 20 个线图，每秒需要更新 1000000 个数据点。

EDIT4:我知道要绘制的点数量很大，但我从未说过这个问题很简单。我不能仅仅忽略某些数据点，因为在实际绘制数据之前无法评估哪些点是重要的（数据以亚毫秒时间分辨率采样）。事实上，我的数据是使用商业数据采集系统获取的，该系统带有数据查看器（用c++编写）。该程序可以毫无问题地可视化多达 60 个线图，甚至包含超过 1000000 个数据点。

EDIT5:我不喜欢当前讨论的方向。我知道对数据进行二次采样可能会加快速度 - 然而，这不是问题。这里的问题是如何让 c / c++ / python / java 接口与 matlab 一起使用，以便通过直接与硬件对话（或使用任何其他技巧/方式）来加速绘图

您是否尝试过将渲染方法更改为的简单解决方案OpenGL ?

opengl hardware;
set(gcf,'Renderer','OpenGL');

警告！ 在这种模式下会有一些东西消失，看起来会有点不同，但通常剧情会运行much更快，特别是如果您有硬件加速器。

顺便问一下，您确定您真的会获得性能提升吗？ 例如，根据我的经验，WPF图形在C#比 Matlab 慢得多，尤其是散点图和圆圈。

Edit：我想到了这样一个事实：实际绘制到屏幕上的点数不可能那么多。基本上，这意味着您需要在屏幕上有像素的地方进行插值。查看这个对象：

classdef InterpolatedPlot < handle
    properties(Access=private)
        hPlot;
    end

    methods(Access=public)
        function this = InterpolatedPlot(x,y,varargin)
            this.hPlot = plot(0,0,varargin{:});
            this.setXY(x,y);
        end        
    end    

    methods
        function setXY(this,x,y)
            parent = get(this.hPlot,'Parent');
            set(parent,'Units','Pixels')
            sz = get(parent,'Position');
            width = sz(3); %Actual width in pixels
            subSampleX = linspace(min(x(:)),max(x(:)),width);

            subSampleY = interp1(x,y,subSampleX);
            set(this.hPlot,'XData',subSampleX,'YData',subSampleY);
        end

    end
end

这是如何使用它的示例：

function TestALotOfPoints()
    x = rand(10000,1);
    y = rand(10000,1);

    ip = InterpolatedPlot(x,y,'color','r','LineWidth',2);

end

另一个可能的改进：另外，如果您的 x 数据已排序，您可以使用interp1q代替interp，这会快得多。

classdef InterpolatedPlot < handle
    properties(Access=private)
        hPlot;
    end

%     properties(Access=public)
%         XData;
%         YData;      
%     end

    methods(Access=public)
        function this = InterpolatedPlot(x,y,varargin)
            this.hPlot = plot(0,0,varargin{:});
            this.setXY(x,y);
%             this.XData = x;
%             this.YData = y;
        end        
    end    

    methods
        function setXY(this,x,y)
            parent = get(this.hPlot,'Parent');
            set(parent,'Units','Pixels')
            sz = get(parent,'Position');
            width = sz(3); %Actual width in pixels
            subSampleX = linspace(min(x(:)),max(x(:)),width);

            subSampleY = interp1q(x,y,transpose(subSampleX));
            set(this.hPlot,'XData',subSampleX,'YData',subSampleY);
        end

    end
end

以及用例：

function TestALotOfPoints()
    x = rand(10000,1);
    y = rand(10000,1);
    x = sort(x);
    ip = InterpolatedPlot(x,y,'color','r','LineWidth',2);

end