a*******1
发帖数: 1554 | 1 你看这样可以不？
S0是初始价格，r是利率，sigma是波动率，T是到期日，K是执行价，hsigma是算vega用
的，h0是算delta用的，N是200，p是10^8，但一般不用取这么大吧，我读书做project
才取10万，这是当时作业的一小部分
function []=ThreeMethodsMilter(S0, r, sigma, T, K, h0, hsigma,p,N)
%% Finite difference
SqrRootp = sqrt(p);
DeltaT = T/N;
SqrtDeltaT = sqrt(DeltaT);
Z= randn(N-1,p);
P = prod([ones(1,p);1+r*DeltaT+sigma*SqrtDeltaT*Z+0.5*sigma^2*DeltaT*(Z.^2-1
)]);
S = P*S0;
C = exp(-r*T)*max(0,S-K);
SPlus = exp(-r*T)*max(0,(S0+h0)*P-K);
SMinus = exp(-r*T)*max(0,(S0-h0)*P-K);
DeltaFD = (SP... 阅读全帖 |
|
r*t
发帖数: 34 | 2 I just write a program of DOPRI5(4), but have a problem with automatic step-
size selection. after calculation of K1~K7, the error is estmated by:
Error = (71/57600*k1 - 71/16695*k3 + 71/1920*k4 - 17253/339200*k5 + 22/525*
k6 - 1/40*k7)*DeltaT (eq. 1)
Then I adjust DeltaT by:
Theta = 0.9*(Tolerance/Error)^(1/5) (eq.2)
(Here I select Tolerance = 1e-7)
and if Theta<1, then DeltaT = DeltaT*Theta
However, I found that the simulation is exemely slow (much slower than
without step-size control), |
|
J***A
发帖数: 1511 | 3 我知道是求导 问题是求导以后是除以零度时的v0 还是现在温度的体积V, 这个很关键
， 即便是纯线性的。 V = V0 a* deltaT 这个绝对不能用V0 来近似的。 即便可以
， 那说明a 或者deltaT足够小， 如果Delta 足够小， 所有的V都可以用V0 代替 不局
限为线性。
除非wiki是错的 否则这个V和V0 失之毫厘谬以千里。 当然如果定义就是V我觉得在线
性关系下， 膨胀关系竟然反比于温度， 我挺难以接受的。
前面用dV/dT就是求导的意思,否则就该用delta ...而且我前面的"科普"里是必须用求
导(温度变化很小的测量比较的意思),这样才能选出所需的线性材料.delta V = ......
.. |
|
J***A
发帖数: 1511 | 4 我们说的东西总不在一个频道， 不知道回帖的时候都可有认真看帖。
我不认为Wiki中用V0代替V是一种近似， 人家明确说除以摄氏零度的体积。 因为人家
认为normalize by 物质的量就可以
常用的固体 膨胀系数接近于常数， 如果按照V来定义， 并不能得出和理想气体一样的
均匀刻度。
理想气体
delta V/V/ deltaT 是1/ V 也是1/ T而不是我常数
而wiki 定义中固体的 deltaV/ V/delta 是常数， 但由这个常数 要积分 得到 温差是
exp （ deltaT ) 这就得泰勒展开求delta T一次项得近似膨胀系数了吧
更简单的算法是
直接把 除以V变成除以V0 ， 因为固体液体温度计V和V0差不多， 这样根本不用积分那
么麻烦， 直接可得线性关系， 也就是普通温度计刻度均匀的原理。
至于为何体膨胀是线膨胀系数三倍， 按happy life 说的， L加Delta L 取立方 在略
去deltaL的二三次方项即可。
顺便和某月光说下， 现在是信息时代了， 生活工作用Google 用wiki 能学到很多东西
。 看版上讨论， 自己wiki 我非物理系... 阅读全帖 |
|
g*z
发帖数: 124 | 5 If you are fighting ground bounce, you should use some capacitors with ESR
and ESL as small as possible, the capacitance value is not that important.
You can start with some X7R 0603 SMD caps(0.1uF or 0.01uF) or even smaller (
0402), put them as close as possible to the power pins. The keyword here is
ESL and ESR.
The gound bounce voltage V=ESL * deltaI/deltaT
for digital circuits, when switching happens, deltaI/deltaT can be huge, little inductance can cause huge voltage bounce. |
|
A****s
发帖数: 932 | 6 另外空间膨胀时同时产生物质，空间是正能量，物质是负能量，总能量为零，所谓免费
的午餐。
因为总能量为零，时间无穷大，按照海森堡测不准原理，DeltaE x DeltaT > h/2，这
个宇宙完全是可以从量子泡沫中自发产生的。 |
|
b*****g
发帖数: 2727 | 7 4年前买房子时候房检说deltaT是10度，但不知道是怎么测的。但是过去三年我用测温
枪不断在测，一直是9度的样子。房子可以冷下来只是时间长。如果有大漏空调应该完
全不制冷。
另外，我压缩机原厂标称8.5磅制冷剂，如果加上管道我胡乱猜应该总共10磅。他加了3
磅，说明原来有7磅。7磅：10磅 离 50psi:75psi 不远。所以加了之后应该没有比标准
低太多。 |
|
s******2
发帖数: 852 | 8 测量室内外温湿度几个基本参数，再和厂家的performance sheet对照高低压，deltaT
，电流（功耗）等，就知道了。你可以根据你的机器型号搜索一下这个performance
data。
但一般来说，ΔT在20以下很多是不正常的，除非相当潮湿 |
|
l**********7
发帖数: 1401 | 9 妈妈今天回国，DELTA从纽约到底特律，等一个小时，然后到上海。我在纽约送她上飞
机后就回办公室了，然后接到我妈妈的电话，说她坐上飞机了，但因为天气原因飞机不
能起飞，不知道要等多久。后来我自己查了DELTA网站，晚了2个小时才起飞，后一班飞
机肯定赶不上了。
现在从底特律到上海的飞机已经起飞了，再打给DELTA，说这班飞机上旅客名单中没有
我妈妈的名字，根据她的CONFIRMATION NUMBER，她被安排到明天的同一班航班，所以
她要在底特律等24个小时。打我妈妈手机又不接。我问DELTAT这种情况他们负责不负责
住宿，他们说因为天气原因不负责。
我现在真的急死了，都不知道妈妈身上带了多少钱，够不够住酒店。我怎么办呀！！！
*******************************************************************************
刚和妈妈联系上了，她在酒店发EMAIL给我的，说她手机没电了。DELTA已经给所有乘客安排了住宿和两顿饭。这下我放心了，谢谢大家。 |
|
g********d
发帖数: 19244 | 10 那个会中文的柜台小姐 简直就是个傻逼
妈妈今天回国，DELTA从纽约到底特律，等一个小时，然后到上海。我在纽约送她上飞
机后就回办公室了，然后接到我妈妈的电话，说她坐上飞机了，但因为天气原因飞机不
能起飞，不知道要等多久。后来我自己查了DELTA网站，晚了2个小时才起飞，后一班飞
机肯定赶不上了。
现在从底特律到上海的飞机已经起飞了，再打给DELTA，说这班飞机上旅客名单中没有
我妈妈的名字，根据她的CONFIRMATION NUMBER，她被安排到明天的同一班航班，所以
她要在底特律等24个小时。打我妈妈手机又不接。我问DELTAT这种情况他们负责不负责
住宿，他们说因为天气原因不负责。
我现在真的急死了，都不知道妈妈身上带了多少钱，够不够住酒店。我怎么办呀！！！
****************************************************************************
***
11月4日更新
刚和妈妈联系上了，她在酒店发EMAIL给我的，说她手机没电了。DELTA已经给所有乘客
安排了住宿和两顿饭。这下我放心了，谢谢大家。
*****... 阅读全帖 |
|
m********e
发帖数: 585 | 11 It's easy: assume two equivalent batteries in parallel, the voltage is the
same as one battery. L*delta(I/2)/delta T *2 = L *deltaI/deltaT
or
... |
|
h********y
发帖数: 2 | 12 I have a multi input coupled 2nd order non-linear system and want to generate
a
S-funciton generic simulink block. The forward numerical method run well in m-
file formate if given a very small step size: deltaT=.00005, but it will blow
up when implemented by using s-function block.
I don't have much experience using m-file to write s-function. I am wondering
if the m-file s-funciton itself have some bugs when impelmenting this kind of
highly non-linear system.
Help help. and Thanks a lot. |
|
h********y
发帖数: 2 | 13 I have a multi input coupled 2nd order non-linear system and want to generate
a
S-funciton generic simulink block. The forward numerical method run well in m-
file formate if given a very small step size: deltaT=.00005, but it will blow
up when implemented by using s-function block.
I don't have much experience using m-file to write s-function. I am wondering
if the m-file s-funciton itself have some bugs when impelmenting this kind of
highly non-linear system.
Help help. and Thanks a lot. |
|
k*****e
发帖数: 229 | 14 只是估算，实际情况没有办法完全定义好，所以恒温确实是有些歧义，抱歉。
最后的大致估计是利用公式：
Qin=Qout=C*m*deltaT
冷水端：Q=UAdeltaT (U:传热系数，A:传热面积，包含了长度L)
两者相等，最后得到估计值。奇怪的是，长度居然和冷水流速成正比，恐怕计算不太正
确。 |
|
c****c
发帖数: 29 | 15 我现在又晕了，比如CIR model，如果求rate的volatility，那么不需要考虑drift的部
分，也就是说只需要考虑sigma*sqrt(deltat)*sqrt(rt)*et. et是std为1的正态分布，
不用管。那么就是说只要用historical data求得r(t+1)/sqrt(rt)的std就可以了？
vol |
|
