主页 / coding / 问题

问题[matlab](coding)

Martin Hope
zajer
Asked: 2024-11-06 02:37:45 +0800 CST
  • 6

我有两个矩阵：

A = [ -1 0  0;
       1 1 -1;
       0 -1 1 ];
B = [-1; 0; 1];

我想解下列方程：

Ax=B

当我使用mldivide函数时，我得到一个 NaN 矩阵

X = mldivide(A,B)
X =

   NaN
   NaN
   NaN

知道该问题有多种解决方案，我手动测试了其中一种，即[1;0;1]返回B：

A*[1; 0; 1]

并且，正如预期的那样，我确信它是多个解决方案之一。所以我的问题是：为什么mldivide返回错误的解决方案？

matlab
Martin Hope
user28134562
Asked: 2024-11-05 15:25:50 +0800 CST
  • 5

我正在尝试通过将一些初始化分离到单独的 .m 文件中来整理我的 MATLAB 应用程序代码。为此，我为每种类型的组件设置了各种文件（例如按钮、图形等的文件）。我正尝试从按钮文件访问主初始化文件中的函数。我的代码在按钮 .m 文件中如下所示：

classdef buttons < handle
    methods
        %initializes the UI
        function buttonCreate(app)
           
            %Create file load 1
            app.fileload1 = uibutton(app.gridLayout, 'push');
            app.fileload1.FontSize = 36;
            app.fileload1.Layout.Row = [8 9];
            app.fileload1.Layout.Column = 1;
            app.fileload1.Text = 'Load 1';

            %proceeds to create the rest of the buttons
        end
    end
end

现在我尝试访问buttonCreate()主初始化文件中的函数initialize.m：

classdef initialize < handle
    properties 
        fig
        gridLayout
        axes

        fileload1
    end
    methods
        %initializes the UI
        function init(app)
            %create canvas
            import buttons.*;
            fig = uifigure;
            fig.Position = [100 100 1920 1080];
            movegui(fig,'center');
            fig.Name = "Audio Editor";

            %Create grid layout
            gridLayout = uigridlayout(fig);
            gridLayout.ColumnWidth = {'1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x'};
            gridLayout.RowHeight = {'1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x', '1x'};

            buttonCreate(app);

        end
    end
    %code for calling and deleting
    methods
        %calls code to create canvas upon app start
        function app = initialize
            init(app)
        end
        %removes the app and deletes app.fig
        function delete(app)
            delete(app.fig);
        end
    end
end

Error in initialize/init (line 41)
            buttonCreate(app);
            ^^^^^^^^^^^^^^^^^^
Error in initialize (line 54)
            init(app)
            ^^^^^^^^^

这导致程序UIFigure仍在创建，但没有按钮，并且终端出现上述错误。

matlab
Martin Hope
Andrew Geistfeld
Asked: 2024-10-29 01:31:15 +0800 CST
  • 6

我正在模拟一个 APM（自动电源管理）系统，其中机器调整其齿轮和 RPM 以对可变工作负载做出反应。我以拖拉机为模型。因此，对于此模型，可变工作负载将增加/减少拖拉机必须拉动的阻力。

我被变速箱难住了。我有确定所需车轮扭矩的代码，我需要进行逆表查找。表格设置为 f(x,y)=z，其中 x=油门位置/转速，y=档位，z=车轮扭矩。在这种情况下，我知道我的油门位置，也知道所需的扭矩，我需要一种方法让它给我相应的档位。

该表如下所示：

桌子

我知道如何使用二维表查找从 x 和 y 中找到 z，但不知道如何使用 x 和 z 来查找 y。

matlab
Martin Hope
qtqt
Asked: 2024-10-19 13:49:21 +0800 CST
  • 5

我使用的是 Matlab 符号。我有一个大小为 (Nx3) 的矩阵 A，其中 A 的每个条目都是正整数。我想创建一个维度为 (max(max(A))x2) 的矩阵 Q。Q 的第 i 行对应于 A 的列索引，因此 A(k,r) = i。因此，如果 A(k_1,r_1) = i 且 A(k_2,r_2) = i，我们可以设置 Q(i,:)=[k_1,k_2]。对于 A 的唯一条目，Q 的相关行的每一列都相等。

有没有办法将任务矢量化？我可以使用循环执行任务，但在 MATLAB 中效率很低（当矩阵 A 有很大的行并且 max(max(A)) 也很大时）。这是一个例子。

A=[  1    10     9
     1     2     4
     3     2    11
     3    12     5
     6    13     4
     6     7    15
     8     7     5
     8    14    16]

在上面的例子中，预期结果是

Q =
     1     1
     2     2
     1     1
     3     3
     3     3
     1     1
     2     2
     1     1
     3     3
     2     2
     3     3
     2     2
     2     2
     2     2
     3     3
     3     3

