作业帮有个方程想用MATLAB数值求解 具体如图
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/21 22:58:08
有个方程想用MATLAB数值求解 具体如图
第一个积分式很奇怪,被积变量到底是I还是y0?
而按照第二个式子显示,y0应该是一个常数,只不过取值暂时未知.难道第一个式子在对一个常数做积分?
再问： y0算中间变量
再答： y0到底是中间“变量”还是中间函数？ 我只能理解为y0 = f(P, I) 这样才有可能代入1式中对I积分。 BTW，这又引入另外一个问题，既然对I积分，I怎么能同时做积分限呢？
再问： 这是变限积分的表达呀
再答： 嚯，连积分变量名都懒得换了。。。得，我写个试试。
再问： 不是懒得换 推导出来就是这样 在说明白点 如果能积出来 先带入积分上限的I值算积分 得I的多项式 再代入同一个I的值 当然 这个式子积不出 只能数值求解
再答： “不是懒得换 推导出来就是这样 ” 这句话说的真好。 变上限积分，最基本的一条是，积分界上的变量，在求积分的时候，视作常数。现在你的被积变量和积分界都是I，至少这是表达上的不严谨。 你别用追问了，光浪费钱，我知道你的意思。出结果再说。 程序在附件里。求解的是当P=0.1e6时关于I的变上限积分的的结果，可见在I=1:11的范围内，达不到要求的20e6，无解。
再问： 非常感谢 但能否得出 I与P的图像关系呢
再答： 你最好给我一个靠谱点的解区间，要不然搜索范围太大。 比如上面，P=0.1e6时，在[1, 11]范围内就找不到对应的I。
再问： 我的意思是取一组P或者一组I直接求另一个参数值 这两组区间是无关的 用一组就行
再答： 在P=0.01e6附近，I是有解的。也就是说，现在我能给出一个I = f(P)的函数的数值解，但仅限于P在1e4附近。结果如下图所示。由于使用了穷举法求解，所以速度比较慢还用了并行加速。
横坐标是P，纵坐标是I。
而按照第二个式子显示,y0应该是一个常数,只不过取值暂时未知.难道第一个式子在对一个常数做积分?
再问： y0算中间变量
再答： y0到底是中间“变量”还是中间函数？ 我只能理解为y0 = f(P, I) 这样才有可能代入1式中对I积分。 BTW，这又引入另外一个问题，既然对I积分，I怎么能同时做积分限呢？
再问： 这是变限积分的表达呀
再答： 嚯，连积分变量名都懒得换了。。。得，我写个试试。
再问： 不是懒得换 推导出来就是这样 在说明白点 如果能积出来 先带入积分上限的I值算积分 得I的多项式 再代入同一个I的值 当然 这个式子积不出 只能数值求解
再答： “不是懒得换 推导出来就是这样 ” 这句话说的真好。 变上限积分，最基本的一条是，积分界上的变量，在求积分的时候，视作常数。现在你的被积变量和积分界都是I，至少这是表达上的不严谨。 你别用追问了，光浪费钱，我知道你的意思。出结果再说。 程序在附件里。求解的是当P=0.1e6时关于I的变上限积分的的结果，可见在I=1:11的范围内，达不到要求的20e6，无解。
再问： 非常感谢 但能否得出 I与P的图像关系呢
再答： 你最好给我一个靠谱点的解区间，要不然搜索范围太大。 比如上面，P=0.1e6时，在[1, 11]范围内就找不到对应的I。
再问： 我的意思是取一组P或者一组I直接求另一个参数值 这两组区间是无关的 用一组就行
再答： 在P=0.01e6附近，I是有解的。也就是说，现在我能给出一个I = f(P)的函数的数值解，但仅限于P在1e4附近。结果如下图所示。由于使用了穷举法求解，所以速度比较慢还用了并行加速。
横坐标是P，纵坐标是I。