🏆本文收录于「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>0fx=xx-b;elseif xx+b<0fx=xx+b;elsefx=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电子书籍、简历模板等海量资料，你想要的我都有，关键是你不来拿哇。