这就是我的慢速代码所做的

max_Q = max(max(A));
Q     = zeros(max_Q,2);
for i = 1 : max_Q
    [~,cols] = find(i == A); 
    Q(i,:) = cols;
end
matlab
Martin Hope
Varga
Asked: 2024-10-19 01:19:20 +0800 CST
  • 5

有一个任务：模拟海面上的波浪（3D 情况）。我遇到的问题是，当模拟大片海洋（1 公里 x 1 公里）的波浪时，程序需要很长时间才能执行。事实上它根本没有执行。因为在等待结果的半个小时里，它从未被计算过。

我的代码示例：

clear,clc
%% Data 
x = 0:1:100; % coordinates х [m]
y = 0:1:100; % coordinates y [m]
g = 9.81; % gravitational constant [m/s^2]
speed = 5; % wind velocity [m/s]
w0 = (g / speed); % norm.frequency [Hz]
dw = 0.1; % frequency step [rad/s]
w = 0.8:dw:11.1; % frequency [rad/s] 
dtt = pi / 18; % angular step [rad]
theta = 0:dtt:pi; % direction angles, angles between the wavevector & coordintae axis [rad]

%% P-M spectrum, Frequency-Angular spectrum & Amplitude
Psi = 8.1e-3 .* ((w/w0).^(-5)) .* exp((-0.74) ./ ((w/w0).^(4))); % P-M spectrum [none]
Phi = ((speed)^(5)/g^(3)) * Psi; % self-similar spectrum [s*m^2]
Sw = Phi / 2; % frequency spectrum [s*m^2]

St = cos(theta).^(4); % angular spectrum [none]
Norm = trapz(dtt, St); % norm.coefficient [none]
Swt = Sw .* St'; % frequency-angular spectrum [s*m^2]

eta0 = sqrt((Swt * dw * dtt) ./ Norm); % amplitude [m]

figure(1);
subplot(2,1,1)
plot(w, Psi);
title('$$\Psi$$($$\omega$$) - P-M spectrum', 'Interpreter', 'LaTex');
xlabel('\omega [rad/s]');
ylabel('\Psi [none]');
grid on;
subplot(2,1,2)
plot(w, Swt); 
title('$$S(\omega , \theta)$$($$\omega$$) - frequency-angular spectrum', 'Interpreter', 'LaTex');
xlabel('\omega [rad/s]');
ylabel('S(\omega,\theta) [s*m^2]');
grid on;

%% Setting the initial phase parameter
phase = 2*pi*rand(length(theta),length(w)); %% random initial phase ranging from 0 to 2pi [rad]

%% Surface Waving [Linear, 3D (eta & x,y)] at different harmonics & random phase (at one moment in time), different directions of the wavevector (multiple angles)
t = 0; % time moment [s]
Kabs = (w.^2) / g; % wavevector module [rad/m]
Kx = Kabs .* cos(theta)'; % projection of the wavevector onto the x-axis [rad/m]
Ky = Kabs .* sin(theta)'; % projection of the wavevector onto the y-axis [rad/m]

eta = zeros(length(x),length(y),length(theta),length(w)); % reserving space for calculation results
tic
for i = 1:length(x)
    for j = 1:length(y)
        eta(i,j,:,:) = eta0 .* cos(w * t - Kx .* i - Ky .* j + phase);
    end
end
toc
% sum(sum(eta,4),3) - double sum of eta over all harmonics (frequencies) and wavevector directions (angles theta),
% where '4' и '3' summation indicator for variable frequency and angle
etaW = sum(eta,4);
etaWA = sum(etaW,3);

figure(2)
surf(x,y,etaWA);
title('\eta(x,y) - surface waving');
xlabel('x [m]');
ylabel('y [m]');
zlabel('\eta [m]');
cbar = colorbar;
cbar.Label.String = '\eta [m]';
grid on
shading flat

我能够使用的代码优化方法之一是创建一个“空”的 4D 数组（零的 4D 数组）eta = zeros(length(x),length(y),length(theta),length(w));，在执行循环后，计算结果将填充到该数组中：

eta = zeros(length(x),length(y),length(theta),length(w)); % reserving space for calculation results
tic
for i = 1:length(x)
    for j = 1:length(y)
        eta(i,j,:,:) = eta0 .* cos(w * t - Kx .* i - Ky .* j + phase);
    end
end
toc

然后我按频率和角度变量总结结果：

etaW = sum(eta,4);
etaWA = sum(etaW,3);

