MATLAB 错误 "Function definitions are not permitted in this context."-6ren

MATLAB 错误 "Function definitions are not permitted in this context."

转载作者：行者123 更新时间：2023-12-01 14:35:34

26

4

<分区>

Possible Duplicate: How to correct “Function definitions are not permitted at the prompt or in scripts”

不足为奇，如果我尝试在 MATLAB 中运行以下 M 脚本，我会收到错误

??? Error: File: kalmanmle.m Line: 47 Column: 1 Function definitions are not permitted in this context.

我不确定这是否可以像我一样运行。或者，我如何在 MATLAB 命令行上运行它？

clear all;

%  State space reprsentation to be forcasted by kalman filter
%   zhi(t+1) = F*zhi(t) + v(t+1)   --> unbobserved varaibles
%   v~N(0,Q)
%   y(t) = A'*x(t) + H'*zhi(t) + w(t)
%   w~N(0,R)

global y;
global x;
global Hvec;
%%----    Enter  Input parameters
F = 0.9;
Q = 0.1;
A = 2;
n = 100;
Hvec = zeros(n,1); %index returns process
indexshock = normal_rnd(0,0.1,n,1);
Hvec(1) = 0;
for i = 2:n,
    Hvec(i) = 0.95*Hvec(i-1) + indexshock(i);
end

%H = 0.3;
R = 0.05;

x = ones(n,1);
zhi = zeros(n,1);
y = zeros(n,1);
zhi(1) = 0;
v = normal_rnd(0,Q,n,1);
w = normal_rnd(0,R,n,1);

H = Hvec(1);
y(1) = A'*x(1) + H'*zhi(1) + w(1);
for i = 2:n,
    H = Hvec(i);
    zhi(i) = F*zhi(i-1) + v(i);
    y(i) = A'*x(i) + H'*zhi(i) + w(i);
end
%% ------------------
%test = [zhi y]

function ret = MyLikelihoodFn(p)
    global y;
    global x;
    global Hvec;
    F = p(1);
    Q = p(2)^2;
    A = p(3);
    R = p(4)^2;
    n = size(y,1);
    P = Q;
    Ezhi = 0;
    Ezhivec = zeros(n,1);
    Ezhivec(1) = Ezhi;
    tmpsum = 0;
    tmp1 = -(n/2)*log(2*pi);
    for i = 2:n,
        yt = y(i);
        xt = x(i);
        H = Hvec(i);
        Ezhi = F*Ezhi + F*P*H*inv(H'*P*H+R)*(yt-A'*xt-H'*Ezhi);
        P = F*P*F' - F*P*H*inv(H'*P*H+R)*H'*P*F' + Q;
        Ezhivec(i) = Ezhi;
        tmpmat = H'*P*H + R;
        tmp2 = -0.5*log(det(tmpmat));
        tmpmat2 = yt - A'*xt - H'*Ezhi;
        tmp3 = -0.5*tmpmat2'*inv(tmpmat)*tmpmat2;
        tmpsum = tmp1+tmp2+tmp3;
    end
    ret = -tmpsum;
endfunction

param = zeros(4,1);
param(1) = 0.2;
param(2) = 0.2;
param(3) = 1;
param(4) = 0.2;

resultparam = fmins('MyLikelihoodFn',param)

actualF = F
F = resultparam(1)
actualQ = Q
Q = resultparam(2)^2
actualA = A
A = resultparam(3)
actualR = R
R = resultparam(4)^2

n = size(y,1);
P = Q;
Ezhi = 0;
Ezhivec = zeros(n,1);
Ezhivec(1) = Ezhi;

for i = 2:n,
    yt = y(i);
    xt = x(i);
    H = Hvec(i);
    Ezhi = F*Ezhi + F*P*H*inv(H'*P*H+R)*(yt-A'*xt-H'*Ezhi);
    P = F*P*F' - F*P*H*inv(H'*P*H+R)*H'*P*F' + Q;
    Ezhivec(i) = Ezhi;
end
test = [zhi Ezhivec Hvec y];
tmp = 1:n;
%plot(tmp,zhi,'-',tmp,Ezhivec,'-',tmp,Hvec,'-',tmp,y,'-');
plot(tmp,zhi,'-',tmp,Ezhivec,'-');

26

4

0

文章推荐： Sesame 服务器对 SPIN 的支持

文章推荐： wpf - 如何在设计时解析 ContentControl ContentTemplateSelector？

文章推荐： json - PagedList 的 Newtonsoft.Json 序列化不包括某些属性

文章推荐： oracle - 使用 Oracle 11g 比较示例表中具有不同行和列的列

