查看: 1897|回复: 3|关注: 0
运用min求列向量的最小值出现下标索引必须为正整数类型...
<h1 style="color:# 麦片财富积分
新手, 积分 6, 距离下一级还需 44 积分
a为一个43个元素的列向量
&&Columns 1 through 10
-221.36 -223.62 -226.47 -229.95 -231.82
&&Columns 11 through 20
-234.22 -236.50 -239.26 -242.26 -246.32
&&Columns 21 through 30
-249.72 -252.91 -256.84 -260.71 -263.85
&&Columns 31 through 40
-267.28 -271.52 -276.35 -280.99 -282.36
&&Columns 41 through 43
-276.90 -247.1121
max(a)=-221.4515
可是求min(a)总是出现
下标索引必须为正整数类型或逻辑类型。如何解决？
论坛优秀回答者
帖子最佳答案
关注者： 667
|此回复为最佳答案
说明你定义了一个名字为min的变量，所以当你写min(a)并不是调用了min函数，而是去访问了你的min变量。换一个变量名就可以了
<h1 style="color:# 麦片财富积分
说明你定义了一个名字为min的变量，所以当你写min(a)并不是调用了min函数，而是去访问了你的min变量。换一 ...
前面代码定义了一个min的向量。。。怪不得
<h1 style="color:# 麦片财富积分
对哦，我也是这样，哈哈
站长推荐 /2
机器视觉和人工智能在医疗设备中的应用及实现
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区查看: 2906|回复: 9
大神求教这个警告是什么原因造成的，如何修改程序？
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
才可以下载或查看，没有帐号？
大神求教这个警告是什么原因造成的，如何修改程序？
用ode45求解系统方程
[t ,Y]=ode45(@weifenfangcheng,[0:(2*pi/200):(288*2*pi),[0.01 0.001 0.01 0.001 0.01 0.001 0.01 0.001])；
当j=2时出现警告：“在t=1.处失败。在时间t处，若不将步长降至允许的最小值（2.）一下，积分公差要求无法满足。”
in ode45 at 308
下标索引必须为正整数类型或逻辑类型。
另外，我是在求解方程里面加了一项函数，在不添加前方程可以运行得到结果，由于多考虑其他因素添加了一项后就出现了错误。请问出现这种错误的原因是什么啊，怎么解决啊？
有没有知道，指点一下
好奇问下, 报错讯息是中文?
求助完整格式：出错代码和出错提示
好奇问下, 报错讯息是中文?
用的是最新的中文版matlab2014b，报错提示是中文的。
本帖最后由 牛小贱 于
09:34 编辑
好奇问下, 报错讯息是中文?
程序如下，主要问题就是我在程序里面未添加Fa这一项之前是可以进行计算得到图，添加Fa之后就出现警告，终止运行。
function fangchengfenchatu
tic
N=256;
M=200;
RelTol = 1e-6;& && && && && & %相对误差Relative tolerance
AbsTol = 1e-6;& && && && && & %绝对误差Absolute tolerance
options = odeset('RelTol',RelTol,'AbsTol',ones(1,8)*AbsTol);&&%需要改变的参数‘ones(1,4)’
ommiga=200:10:2500;
L=length(ommiga);
for j=1:L
j
&&[t,Y] = ode45( @fangcheng001,[0:(2*pi/M):(288*2*pi)],[0.01 0.001 0.01 0.001 0.01 0.001 0.01 0.001],options,ommiga(j));
&&y1=Y(:,1);
y11=y1((end-N*M):end);
for i=1:N
&&
& &Point1(i,j)=y11(1+(i-1)*M);
end
&&
end
ommiga = ommiga';
plot(ommiga,Point1,'.k','MarkerSize',1)
toc复制代码function dy=fangcheng001(t,y,ommiga)&&
%ommiga=450;%转速
m1=4;
m2=32.1;
c1=1050;
c2=2100;
c=0.00011;
R=0.025;
L=0.012;
e=0.00005;
miu=0.018;
k=2.5*10^7;
s=miu*ommiga*R*L*(R/c)^2*(L/2/R)^2;
M=s/(m1*ommiga^2*c);
%求解Fa;
xi=1.2e-3;
Rt=0.5;
Rb=0.37;
B1=35/360*2*
B2=40/360*2*
xii=0.83;
V=30;
p0=11.8;
C=V^2*sin(B1)*p0*(cos(B1)+xii*B2);
A1=(Rt^2-Rb^2)^2*pi*C*Rt/(Rt^2-Rb^2+2*Rt*xi)^2;
A3=3*(Rt^2-Rb^2)^2*pi*C*Rt^3/(Rt^2-Rb^2+2*Rt*xi)^4;
E=sqrt(y(5)^2+y(7)^2)/
Fa=A1*E+A3^3*E^3;
%无量纲方程
dy=[ y(2);
& & -c1/(m1*ommiga)*y(2)-k/(m1*ommiga^2)*(y(1)-y(5))+fx(y(1),y(3),y(2),y(4))*M;
& & y(4);
& & -c1/(m1*ommiga)*y(4)-k/(m1*ommiga^2)*(y(3)-y(7))-9.8/(c*ommiga^2)+fy(y(1),y(3),y(2),y(4))*M;
& & y(6);
& & -c2/(m2*ommiga)*y(6)-2*k/(m2*ommiga^2)*(y(5)-y(1))+e*cos(t)/c+Fa*cos(t)/(c*m2*ommiga^2);
& & y(8);
& & -c2/(m2*ommiga)*y(8)-2*k/(m2*ommiga^2)*(y(7)-y(3))+e*sin(t)/c-9.8/(c*ommiga^2)+Fa*sin(t)/(c*m2*ommiga^2);];
end复制代码
程序如下，主要问题就是我在程序里面未添加Fa这一项之前是可以进行计算得到图，添加Fa之后就出现警告，终 ...
fx,fy表示什么
新人报到，向前辈致敬！
本帖最后由 牛小贱 于
09:34 编辑
fx,fy表示什么
忘记写了&&不好意思
%油膜力程序
function out=fx(y1,y3,y2,y4)
alpha=atan((y3+2*y2)/(y1-2*y4))-pi/2*sign((y3+2*y2)/(y1-2*y4))-pi/2*sign(y3+2*y2);
G=2*(pi/2+atan((y3*cos(alpha)-y1*sin(alpha))/(sqrt(1-y1^2-y3^2))))/(sqrt(1-y1^2-y3^2));
V=(2+(y3*cos(alpha)-y1*sin(alpha))*G)/(1-y1^2-y3^2);
S=(y1*cos(alpha)+y3*sin(alpha))/(1-(y1*cos(alpha)+y3*sin(alpha))^2);
K=sqrt((y1-2*y4)^2+(y3+2*y2)^2)/(1-y1^2-y3^2);
out=K*(3*y1*V-sin(alpha)*G-2*cos(alpha)*S);
end
function out=fy(y1,y3,y2,y4)
alpha=atan((y3+2*y2)/(y1-2*y4))-pi/2*sign((y3+2*y2)/(y1-2*y4))-pi/2*sign(y3+2*y2);
G=2*(pi/2+atan((y3*cos(alpha)-y1*sin(alpha))/(sqrt(1-y1^2-y3^2))))/(sqrt(1-y1^2-y3^2));
V=(2+(y3*cos(alpha)-y1*sin(alpha))*G)/(1-y1^2-y3^2);
S=(y1*cos(alpha)+y3*sin(alpha))/(1-(y1*cos(alpha)+y3*sin(alpha))^2);
K=sqrt((y1-2*y4)^2+(y3+2*y2)^2)/(1-y1^2-y3^2);
out=K*(3*y3*V+cos(alpha)*G-2*sin(alpha)*S);
end复制代码
y1=Y(:,1);&&
y11=y1((end-N*M):end);
y1只有一个值，所以y11=y1((end-N*M):end)相当于y11=y1((1-256*200):1)=y1(-51199:1)，
y1下标索引只能是正值或者逻辑数值
<p id="rate_397" onmouseover="showTip(this)" tip="回复具有一定的启发引导性&体能 + 2 点
" class="mtn mbn">
y1=Y(:,1);&&
y11=y1((end-N*M):end);
y1只有一个值，所以y11=y1((end-N*M):end)相当于y11=y1((1-256*200 ...
非常感谢您的指导，我没有添加Fa之前是可以得到结果的。还有就是提示“当j=2时出现警告：“在t=1.处失败。在时间t处，若不将步长降至允许的最小值（2.）一下，积分公差要求无法满足。”应该怎么处理啊？
Powered by正整数类型或逻辑类型_中华文本库
型数据(用%lf或%le) 用于输入短整型数据(可用%hd,%ho,%hx) 指定输入数据所占宽度(列数),域宽应为正整数 表示本输入项在读入后不赋值给相应的变量 逻辑运算...
(1)&出厂价&只能输入3位整数和2位小数(整数部分可以不足3位);(2)&单位&...基本操作 1.将&CD.MDB&数据库中的&类型&及&出版单位&表导入到当前数据库中...
丌同,占用空间也丌同 数值型数据 (主要数据类型) 日期型 字节型 货币型 逻辑型 字符串型 对象型 变体型 一、 数值数据类型 数值类型分为整数型和实数型两大...
[j+1]对换; 其中 n 为正整数,则最后一行的语句频度在最坏情况下是( ) A...4.抽象数据类型的定义仅取决于它的一组__逻辑特性 _,而与_ 在计算机内部如何...
Matlab 常见错误及解决办法归纳 1.Subscript indices must either be real positive integers or logicals 中文解释:下标索引必须是正整数类型或者逻辑类型 出错原因:...
五、习语言函数分类 习语言函数可分为:系统函数、文件操作函数、字符串处理 ...(); 21、整数类型 输入整数( 无参数 ); 22、双精度类型 输入浮点数( 无...
逗号运算符等的基本运用和规则 语言的关系运算符和逻辑运算符, 2、学会使用 C...因此实型变量中 允许存放整型数 程序中,无论是整数还是实数, C)在 C 程序中...
SQL SERVER 常用数据类型 一、 整数数据类型 整数数据类型是最常用的数据类型之...四、逻辑数据类型 BIT: BIT 数据类型占用 1 个字节的存储空间,其值为 0 或...
(1)&出厂价&只能输入3位整数和2位小数(整数部分可以不足3位,小数部分也 ...一、 基本操作 1.在&库存&表中增加&单位&字段,其数据类型为&文本&,字段大小...
计算机中可以包含多个CPU, 以实现多个操作的并行处...两种工作模式均要求计算机网络的拓扑结构必须为总线型...运算器可以有多个,如整数运算器和浮点运算器等 45...}