紧急!求助!我正在利用matlab建立三维copula函数模型,正利用极大似然估计方法估计copula函数中的参数,但是卡在这里了。这个似然函数非常的长

下面是我的程序:function siranhanshu=myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)h=(theta^2*(1/theta + 1)*(1/theta + 2))/(U(:).^(theta + 1)*V(:).^(theta +1)*W(:).^(theta + 1)*(1/U(:).^theta + 1/V(:).^theta + 1/W(:).^theta - 2)^(1/theta +3))*Xpdf*Ypdf*Zpdf;

k=prod(h);

siranhanshu=0-log(k);

然后我再命令窗口键入

%在命令窗口输入

[theta,fval]=fminsearch(@(theta) myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf),5.000000001)

其中U,V,W,Xpdf,Ypdf,Zpdf是程序中也求出来的一列多行的矩阵

运行以后报错

Error using  ^

Inputs must be a scalar and a square matrix.

To compute elementwise POWER, use POWER (.^) instead.

Error in myfun (line 3)

h=(theta^2*(1/theta + 1)*(1/theta + 2))/(U^(theta + 1)*V^(theta

+1)*W^(theta + 1)*(1/U^theta + 1/V^theta + 1/W^theta - 2)^(1/theta

+3))*Xpdf*Ypdf*Zpdf;

Error in @(theta)myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)

Error in fminsearch (line 191)

fv(:,1) = funfcn(x,varargin{:});

Error in P_ET0_IR (line 131)

[theta,fval]=fminsearch(@(theta)

myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf),5.000000001)

而且警告Warning: The value of local variables may have been changed to match the

globals.  Future versions of MATLAB will require that you declare

a variable to be global before you use that variable.

> In myfun at 2

In @(theta)myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)

In fminsearch at 191

In P_ET0_IR at 131

新手上路,能否给予指点!谢谢!

matlab怎么根据公式建立三维模型,紧急!求助!用matlab建立三维copula模型,正用极大似然估计方法求copula函数参数...相关推荐

  1. matlab gui输入函数,紧急求助!关于matlab中GUI用户图形界面通过edit输入调用函数问题!...

    是可以的,可以见张志涌<精通matlab6.5版>的一个例子:(图片中的绘图指令输入框就是你说的命令输入框,可以参考这个例子,实现你的功能,主要是eval的使用,用法你可以查看eval的帮 ...

  2. matlab用regress方法求ln函数_高考数学48条秒杀型公式与方法,想要120分一定要掌握...

    来源:高考直通车综合自网络 除了课本上的常规公式之外,掌握一些必备的秒杀型公式能够帮你在考试的时候节省大量的时间,小车这次的分享就是48条秒杀公式,直接往下看! 高考数学48条秒杀公式 1.适用条件: ...

  3. matlab用regress方法求ln函数_数学篇|高中数学48条秒杀型公式与方法,一定要掌握!...

    「 致于学教育 」 高中数学48条秒杀型公式 1.适用条件:[直线过焦点],必有ecosA=(x-1)/(x+1),其中A为直线与焦点所在轴夹角,是锐角.x为分离比,必须大于1. 注上述公式适合一切圆 ...

  4. matlab每两个数插0,【求助】matlab里两个相等的数相减为什么不等于0呢

    程序比较长,数据也很多,截其中一部分吧,具体数据就不写出来了 就是很简单的矩阵运算,求解一个四元方程,得到的解回带回去就不等于0. b(1,1)=-U1(1,5); b(2,1)=-U12(1,5); ...

  5. MATLAB对数图拟合,求助,MATLAB中对数分布拟合问题

    在用蒙特卡洛方法算不确定性,卡壳了.. 思路是这样: 1.利用数据构建概率密度函数,选择拟合优度最大的 2.输入参数进行蒙特卡洛模拟 3.得到数据和概率分布图 之前用水晶球软件算第一步很简单,但是我的 ...

  6. 数理统计仿真实验:大数定律、中心极限定理、矩估计与极大似然估计(含MATLAB代码)

    目录 数理统计仿真实验(Computational Practice) 大数定律(the Law of Large Numbers) 二项分布(Binomial Distribution) 泊松分布( ...

  7. Matlab绘制普朗克公式曲线1/2

    Matlab绘制普朗克公式曲线 1.确定波长和温度区间. 2.运行程序,结果如下. 3.部分源代码. 1900年德国物理学家M.普朗克在量子论基础上建立了关于黑体辐射的正确公式,即普朗克公式. 它给出 ...

  8. 服务器上如何连接显示器和键盘,紧急求助,一个电脑主机能接两个显示器和键盘么?...

    紧急求助,一个电脑主机能接两个显示器和键盘么? 常见问题 可以,不过需要你买一个连接设备.如果你的主机配置的不是很低,就可以这样做.这叫一拖二.最多可以一拖五.具体参数看下面1. 什么是PC插卡套件? ...

  9. matlab传函零极点形式,2013实验一 MATLAB 中控制系统模型的建立与仿真

    实验一 MATLAB 中控制系统模型的建立与仿真 一. 实验目的 (1)熟悉MATLAB 控制系统工具箱中线性控制系统传递函数模型的相关函数. (2)熟悉SIMULINK 模块库,能够使用SIMULI ...

  10. 数据库安全小白紧急求助论坛大神,困扰许久

    数据库安全小白紧急求助论坛大神,困扰许久 mircosoft SQL server 2008 R2怎样才能逃过甚至杜绝外来IP不断尝试SA密码? 求助各位论坛大神,封闭1433端口,设置数据库访问策略 ...

最新文章

  1. 业界丨AI公司魅力榜:谷歌最慷慨排第二,最抠门的你应该能想到
  2. 大学计算机html,编程基础(C+VB+HTML)(辅)19级计算机
  3. 前端JS通过Ajax下载后端返回的Excel文档
  4. Window下UDP(socket)接和收数据案例
  5. write up: web 前女友(SKCTF)
  6. django异常日志_【python小随笔】Django+错误日志(配置Django报错文件指定位置)...
  7. 车林通购车之家--购车计算器模块--保险
  8. django-普通的cookie操作
  9. 首次!阿里达摩院将Pure Transformer 应用于目标重识别ReID!
  10. 一个JS多个数组取交集算法
  11. 隐马尔可夫模型(HMM)简介
  12. 刷屏的北京雾霾,2018 年北上广深空气质量分析
  13. 基于支持向量机的新闻分类
  14. Flutter 闪屏页实现
  15. 《安全评估报告》7条回答范例
  16. XMind 2022 for Mac(思维导图软件)中文免费版
  17. 干货全拿走-用Excel制作小市值轮动价值投资选股器
  18. Kotlin 之 lateinit关键字 与lazy
  19. 软件工程师的职业建议
  20. 计算物理学(数值分析)上机实验答案5、常微分方程初值问题的数值解法

热门文章

  1. Windows系统清理随记
  2. 嵌入式linux全套教学视频(国嵌版)
  3. activiti工作流在线表单设计功能(activiti + ueditor + Ueditor Web Form Design扩展 )
  4. 手把手带你免费申请《软件著作权》 超详细计算机软件著作权申请教程 文末送模板
  5. 89C51流水灯代码示例,带proteus项目
  6. 数字信号处理-04- FPGA常用运算模块-除法器
  7. ASP.NET Web Pages - WebGrid 帮助器
  8. java虚拟机执行过程
  9. 计算机真有趣作文,真有趣作文10篇
  10. 自然语言处理 4.语义分析