python - 如何获取 LaTeX 文件中的所有 `\begin{definition}...\end{definition}` block ？
我刚用 Latex 写完微积分的总结。现在的主要问题是文件中包含许多我现在并不真正需要的东西。 .tex 文件包含许多我需要用心学习的定义和定理。定义在 tex 文件中有自己的定义，因此文件中的任
css - 可以隐藏吗(显示: none;) the definition term (
) in a definition list (
)?
所以我有一个菜单项列表，我想弄清楚我是否应该使用具有类属性的跨度或每个元素的特征的定义列表。以下是我正在考虑的两个选项: 选项 1) // HAML Markup %article.menu-item
C++ 部分模板特化，需要语句 : error: out-of-line definition of 'foo' from class Bar without definition
考虑下面的代码，它试图实现 class Bar 的部分特化。 .在第一种情况下，foo成员函数是内联定义的，在第二种情况下是外联的。行外定义产生了一个我无法弄清楚的编译错误:error: out-of
javascript - 转到 Typescript Source Definition 而不是 Visual Studio Code 中的 Compiled Definition
我正在使用 Visual Studio Code 开发一个 typescript 项目，包括多个结构如下的 npm 包: 源代码:/src/index.ts 编译后的代码:/dist/... 当我右键
definition - 神经网络定义
我正在编写一个神经网络类，遇到了两个我不理解的概念。谁能告诉我bias 和momentum 是什么以及做什么最佳答案偏差是给予神经元的恒定输入。例如在普通的前馈网络中，您可能有 2 个输入单元、2
definition - 重命名特征文件中的步骤时如何让SpecFlow更改步骤定义？
假设我在功能文件中有一个场景，如下所示 Given I log in as "super" user When I click on login Then Home page is displayed
definition - 遗留数据库一词是什么意思？
关闭。这个问题是opinion-based 。目前不接受答案。已关闭去年。已锁定。这个问题及其答案是locked因为这个问题是题外话，但却具有历史意义。目前不接受新的答案或互动。我读了很多这个词
definition - 生产的具体定义是什么？
就像几乎任何一个已经编程了一段时间的人一样，我熟悉“生产代码”这个术语，并且对其含义有一个模糊的认识。然而，有人可以提供一个半严格的定义吗，因为维基百科和谷歌似乎不能？在生产中似乎存在很多灰色地带，例
definition - 伪代码:一个明确的定义？
以下代码是我认为符合伪代码条件的示例，因为它不以任何语言执行，但逻辑是正确的。 string checkRubric(gpa, major) bool brake = false nu
json - 如何从示例JSON生成Swagger #definition
从宠物商店示例中获取以下#definition。给定#definition部分，可以生成JSON结构例如给定一个较大的复杂JSON文件，是否可以做一些相反的事情？给定下面的JSON结构，我可以获
definition - 什么是多态，简单解释一下？
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。想改进这个问题？将问题更新为 on-topic对于堆栈溢出。 7年前关闭。 Improve this qu
definition - 什么是托管可扩展性框架 (MEF)？
很难说出这里问的是什么。这个问题是模棱两可的、模糊的、不完整的、过于宽泛的或修辞的，无法以目前的形式得到合理的回答。如需帮助澄清这个问题以便重新打开它，visit the help center .
definition - "surfacing"用作技术术语时是什么意思？
Joel Spolsky 在 SO 播客中提到了“浮出水面”。这是什么意思？它是类似于“暴露”的东西吗，比如“暴露接口(interface)”？最佳答案这里引用自 podcast 51 的文字记
definition - "semantic"这个词在计算机科学背景下意味着什么？
我不断遇到这个词的用法，但我从来不理解它的用法或所传达的含义。像...这样的短语 "add semantics for those who read" "HTML5 semantics" "sema
definition - "routine"的技术定义是什么？
我正在学习 lisp 语言(做 lisp 例程)，在一般情况下我知道什么是例程，但在技术上下文中我可以谈论它，因为我现在开始学习例程。那么，例行公事的真正定义是什么？(我已经用谷歌搜索过这个，但没有找
Haskell前奏: definition of seq
在 definition of the haskell prelude我们看到 ... 是为无法在 Haskell 中实现的表达式保留的。例如，现在 IO monad 无法在 haskell 中实现。
文件编译时出现multiple definition of 'xxxxxx'的具体解决方法
问题：原来有单个文件tcpclient.c，运行gcc -o tcpclient tcpclient.c可以顺利完成编译，并能与下载到目标板中的tcpserver成功通讯；现在把tcpclien
Emacs clojure : go to definition
我使用 ergoemacs-mode、clojuremode 和自动完成自动完成效果很好。我应该按什么才能进入函数的定义？我习惯了 Cursive，但我想使用免费的 Emacs，我需要去定义功能。
definition - Nearley Tokenizers 与规则
我对 nearly.js 很陌生，我想知道与规则相比，分词器/词法分析器做了什么，根据网站: By default, nearley splits the input into a stream of
definition - 如何在 Forth 定义中创建单词
我正在使用 Gforth ，我想在定义中创建一个词。在 Gforth 的 cmd 行中，我可以输入: create foo ok 或者更具体地说，我定义了一个数组函数，它期望堆栈上的大小并创建一个带有

首页

博学

6Ren·AI

商城

MATLAB 错误 "Function definitions are not permitted in this context."