求二元函数

z=0.2323*x^2-0.2866^2+2*(-0.5406)*a0^2+1.0203*a0^2*x^2/((x^2+y^2)^0.5*tanh(2*(x^2+y^2)^0.5)-x^2*(0.5733-u0)^2)

的最大值,变量x和y都是在0.2附近。求z的最大值,以及x,y的取值。

先用diff命令求z关于x,y的偏导数得到q(1)和q(2)两个方程,

再用fsolve求解方程组,想得到各驻点的坐标,但一直不成功,只得到x=0,y=0;

程序如下:

function q=myfun(p)

x=p(1);

y=p(2);

q(1)=(1046382468744793*x)/2251799813685248 + (2940675037869895*x)/(144115188075855872*(tanh(2*(x^2 + y^2)^(1/2))*(x^2 + y^2)^(1/2) - (2960289404442881*x^2)/9007199254740992)) + (2940675037869895*x^2*((2960289404442881*x)/4503599627370496 + 2*x*(tanh(2*(x^2 + y^2)^(1/2))^2 - 1) - (x*tanh(2*(x^2 + y^2)^(1/2)))/(x^2 + y^2)^(1/2)))/(288230376151711744*(tanh(2*(x^2 + y^2)^(1/2))*(x^2 + y^2)^(1/2) - (2960289404442881*x^2)/9007199254740992)^2);

q(2)=(2940675037869895*x^2*(2*y*(tanh(2*(x^2 + y^2)^(1/2))^2 - 1) - (y*tanh(2*(x^2 + y^2)^(1/2)))/(x^2 + y^2)^(1/2)))/(288230376151711744*(tanh(2*(x^2 + y^2)^(1/2))*(x^2 + y^2)^(1/2) - (2960289404442881*x^2)/9007199254740992)^2) - (1290927876507717*y)/2251799813685248;

main    % 主程序

x0 = [0.2;0.2];   % 初值

options = optimoptions('fsolve','Display','iter');

[x,fval] = fsolve(@myfun,x0,options),

谢谢。

求二元函数最大值matlab,利用matlab, 二元函数求最大值相关推荐

  1. 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序

    环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. YID:856064261479453 ...

  2. 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强

    环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...

  3. 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数

    环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...

  4. matlab求多元函数的极小值,[转载]利用MATLAB求多元函数的极值(2)

    利用MATLAB求多元函数的极值分两种情况,(1)无约束条件:(2)有约束条件. (2)有约束条件下求极小值的方法: 假设多变量非线性函数的数学模型为 min f(x) c(x)<=0 ceq( ...

  5. 利用matlab求零输入响应波形,实验3 利用matlab求LTI连续系统的响应

    实验3 利用matlab求LTI连续系统的响应 一. 实验目的: 1. 了解LTI系统的冲激响应h(t)及matlab实现: 2. 了解LTI系统的阶跃响应g(t)及matlab实现: 3. 了解LT ...

  6. 多元函数的极值matlab,利用MATLAB求多元函数的极值(2)

    利用MATLAB求多元函数的极值分两种情况,(1)无约束条件:(2)有约束条件. (2)有约束条件下求极小值的方法: 假设多变量非线性函数的数学模型为 min f(x) c(x)<=0 ceq( ...

  7. matlab std函数_如何利用Matlab进行小波分析

    小波分析(wavelet analysis), 或小波转换(wavelet transform)是指用有限长或快速衰减的.称为母小波(mother wavelet)的振荡波形来表示信号.该波形被缩放和 ...

  8. j散度matlab,利用Matlab绘制梯度图、散度图、旋度图。.doc

    利用Matlab绘制梯度图.散度图.旋度图..doc 题 目电磁场理论实验姓 名学 号班 级任课老师实验日期2013年 10月 19日 一.实验目的: 1.利用Matlab绘制梯度图: 2.利用Mat ...

  9. 二维静电场有限元matlab,利用MATLAB计算电磁场有关分布

    <利用MATLAB计算电磁场有关分布>由会员分享,可在线阅读,更多相关<利用MATLAB计算电磁场有关分布(17页珍藏版)>请在人人文库网上搜索. 1.电磁场实验报告实验一 模 ...

  10. 物联网控制原理与技术--基于Matlab/利用MATLAB进行频域分析(伯德图)的应用(超详细/设计/实验/作业/练习)

    目录 课程名:物联网控制原理与技术 内容/作用:设计/实验/作业/练习 学习:利用MATLAB进行频域分析(伯德图) 一.前言 二.环境与设备 三.原理 四.内容 五.结果与分析 课程名:物联网控制原 ...

最新文章

  1. WCF系列(四) -- 也谈序列化(上)
  2. 美国空军开发新型机载网络技术
  3. 蛤玮打扫教室(区间覆盖)
  4. struts2和hibernate的简单新闻发布系统_企业做不好新闻软文营销,原因都在这里!...
  5. 项目管理思考——我适合做项目经理吗
  6. React Fiber 了解一下
  7. mysql是否需要设置外键_数据库到底需不需要设置外键?
  8. tomcat8打开闪退_直接双击启动tomcat中的startup.bat闪退原因及解决方法
  9. eclipse 语言包在线更新地址
  10. 2017年多校联合训练 第五场(吉如一)
  11. 2018西门子逻辑控制赛项6部10层电梯PLC
  12. 怎样去除抖音短视频的水印,快速去除视频水印的方法
  13. android自动适应横屏,Android屏幕适配(一)--自定义View屏幕适配
  14. 磁盘被写保护怎么解除
  15. ffmpeg MP4 提取 MP3文件
  16. 基于yolov3的目标检测与LPRnet字符识别的车牌识别(CCPD2020新能源车牌数据集)
  17. springboot项目Banner配置
  18. 新增商品(商品维护模块)
  19. On Premises+Cloud:平滑实现IT应用云转型
  20. MySQL中LAG()函数和LEAD()函数的使用

热门文章

  1. python request 报错 #No JSON object could be decoded
  2. SKUA-GOCAD建模常见问题汇总
  3. 飞猪大数据:85后带来春节旅行新气氛
  4. 深度 | 巨头都在追逐的眼球追踪技术,究竟能带来什么?
  5. 国内技术管理人员批阅google的“春运交通图”项目
  6. 人工智能实验---MNIST图像识别
  7. HTML figure 标签 figcaption 标签
  8. MySQL详细安装步骤
  9. linux脚本出错仍执行后续脚本,crontab执行脚本出错
  10. 浅谈ipad阅读类应用设计