公告：本站正式转型为非交互式静态网站！
转型：本站将通过笔记和博客的形式继续为大家服务，关于 Mathematica 问答服务请移步至QQ群：365716997。
联系：如有问题请联系QQ群管理员，或发送邮件至：lixuan.xyz@qq.com。
感谢：最后非常感谢大家多年来的支持与帮助！
参考：《互联网跟帖评论服务管理规定》《中华人民共和国网络安全法》《网络信息内容生态治理规定》《互联网用户账号信息管理规定》

—— 2022-11-27

欢迎来到 Mathematica 问答社区

● 提问时请贴上文本代码

● 语法高亮：在编辑器中点击

● 被禁止的话题：广告，破解

● 请阅读：《提问的智慧》

● 备用域名：mma.ooo

● 支持LaTex数学公式：
行内公式标识符：\$ 或“$\backslash ($”+“$\backslash )$”，
行间公式标识符：\$\$ 或 “$\backslash [$”+“$\backslash ]$”

● 社区建议QQ群：365716997

538 问题

531 回答

1.5k 留言

229 用户

标签热度

画图 函数 方程 列表操作 微分方程 绘图 算法 化简 wolfram语言入门 编译 方程；函数；符号运算； 符号运算 数值算法 bug 插值函数 表达式 数据筛选 符号，匹配 定义函数 文本 运行 积分 并行 矩阵 随机 动态 拟合 精度 定积分 数据的导入与导出 迭代 替换 数值积分 颜色 图像处理 为什么这个程序的运行结果一直不对呀？ 模式 线性规划 爬虫 极限 读入文件 数组操作 程序包 规则 显示调用 计算 图论-矩阵-文件导出 并行运算 偏微分方程 大数据 excel 输出格式 字符串操作 安装 性能优化 纯函数 colorfunction 选项 manipulate 求最值 基础数学 空间曲面 解方程 变换 字符画 计算控制 g6-图论-矩阵导出 导出 坐标轴 软件打开 隐函数 化简；符号运算 ；列表 子群 内部数据获得 分段函数 多项式 列表运算 网络 统计 字体大小调节 逻辑运算化简 数组赋值 数值 相似矩阵 柱状图 netlink 数值算法，超越方程 分布式计算 复变函数 图像 鼠标 动画 聚类分析 椭圆 未知错误 ndsolve ， 曲线图 傅里叶变换

分类

求目标函数在非线性下的最小值

0 投票
2.4k 浏览
f[x1_, x2_, x3_, x4_] := (((x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) - 8)/
242)^2 + ((6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
7 x4*E^(x4/4) - 34.6)/
465)^2 + ((10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4)) - 114)/
286)^2 + (((1 - (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^
x2) (1 - (1 - 1/E)^x3) (1 - (1 - 1/E)^x4)) - 0.04)/0.21)^2
g[x1_, x2_, x3_, x4_] :=
1 <= x1 <= 10 && 1 <= x2 <= 10 && 1 <= x3 <= 10 && 1 <= x4 <= 10 &&
x1 \[Element] Integers && x2 \[Element] Integers &&
x3 \[Element] Integers && x4 \[Element] Integers

h[x1_, x2_, x3_, x4_] :=
Minimize[f[x1, x2, x3, x4], g[x1, x2, x3, x4], {x1, x2, x3, x4}]
t[{x1_, x2_, x3_, x4_}] := (x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) <=
   250 && (6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
     7 x4*E^(x4/4)) <=
   500 && (10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
     2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
     0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
     2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4))) <=
   400 && (1 - (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^x2) (1 - (1 - 1/E)^
         x3) (1 - (1 - 1/E)^x4)) <= 0.25

 

 

 

我的目的是求f[x1,x2,x3,x4]在满足g[x1,x2,x3,x4]和t[x1,x2,x3,x4]下的最小值
我打算的做法是先求出f的一个最小值判断是否满足t，如果满足就输出，不满足就把这个点去掉再求最小值，求指教这个怎么实现？（我不知道能不能不用循环直接算的，但是Minimize的限制语句好像只能是多项式函数）

