24小时热门版块排行榜    

查看: 4196  |  回复: 2

bettyzhy1

铁虫 (小有名气)

[求助] 错误使用 == 矩阵维度必须一致。

能否帮忙看看这个程序怎么回事,求助感谢!!!
命令窗口显示以下部分内容:
max=5.7899,5.7164,5.5852,5.9084,5.5288,5.6218,5.6064,5.6153,5.5596,错误使用  ==
矩阵维度必须一致。

出错 FBG_Tmax (line 35)
    p1=find(y==max(y(t>(i-1)&t<i)));
我的程序如下:(就是绘制了一个正弦波形处来,需要找出波峰和波谷点求出其平均值,然后利用这两个值再解一个二元一次方程组,问题应该在于第二部分求最大值位置)
y=xlsread('E:\Postgraduate\摩擦自激振动实验\实验数据\2017年8月17日\2017_8_17已转\52.xlsx');
y0=xlsread('E:\Postgraduate\摩擦自激振动实验\实验数据\2017年8月17日\2017_8_17已转\52_s.xlsx');

Fs = 4000;                    % Sampling frequency
T = 1/Fs;                     % Sample time
L = length(y);                % Length of signal
t = (0:L-1)*T;

y1=y(:,9);
y2=y(:,10);
y01=y0(:,9);
y02=y0(:,10);
y=11.036999*(y1-(mean(y01)))-(y2-(mean(y02)));
figure
plot(t,y);
xlabel('Vibration Time [s]','Fontsize',12);
ylabel('Tension of Spring [N]','Fontsize',12);
legend('Tension','Location','southeast');
%%
fprintf('max=');
a=0;
for i=1:10;
    p1=find(y==max(y(t>(i-1)&t<i)));
    text(t(p1),y(p1),['(',num2str(t(p1)),',',num2str(y(p1)),')'],'color','r');%在图中将最大值坐标表示出来
    fprintf([num2str(y(p1)),',']);%输出最大值
    a=y(p1)+a;
end
m=a/10;%最大值平均值
fprintf('\nm=%f\nmin=',m);
%%
b=0;
for j=1:10;
    p2=find(y==min(y(t>(j-1)&t<j)));
    text(t(p2),y(p2),['(',num2str(t(p2)),',',num2str(y(p2)),')'],'color','r');
    fprintf([num2str(y(p2)),',']);
    b=y(p2)+b;
end
n=b/10;
fprintf('\nn=%f\n',n);
%%

A=[1,-2;1,2];
b=[2*m;-2*n];
x=A\b;
fprintf('x=%f\n',x)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bettyzhy1

铁虫 (小有名气)

有时候还会出现错误使用 horzcat
串联的矩阵的维度不一致。
2楼2017-08-21 11:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
月只蓝: 金币-200, 违规存档, 多次恶意灌水。 2017-08-28 09:57:24
月只蓝: 金币-200, 多次恶意灌水。 2017-08-28 09:58:21
月只蓝: 金币-200, 多次恶意灌水。 2017-08-28 09:58:45
月只蓝: 金币-200, 多次恶意灌水。 2017-08-28 09:59:06
3楼2017-08-22 11:04:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bettyzhy1 的主题更新
信息提示
请填处理意见