这样就提前为结果准备好了地方。这有点帮助。例如，x = 0:1:100; y = 0:1:100;使用此方法对一个区域（100 mx 100 m）的代码执行时间为 0.7 秒（不使用则为 3.9 秒）。对于一个区域x = 0:1:500; y = 0:1:500;（500 mx 500 m），执行时间约为 19 秒（不使用则...我不知道，我没有等待代码执行，但结果发现它很长）。但是，对于一个区域x = 0:1:1000; y = 0:1:1000;（1000 mx 1000 m），我很长时间都没有得到想要的结果（感觉根本得不到）。

在我的情况下，还有其他方法可以实现所需的结果并优化我的代码，以便它可以应对如此规模的数据（同时，无需改变数组中的步骤）？

注意：我的电脑有 16 GB 的 RAM。我的第二台电脑（我最初用来进行计算的那台电脑）在执行程序时完全挂断了，所以我不得不“手动”重启它。所以我猜它的 RAM 更少了。

matlab
Martin Hope
user2587726
Asked: 2024-09-20 20:33:09 +0800 CST
  • 5

我正在进行一些线性建模并绘制输出。我希望在标题 title doc中包含 r^2 值，但以 2 作为上标。

由于某种原因，在绘制线性模型输出fitlm doc时，文档中的代码行未按预期生成标题。如何解决这个问题？

clear
clc

%Make test data
%----------------------------
x1 = 1:0.01:3;
r = -1 + (1+1)*rand(1,201);
x1 = x1 + r;

y1 = 1:0.01:3;
r = -1 + (1+1)*rand(1,201);
y1 = y1 + r;

x1 = x1';
y1 = y1';
label1 = zeros(length(y1),1);

Tclust1 = [label1,x1,y1];
%----------------------------
x2 = 4:0.01:6;
r = -1 + (1+1)*rand(1,201);
x2 = x2 + r;

y2 = 10:0.01:12;
r = -1 + (1+1)*rand(1,201);
y2 = y2 + r;

x2 = x2';
y2 = y2';
label2 = ones(length(y2),1);

Tclust2 = [label2,x2,y2];

%----------------------------
x3 = 5:0.01:7;
r = -1 + (1+1)*rand(1,201);
x3 = x3 + r;

y3 = 11:0.01:13;
r = -1 + (1+1)*rand(1,201);
y3 = y3 + r;

x3 = x3';
y3 = y3';
label3 = label2+1;

Tclust3 = [label3,x3,y3];
%----------------------------
T = [Tclust1;Tclust2;Tclust3];

%Prep for LM -------------------------------------------------------------
%Classify by labels
Label_idx = zeros(length(T(:,1)),1);

for qq = 1:length(Label_idx)  
    if T(qq,1) > 0
        Label_idx(qq) = 1;
    end
end

for qq = 1:length(Label_idx)
    if T(qq,1) > 1
        Label_idx(qq) = 2;
    end
end

%dummy variables
for qq = 1:length(Label_idx)
    if Label_idx(qq) == 0
        Label_1(qq) = 1;
    else
    Label_1(qq) = 0;
    end
end

Label_1 = Label_1';

for qq = 1:length(Label_idx)
    if Label_idx(qq) == 2
        Label_2(qq) = 1;
    else
    Label_2(qq) = 0;
    end
end

Label_2 = Label_2';

%Do the Linear modelling -------------------------------------------------

%Input predictor
pred = (T(:,2));
%Input comparator
comp = (T(:,3));

%interaction parameters
int_lo = pred.*Label_1;
int_hi = pred.*Label_2;

%predictors
    X = [pred,Label_1,Label_2,int_lo,int_hi];
    LR_mdl = fitlm(X,comp);

%plot ==========================================================
plot(LR_mdl)
arrayfun( @(x) set(x,'DisplayName',strrep(x.DisplayName,'*',' × ')), findall(gcf,'type','line') )
% title('title^2')
% title('^{super} normal _{sub}')
title('text ^{superscript} ')