For[h[x1, x2, x3,
  x4], (x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) >
   250 && (6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
     7 x4*E^(x4/4)) >
   500 && (10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
     2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
     0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
     2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4))) >
   400 && (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^x2) (1 - (1 - 1/E)^
       x3) (1 - (1 - 1/E)^x4) > 0.25
 , ,
 If[t[x1, x2, x3, x4],
  Print[x1, x2, x3, x4], {Print[x1, x2, x3, x4],
   g[y1, y2, y3, y4] =
    Complement[
     g[y1, y2, y3, y4], {y1 = x1 && y2 = x2 && y3 = x3 && y4 = x4}]}]]

这个不对

 

 

抱歉刚才发布的打错了

t函数是t[{x1_, x2_, x3_, x4_}] := (x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) <=
   250 && (6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
     7 x4*E^(x4/4)) <=
   500 && (10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
     2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
     0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
     2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4))) <=
   400 && (1 - (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^x2) (1 - (1 - 1/E)^
         x3) (1 - (1 - 1/E)^x4)) <= 0.25
用户: zzz7857638 (11 分)
修改于 用户:zzz7857638

1个回答

+1 投票

你这问题完全没必要使用Minimize。

首先修改f和t的形式，仅仅是把输入变量的形式变了，加个括号。

Clear[f, t]
(*修改t和f的输入形式*)

t[{x1_, x2_, x3_, x4_}] := (x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) <= 
   250 && (6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) + 
     7 x4*E^(x4/4)) <= 
   500 && (10^(-5)*1000^1.5*(x1 + E^(x1/4)) + 
     2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) + 
     0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) + 
     2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4))) <= 
   400 && (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^x2) (1 - (1 - 1/E)^
       x3) (1 - (1 - 1/E)^x4) <= 0.25
f[{x1_, x2_, x3_, 
   x4_}] := (((x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) - 8)/
     242)^2 + ((6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) + 
       7 x4*E^(x4/4) - 34.6)/
     465)^2 + ((10^(-5)*1000^1.5*(x1 + E^(x1/4)) + 
       2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) + 
       0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) + 
       2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4)) - 114)/
     286)^2 + (((1 - (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^
             x2) (1 - (1 - 1/E)^x3) (1 - (1 - 1/E)^x4)) - 0.04)/
     0.21)^2

其次，根据g的限定，也就10000种可能，直接暴力枚举所有可能。

all = Tuples[Range[10], 4];
possible = Select[all, t];
data = SortBy[{f@#, #} & /@ possible, First];
First@data

 

用户: 苹果 (2.2k 分)
显示 1 上一条评论
f[{x1_, x2_, x3_, x4_}] 里面的几个变量为什么要加括号呢？
用户: wonderlands0 (451 分)
修改于 用户:wonderlands0
方便后续操作。例如a={1,2}; f[x_,y_]:=x+y; f@@a可以，但是不方便。
用户: 苹果 (2.2k 分)
t[{x1_, x2_, x3_, x4_}] := (x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) <=
   250 && (6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
     7 x4*E^(x4/4)) <=
   500 && (10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
     2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
     0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
     2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4))) <=
   400 && (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^x2) (1 - (1 - 1/E)^
       x3) (1 - (1 - 1/E)^x4) > 0.75
f[{x1_, x2_, x3_,
   x4_}] := (((x1^2 + 2 x2^2 + 3 x3^2 + 2 x4^2) - 8)/
     242)^2 + ((6 x1*E^(x1/4) + 6 x2*E^(x2/4) + 8 x3*E^(x3/4) +
       7 x4*E^(x4/4) - 34.6)/
     465)^2 + ((10^(-5)*1000^1.5*(x1 + E^(x1/4)) +
       2.3*10^(-5)*1000^1.5*(x2 + E^(x2/4)) +
       0.3*10^(-5)*1000^1.5*(x3 + E^(x3/4)) +
       2.3*10^(-5)*1000^1.5*(x4 + E^(x4/4)) - 114)/
     286)^2 + (((1 - (1 - (1 - 1/E)^x1) (1 - (1 - 1/E)^
             x2) (1 - (1 - 1/E)^x3) (1 - (1 - 1/E)^x4)) - 0.04)/0.21)^2
all = Tuples[Range[10], 4];
possible = Select[all, t];
data = SortBy[{f@#, #} & /@ possible, First];
First@data


First::first: "{} has a length of zero and no first element"  就算不出来了
用户: zzz7857638 (11 分)
就那个t 的最后改动了一小下
用户: zzz7857638 (11 分)
@苹果: 明白了，谢谢。
用户: wonderlands0 (451 分)
 | Snow Theme by Q2A Market
Powered by Question2Answer
...