假设忽视了象征性的PDF

2视图(30天)
Lea
Lea 2023年2月26日
评论道: 保罗 2023年2月27日
我要定义 随机变量的概率密度函数z。 累计密度函数。
作为 , 如果 的概率密度函数是z。
然而,Matlab忽略这个假设(见下面的代码)。为什么会这样,我怎么能修复它吗?
感谢任何帮助:)
%设置符号
信谊x yβdb直流d
信谊z (z) f (z)% z, PDF和运作
信谊理查德·道金斯(z) rp (z)%其他函数依赖于z
%假设参数
假设(x > 0)
假设(y > 0)
假设(β* x < y)
假设(β< 1 &β> 0)
假设(F (z) < = 1 & F (z) > = 0)% z的运作
% z的PDF
f (z) = = diff (f (z), z)
ans =
假设(f (z) > = 0)
假设(int (f (z), z, 0, 1) = = 1)
警告:忽略假设常数。
% * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
Vb =日志(x-db) +β* int ((z *日志(rp (z) * db) * f (z)), z, 0, 1)
+β* int (((1 - z) *日志(rd (z) * db) * f (z)), z, 0, 1)
Vb =
% * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
%区分Vb w.r.t. db和设置为0
diff_Vb_db = diff (Vb, db) = = 0
diff_Vb_db =
隔离(diff_Vb_db db)
ans =
3评论
沃尔特·罗伯森
沃尔特·罗伯森 2023年2月26日
假设(int (f (z), z, 0, 1) = = 1)
你有含蓄地告诉MATLAB, f是一个变量的函数。一个变量的函数的定积分是一个常数。常数= =常数symtrue或symfalse,假设()symtrue或symfalse被忽略。
如果你宣布f是多个变量的函数,然后int()不会被认为是一个常数。

登录置评。

接受的答案

保罗
保罗 2023年2月26日
我不知道这样的方式把一个假设在一个积分。
在这个问题上,似乎f (z)的实际形式并不重要,所以我想你可以分配f (z)的任何函数满足的积分约束。
或者,这些操作工作,但我不确定他们会推广。我不得不玩有点表情到潜艇命令的正确形式的工作
信谊f (z)
信谊dbβx真实
assumeAlso (db ~ = 0);
assumeAlso (db ~ = x);
eqn =β* int (f (z) * (z 1) / db, z, 0, 1) +β* int (f (z) * z / db, z, 0, 1) + 1 / (db - x) = = 0
eqn =
eqn =简化(eqn)
eqn =
eqn =扩大(eqn)
eqn =
eqn =潜艇(eqn, int (f (z), z, 0, 1, 1)
eqn =
eqn =β/ db -简化(解决(eqnβ)/ db) = = 0
eqn =
12个评论
保罗
保罗 2023年2月27日
我同意 假设 不允许假设积分是一个特定的常量,这也是为什么上面的原始解决方案发布潜艇用于积分(在一些操作表达式为潜艇会工作)。但是,我想知道这是预期的行为基于文档页面链接,说:
“工具箱不支持假设象征功能。万博1manbetx符号变量和假设 表达式 相反。”(emphasis added)
让我们来测试这用一个简单的情况
信谊x y
symType (x * y)
ans =“表情”
我们看到x * y是一个表达式,所以它应该可假定的
假设(x * y = = 1);
假设
ans =
它是。和假设收益率预期的结果
简化(x * y)
ans =
1
现在尝试int
信谊f (z)
symType (int (f (z), z, 0, 1))
ans =“表情”
还一个表达式。但这个表达式不是可假定的
假设(int (f (z), z, 0, 1) = = 1)
警告:忽略假设常数。
假设
ans =
有什么区别x * y, int(…)表达式作为假设()而言?

登录置评。

更多的答案(0)

s manbetx 845


释放

R2021a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!