figure(2)
scatter(x1,y1)
title('text ^{superscript} ')
`

请注意，当仅绘制部分原始数据（而非线性模型）时，代码行确实会为图 2 生成所需的上标。

matlab
Martin Hope
Davide Rigato
Asked: 2024-09-20 15:21:34 +0800 CST
  • 6

我正在使用 MATLAB。我有一张图形边缘的黑白图像。我尝试使用函数“find”按顺序查找边缘的坐标（1）。尽管该函数不会按顺序查找 1，但它只会逐行检查是否有 1。有没有办法按顺序找到 1 的序列（如参数化）？谢谢。

这里您可以找到边缘的图像。我必须按顺序找到 1。

边缘

我只是尝试使用 find 函数：

[row_edge, col_edge] = find(frame_edge==1);

matlab
Martin Hope
Fabi
Asked: 2024-09-12 21:06:01 +0800 CST
  • 6

我在 simulink 中有一个模型，它控制一个实时应用程序。我有一个预定义的测试计划，它基本上是一个表格，使用状态流图逐行处理。在某些步骤中，记录非常小的时间步骤非常重要，以查看应用程序如何处理输入的详细信息。对于测试计划中的其他步骤，记录这么多步骤并不重要，因为我必须再次调节我的测试装置或类似的东西。“可变采样率”可以是模型基本采样率的倍数。有没有办法用每行不同的采样时间记录这些数据？测试计划可能如下图所示。命令中的每一行都有匹配的采样时间。在我的状态流图中，我使用每一行中的值将某些内容发送到我的测试装置。在反馈循环中，如果满足要求，它将转到测试计划中的下一行。即达到目标值。

测试计划：

测试计划

Stateflow图流程： Stateflow 图表流程

我在测试计划本身中定义了每行的采样时间，并通过处理整个测试计划的 staeflow 图表给出了值。我尝试使用 matlab 函数块将变量写入基础工作区。然后我在工作区块中将此变量用作采样时间。奇怪的是，它改变了第一条处理线的工作区块的采样时间，但在整个模拟过程中保持这种状态。我也尝试对速率转换块进行同样的操作。

matlab
Martin Hope
Lucy Yang
Asked: 2024-08-29 22:04:24 +0800 CST
  • 5

从下图可以看出以下两个问题：

  1. 图例字体大小不一致：虽然图例 1 和图例 2 都设置为字体大小 6，但它们似乎以不同的大小显示。
  2. 图例标记过大：调整 legend2.ItemTokenSize 对图例中的标记大小影响不大，标记仍然显得太大。

如何解决上述问题？如能提供帮助，我将不胜感激。非常感谢！

图1

这是的代码legend1：

    legend1=legend(legend_entries, ...
                          'FontName', 'Times New Roman', 'FontSize', 6,...
                          'Location', 'NorthEastoutside',...
                          'Box', 'on', 'NumColumns', 1 );
    
    set(legend1, 'Position', [0.84, 0.41, 0.1, 0.2]); 

    
    legend1.ItemTokenSize = [10, 10];

    
    legend('AutoUpdate', 'off');

这是的代码legend2：

if nargin > 1 && ~isempty(set_legend2)
        
        new_ax = axes('Position', get(gca, 'Position'),...
                      'Units', get(gca, 'Units'), 'Visible', 'off');
        legend2=legend(new_ax,set_legend2,plot_info.leg,...
                              'FontName', 'Times New Roman', 'FontSize', 6, ...    
                              'Location', 'northwest', ...
                              'Box', 'off', 'NumColumns', 2 );
   
        set(legend2, 'Position', [0.55, 0.75, 0.1, 0.1]); 
         
        % 调整图例大小和位置
        legend2.ItemTokenSize = [10, 10]; 
    end
matlab
Martin Hope
Annika
Asked: 2024-08-20 22:20:32 +0800 CST
  • 6

我有以下问题：我从图像中提取了一个菱形并将其转换为二进制图像。我需要在该菱形内找到最大的垂直距离。

使用我当前的代码，我可以检测到最大距离（以红色显示），无论方向如何：

带有水平线的菱形图像

但我只对菱形上角点和下角点之间的最大距离感兴趣（所以想象红线旋转约 90 度），这是第二大距离。

我的代码如下（受此帖子启发：找到 Canny 边缘输出的最大欧几里得距离并在点之间画一条线）

clc; clear all;

im = imread('picture');
sharp = imsharpen(im);
grayImage1 = im2gray(sharp);
BinaryImage1 = imbinarize(grayImage1, 0.7);
imshow(BinaryImage1);
title('Extract rhombus');

[posX1, posY1, P1] = impixel;
Rhombus1 = bwselect(BinaryImage1, posX1, posY1);
BW1 = edge(Rhombus1,'Canny',0.8);

hold on;
[ y, x] = find( BW1);
points = [ x y];
[d,idx] = pdist2( points, points, 'euclidean', 'Largest', 1);
idx1 = idx( d==max(d));
[~,idx2] = find(d==max( d));
p={};

for i=1:length(idx1)
   p{end+1} = [ points(idx1(i),1), points(idx1(i),2), points(idx2(i),1), points(idx2(i),2)];
end

figure; 
imshow(BW1);
hold on;
for i=1:numel(p)
    line( [ p{i}(1), p{i}(3)], [p{i}(2), p{i}(4)]);
    hdl = get(gca,'Children');
    set( hdl(1),'LineWidth',2);
    set( hdl(1),'color',[1 0 0]);
end

hold off;

tmp = p{1,1};
dx = tmp(1,1) - tmp(1,3);
dy = tmp(1,2) - tmp(1,4);
Distance = sqrt((tmp(1,3)-tmp(1,1))^2+(tmp(1,4)-tmp(1,2))^2)

手动单击点在这里没有其他选择，我需要它自动化。

我非常感谢您的帮助！

matlab