求matlab高手看看程序有什么问题?（新手,财富很少,呵呵）clear all;clc%输入初步确定的各连杆尺寸值l1=0.06;l2=0.6;l3=0.3;l7=0.4;xg=0.8;yg=0.55;w1=3.5*pi;if l7+l3+l2-l1>=(xg^2+yg^2)^0.5i=0; for t=0:0.01:pi; th1=w1*t;%求
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/19 14:41:31
求matlab高手看看程序有什么问题?（新手,财富很少,呵呵）clear all;clc%输入初步确定的各连杆尺寸值l1=0.06;l2=0.6;l3=0.3;l7=0.4;xg=0.8;yg=0.55;w1=3.5*pi;if l7+l3+l2-l1>=(xg^2+yg^2)^0.5i=0; for t=0:0.01:pi; th1=w1*t;%求
求matlab高手看看程序有什么问题?（新手,财富很少,呵呵）
clear all;
clc
%输入初步确定的各连杆尺寸值
l1=0.06;
l2=0.6;
l3=0.3;
l7=0.4;
xg=0.8;
yg=0.55;
w1=3.5*pi;
if l7+l3+l2-l1>=(xg^2+yg^2)^0.5
i=0;
for t=0:0.01:pi;
th1=w1*t;
%求解位置方程组（1）、（2）
f=@(x)([l2*cos(x(1))+0.325-l7*sin(x(2))+l1*cos(th1)-xg;
l2*sin(x(1))+0.233+l7*cos(x(2))+l1*sin(th1)-yg]);
juzhen=fsolve(f,[1,1]);
th2=juzhen(1);th4=juzhen(2);
%求杆7的角速度
w2=(-l1*w1*cos(th1-th4))/(l1*cos(th2-th4));
w4=(l1*sin(th1)*w1+l2*sin(th2)*w2)/(-l7*cos(th4));
w24=[w2 w4];
%求杆7的角加速度
a2=(l1*(w1^2)*sin(th1-th4)+l7*(w4^2)-l2*(w2^2)*sin(th4-th2))/(l2*cos(th2-th4));
a4=(l1*w1^2*cos(th1)+l2*a2*sin(th2)+l2*w2^2*cos(th2)-l7*w4^2*sin(th4))/(-l7*cos(th4));
a24=[a2 a4];
%求筛体的加速度
a=l7*a4;
%把角位移、角速度、角加速度和筛体加速度放在一个数组中；
i=i+1;
TH(i,:)=juzhen;
W(i,:)=w24;
A(i,:)=a24;
AD(i,:)=a;
end
%绘制筛体整体加速度曲线图
t=0:0.01:1.5*pi;
th1=w1*t;
figure(1)
plot(t,W(:,2));
xlabel('t/s');
ylabel('杆7角速度');
title('杆7的角速度')
grid on;
figure(2)
plot(t,A(:,2));
xlabel('t/s');
ylabel('杆7角加速度')
title('杆7的角加速度')
grid on;
figure(3)
plot(th1/w1,AD(:,1))
xlabel('t/s')
ylabel('a5（m/s2)')
title('筛体的加速度')
grid on;
else
disp('系统存在死点,请重新输入正确的杆长')
end
运行后显示：
Error using ==> plot
Vectors must be the same lengths.（向量必须是相同的长度）
Error in ==> jiexifa at 45
plot(t,W(:,2));
感觉问题出在把数据写入数组那里,但是不知道为什么错,怎么改?
求matlab高手看看程序有什么问题?（新手,财富很少,呵呵）clear all;clc%输入初步确定的各连杆尺寸值l1=0.06;l2=0.6;l3=0.3;l7=0.4;xg=0.8;yg=0.55;w1=3.5*pi;if l7+l3+l2-l1>=(xg^2+yg^2)^0.5i=0; for t=0:0.01:pi; th1=w1*t;%求
whos W t
Name Size Bytes Class Attributes
W 315x2 5040 double
t 1x472 3776 double
显然长度不同
t=0:0.01:1.5*pi;这句是多余的