🏆本文收录于「Bug调优」专栏,主要记录项目实战过程中的Bug之前因后果及提供真实有效的解决方案,希望能够助你一臂之力,帮你早日登顶实现财富自由🚀;同时,欢迎大家关注&&收藏&&订阅!持续更新中,up!up!up!!
问题描述
核心是解三自由度微分方程组,解出θ关于时间t的变化图像.设定一个功率恒定变加速度的工况,图中时间t是关于速度v的函数,取反函数v=@(t)finverse(t)代入 wm_func中求得角速度.图一是主程序,图二是参数和方程.系统提示我未定义与 ‘double’ 类型的输入参数相对应的函数 ‘finverse’,如图三所示,想问一下我这个错在哪里有没有修改办法或者建议
clear;clc;close all;func
%初始参数,该参数变化,myfun程序要对应调整
syms t;
syms v;
t=@(v)(505.2*ln(49897/(55000-398.7*v))+1460*(12.8-v)/398.7);
v=@(t)finverse(t);
wm_func = @(t)(55*(finverse(t)*3.6)*9550/(3600*0.4064/(0.93*14)))*pi*2/60;
zi=4;%内转子极对数
zo=14;%外转子极对数k
z1=32;%小齿轮齿数
z2=99;%大齿轮齿数
ms=4.62e-3;
I=z2/z1;%机械齿轮速比
G=(zo+zi)/zi;%磁齿轮速比
%一、计算过程
tmax=5;
tspan = 0:1e-4:tmax;
% tmax=900*2*pi/wm;% 计算时长,按输入轴的周期数计算
bc=2e-5;
y0=[0;0;0;wm_func(0);wm_func(0)/G;wm_func(0)/(G*I)];
% y0=[0;0;0;wm_func;wm_func/G;wm_func/(G*I)];
[t,y] = ode45(@(t,y)myfun3(t,y,wm_func), tspan, y0);
% [t,y] = ode45(@(t,y)myfun3(t,y), tspan, y0);
% [t,y]=odeRK4sys('myfun3',tmax,bc,y0);
figure;plot(t,y(:,1)/G-y(:,2));
title('时域响应图');
xlabel('时间 (s)');
ylabel('扭转振幅(rad)');%磁齿轮的相对弹性角度
figure;plot(t,y(:,2)/I-y(:,3));
title('时域响应图');
xlabel('时间 (s)');
ylabel('扭转振幅(rad)');%机械齿轮的相对弹性转角;
figure;plot(t,y(:,1),t,y(:,2),t,y(:,3));
title('时域响应图');
xlabel('时间 (s)');
ylabel('转速(rad/s)');
function dydt=myfun3(t,y,wm_func)
wm = wm_func(t);
zi=4;%内转子极对数
zo=14;%外转子极对数
z1=32;%小齿轮齿数
z2=99;%大齿轮齿数
arf0=20/180*pi;%压力角
ms=4.62e-3;
r2=ms*z1/2;%小齿轮分度圆
r3=ms*z2/2;%大齿轮分度圆
rb2=r2*cos(arf0);%小齿轮基圆
rb3=r3*cos(arf0);%大齿轮基圆
I=z2/z1;%机械齿轮速比
G=(zo+zi)/zi;%磁齿轮速比
% Tm =(1460*9.8*0.012+1.04*1460*2.67+0.32*2.25*(10+2.67*t*3.6)^2/21.15)*0.4064/(0.93*14);%起步加速
Tm=9550*55/(60*wm/(2*pi)); ;
TL=Tm*I*G;%变速器负载
%3 转动惯量
IM=0.081197;%电机转子转动惯量kgm2
I0=0.181197;%磁齿轮低速级转动惯量kgm2
I1=0.013663;%磁齿轮高速级转动惯量kgm2
I2=0.0230596;%机械小齿轮转动惯量
I3=1.9825692;%机械大齿轮转动惯量
IL=21;%整车等效转动惯量
%4 刚度与阻尼
%齿轮的阻尼和刚度
w1=wm*G;%机械小齿轮输入转速
kp=10.51e8;
kb=0.25e8;%两个齿轮的刚度曲线为正弦曲线,
km=kp+kb*sin(w1*t*z1);
sg=0.16;%齿轮阻尼比
cm=2*sg*sqrt(kp*I2*I3/(I2*rb3^2+I3*rb2^2));
b=2e-4;
xx=rb2*y(2)-rb3*y(3);
dxx=rb2*y(5)-rb3*y(6);
if xx-b>0
fx=xx-b;
elseif xx+b<0
fx=xx+b;
else
fx=0;
end
F23=km*fx+cm*dxx;
T0=1.7*Tm*G;
Tc1=0.2*T0*sin(wm*t);
Tc=0.1*(zi*y(4)-(zi+zo)*y(5));
T21=T0*sin(zi*y(1)-(zi+zo)*y(2))+Tc+Tc1;
y(7)=(Tm-T21/G)/(I0+IM);
y(8)=(T21-r2*F23)/(I2+I1);
y(9)=(-TL+r3*F23)/(I3+IL);
dydt=[y(4)
y(5)
y(6)
y(7)
y(8)
y(9)];
end
如上问题有来自我自身项目开发,有的收集网站,有的来自读者,如有侵权,立马删除。
解决方案
如下是上述问题的解决方案,仅供参考:
在 MATLAB 中,错误信息 “未定义与 ‘double’ 类型的输入参数相对应的函数 ‘finverse’” 通常表示 MATLAB 无法找到名为 finverse
的函数,该函数应该接受 ‘double’ 类型的输入参数。在您提供的代码中,问题出现在以下这行:
v=@(t)finverse(t);
这里,您尝试定义一个匿名函数 v
,它使用 finverse
作为其操作。然而,finverse
并没有在您的代码中定义,或者没有被正确引用。
可能的解决方案:
- 定义
finverse
函数:如果finverse
是您自定义的函数,确保它在 MATLAB 的搜索路径中,并且具有正确的定义。例如:
function tinv = finv(t)
% 这里是 finverse 函数的实现
tinv = ... % 根据 t 计算逆函数的输出
end
- 使用 MATLAB 内置函数:如果
finverse
应该是 MATLAB 的内置函数,确保函数名正确无误。例如,如果您想使用fzero
函数来找到函数的零点,应该这样写:
v = @(t) fzero(@(x) yourFunction(x) - t, initialGuess);
其中 yourFunction
是您要查找零点的函数,initialGuess
是一个初始猜测值。
-
检查变量作用域:如果您在其他地方定义了
finverse
,请确保它可以在当前作用域中访问到。 -
检查拼写和大小写:MATLAB 是区分大小写的,确保函数名的拼写和大小写都是正确的。
-
使用
feval
或eval
:如果finverse
是一个字符串形式的函数名,您可以使用feval
或eval
来调用:
v = @(t) feval('finverse', t);
或者
v = @(t) eval(['finverse(' num2str(t) ')']);
-
检查 MATLAB 版本:确保您使用的 MATLAB 版本支持您尝试调用的函数。
-
查看 MATLAB 函数索引:如果您不确定
finverse
的正确用法,可以查看 MATLAB 的函数索引或使用helpwin finv
(如果finverse
是指fmin
的逆函数,这可能是一个错误)。
特别注意:
- 您提供的代码中
t
是关于v
的函数,但您试图定义v
为t
的函数的反函数。这需要t(v)
是一个单调函数,以便存在唯一的反函数。 - 如果
t(v)
不是单调的,或者您没有正确定义finverse
,MATLAB 将无法计算v(t)
。
请根据上述建议检查和修改您的代码。如果 finverse
是您自定义的函数,请确保它在 MATLAB 的搜索路径中,并且具有正确的定义。如果 finverse
是内置函数或您打算使用特定的 MATLAB 函数,请确保正确引用它。
希望如上措施及解决方案能够帮到有需要的你。
PS:如若遇到采纳如下方案还是未解决的同学,希望不要抱怨&&急躁,毕竟影响因素众多,我写出来也是希望能够尽最大努力帮助到同类似问题的小伙伴,即把你未解决或者产生新Bug黏贴在评论区,我们大家一起来努力,一起帮你看看,可以不咯。
若有对当前Bug有与如下提供的方法不一致,有个不情之请,希望你能把你的新思路或新方法分享到评论区,一起学习,目的就是帮助更多所需要的同学,正所谓「赠人玫瑰,手留余香」。
☀️写在最后
ok,以上就是我这期的Bug修复内容啦,如果还想查找更多解决方案,你可以看看我专门收集Bug及提供解决方案的专栏「Bug调优」,都是实战中碰到的Bug,希望对你有所帮助。到此,咱们下期拜拜。
码字不易,如果这篇文章对你有所帮助,帮忙给bug菌来个一键三连(关注、点赞、收藏) ,您的支持就是我坚持写作分享知识点传播技术的最大动力。
同时也推荐大家关注我的硬核公众号:「猿圈奇妙屋」 ;以第一手学习bug菌的首发干货,不仅能学习更多技术硬货,还可白嫖最新BAT大厂面试真题、4000G Pdf技术书籍、万份简历/PPT模板、技术文章Markdown文档等海量资料,你想要的我都有!
📣关于我
我是bug菌,CSDN | 掘金 | InfoQ | 51CTO | 华为云 | 阿里云 | 腾讯云 等社区博客专家,C站博客之星Top30,华为云2023年度十佳博主,掘金多年度人气作者Top40,掘金等各大社区签约作者,51CTO年度博主Top12,掘金/InfoQ/51CTO等社区优质创作者;全网粉丝合计 30w+;硬核微信公众号「猿圈奇妙屋」,欢迎你的加入!免费白嫖最新BAT互联网公司面试真题、4000G PDF电子书籍、简历模板等海量资料,你想要的我都有,关键是你不来拿哇。