L*M
发帖数: 8 | 1 Anyone who has Mathematica, maple, or other tools, lease help me solve
following equations:
(Yc-Y0)**2 + (Xc-X0)**2 = r**2
(Y1-Y0)**2 + (X1-X0)**2 = (r-r1)**2
(Y2-Y0)**2 + (X2-X0)**2 = (r-r2)**2
r, X0, Y0 are variables. r1, r2, X1, Y1, X2, Y2, Xc, Yc are paremeters. |
|
i******r
发帖数: 164 | 2 一个二元且有限的concave function：f(x，y)，在(x0,y0)处的左右偏导不相等：
f对x求偏导，从左逼近x0为a，从右逼近x0为b，a>b(因为f是concave)
f对y求偏导，从左逼近y0为c，从右逼近y0为d，c>d。
a、b、c、d存在且有限
那么在(x0,y0)处的subgradient,定义为(s1,s2），是不是
b<=s1<=a, d<=s2<=c
谢谢。
|
|
R*********r
发帖数: 1855 | 3 上面所说的所谓不是“足够好的”分布只有一类：
f(x)=\sum a_n \delta(x-n x0)
其中a_n<>0所对应的n的最大公约数是1
对于这种分布，E(s)=(x0[s/x0]-x0/2)/E1+E2/2E1^2
其中取整为向上取整，而且等式是严格成立的。
对其它的任何分布都有E(s)=s/E1+E2/2E1^2+r(s)
r(s)是指数衰减的。
比如f(x)=\delta(x-1)/2+\delta(x-pi)/2就已经属于足够好的分布范畴，但是f(x)=\
delta(x-1)/2+\delta(x-3)/2不是。 |
|
G******i
发帖数: 163 | 4 This does not necessarily hold, if we only assume the conditions you gave.
> Add a small term k*a*(1-a) to x0, for k small enough, I claim the
inequality
is still true, as follows,
a*U(x0+k*a*(1-a),y2)+(1-a)*U(x0+k*a*(1-a),y1)>U(x0,a*y2+(1-a)*y1). |
|
H****h
发帖数: 1037 | 5 椭圆上点(x1,y1)的外法向量方向是(x1/a^2,y1/b^2).
(x1,y1)和(x0,y0)的联线方向是(x0-x1,y0-y1).
你需要两者方向一致，于是有
(x1/a^2)/(x0-x1)=(y1/b^2)/(y0-y1)>0
所以有：(x1/a^2)(y0-y1)=(y1/b^2)(x0-x1).
得PQ |
|
r**g
发帖数: 120 | 6 BTW,
if a is also a non-negative vector, then the minimizer x0 is:
x0 = k A^{-1} a
where the constant k is the solution to the equation
x0' A x0 = c.
is |
|
c****s
发帖数: 51 | 7 一个椭圆,假设方程是:
A + B*(X-x0)*(X-x0) + C*(Y-y0)*(Y-y0) + D*(X-x0) + E*(Y-y0) + F*(X-x0)*(Y-y0
) = 0 (由于F不为零,所以X和Y之间有correlation)
请问这个correlation怎么样用A, B, C, D, E, F 这些参数来表示. |
|
r*******s
发帖数: 303 | 8 define f(t) = int_0^t k(s) ds
define X(t) = exp(-f(t))X0(t)
then dX(t) = -exp(-f(t)) df(t) X0(t) + exp(-f(t)) dX0(t)
============-X(t) df(t) + v(t) dw(t)
==> dX0(t) = exp(f(t))v(t) dw(t)
X0(t) = X0(0) + int_0^t exp(f(s))v(s) dw(s) |
|
S*****H
发帖数: 90 | 9 Outer search: binary search for F(x)=1/2. Starting with interval [a,b], try
x0=(a+b)/2.
The inner search is to calculate F(x0). For each array, inner search is
needed to count the frequency of <=x0. The sum of the frequency divided by
nk is F(x0). |
|
h*l
发帖数: 19 | 10 1. 给定一个坐标点集,已知连线是凸封闭的,
求实际连线序列.
答: 求给定点集的凸包的算法是熟知的,如查
"ALGORITHM DESIGN MANUAL”, by … at
SUNY. 所得凸包即连线序列.
2. 给定两个凸封闭连线序列,判定有无重叠.
答: 两个凸多边形不重叠当且仅当存在其中
某一多边形的某一边所在直线使二多边形在其
两侧, i.e, at least one direction can
see through. (^_^) , 证略.
这样, 设该序列点集为D1: (x11,y11),
(x12,y12)….和 D2:(x21,y21),(x22, y22)
… 则循环求其中是否存在”边”(x0,y0),
(x1,y1)
使得任给(x1k,y1k)属于D1, (x2k,y2k)属于D2有
| x1k y1k 1 |
| x0 y0 1 |
| x1 y1 1 | 皆 <= 0 而
| x2k y2k 1 | 皆 >= 0
| x0 y0 1 |
| x1 y1 1 |
或反之. 如存在, 则不重叠. (当仅有
公共边界点也算重叠时,则除(x0,y |
|
b*********h
发帖数: 46 | 11 知道了，假定x，y都是joint gaussian，各自的pdf
～exp(-1/2(x-x0)'*inv(Cx)*(x-x0))
和exp(-1/2e'*inv(Ce)*e), 假定x，y相互独立，joint pdf
～exp(1/2(x-x0)'*inv(Cx)*(x-x0)-1/2e'*inv(Ce)*e),
取log并最大化，就是最小化(Ax－y)'*inv(Ce)*(Ax－y)'
+x'*inv(Cx)*x, 原文错了，不是Cy是Ce，误差的covariance
数学太差了，居然搞了好几个小时。。。 |
|
l**i
发帖数: 8144 | 12 简单地讲 就是谋事在人 成事在天
认为社会发展沿着某个预定方向发展的理论 都是傻逼理论
--------------------------------------
Anthropic principle
From Wikipedia, the free encyclopedia
In astrophysics and cosmology, the anthropic principle is the philosophical
argument that observations of the physical Universe must be compatible with
the conscious life that observes it. Some proponents of the argument reason
that it explains why the Universe has the age and the fundamental physical
constants necessary to accommodate conscious life. As a ... 阅读全帖 |
|
T*******x
发帖数: 8565 | 13 这个其实很简单：假设G和H距离为a，根据x方向上的一致收敛性，存在X0，Y0，当x>X0
，y>Y0时，f(x,y)离G距离不到a/2，根据y方向的一致收敛性，存在X1，Y1，当x>X1，y
>Y1时，f(x,y)离H距离不到a/2，找一个点(x,y)同时满足x>X0，y>Y0，和x>X1，y>Y1，
那么f(x,y)既和G距离不到a/2，又和H距离不到a/2，矛盾。
y |
|
发帖数: 1 | 14 You are given an object with the following state and method:
{
private constant int x0 = ?
private constant int v0 = ?
private int i = 0
public boolean test(int x)
{
++i
return x == x0 + i * v0
}
}
Find out the value of x0 and v0,
by invoking test(x) only, for a finite number of times. |
|
c********k
发帖数: 134 | 15 二手交易风险自负！请自行验证是否合法和一手卡！:
我想卖的物品:
$25x20 GAP (x0.82)
$25x7 Old Navy (x0.81)
$25x12 JC Penney (x0.82) |
|
x*******2
发帖数: 5333 | 16 二手交易风险自负！请自行验证是否合法和一手卡！:
我想卖的物品:
$250 GAP (x0.83)
$25x10 Old Navy (x0.83)
$500 JC Penney (x0.85)
$300 options .85
单张面值:
$25
可接受价格(必须明码标价！):
as above
物品新旧要求:
new in hand
邮寄方式要求:
CODE ONLY - FREE
First Class Mail + Delivery Confirmation
买卖双方谁承担邮寄损失(Required if not code only):
before me, after you
付款方式说明:
non-cc paypal or BOA
其他补充说明:
广告的有效期:
物品来源(Required for All Cards!):
EPG
我的联系方式:
mitbbs
Warranty期限:
能否证明是合法的一手卡?(Required for All Cards!):
yes. |
|
f**********a
发帖数: 2885 | 17 ebay
300块x0.96
safeway
125块x0.92
radio shark
125块x0.86
dell
600x0.92
barnes noble
250x0.86
baby r rus
750x0.9
sally beauty
50x0.80
谢谢 |
|
c********k
发帖数: 134 | 18 二手交易风险自负！请自行验证是否合法和一手卡！:
我想卖的物品:
$250 Barnes & Nobel (x0.82)
$200 JC Penney (x0.82)
$450 Foot Locker (x0.82) |
|
c********k
发帖数: 134 | 19 二手交易风险自负！请自行验证是否合法和一手卡！:
我想卖的物品:
$25x20 GAP (x0.82)
$25x7 Old Navy (x0.81)
$25x12 JC Penney (x0.82) |
|
x*******2
发帖数: 5333 | 20 [出售] GAP, Old Navy, JCPenney
发信站: BBS 未名空间站 (Fri Jun 18 14:16:26 2010, 美东)
二手交易风险自负！请自行验证是否合法和一手卡！:
我想卖的物品:
$150 GAP (x0.83)
$25x10 Old Navy (x0.83)
$500 JC Penney (x0.85)
$300 options .85
单张面值:
$25
可接受价格(必须明码标价！):
as above
物品新旧要求:
new in hand
邮寄方式要求:
CODE ONLY - FREE
First Class Mail + Delivery Confirmation
买卖双方谁承担邮寄损失(Required if not code only):
before me, after you
付款方式说明:
non-cc paypal or BOA
其他补充说明:
广告的有效期:
物品来源(Required for All Cards!):
EPG
我的联系方式:
mitbbs
Warranty期限:
能否证明是合法的一手卡?(Required |
|
n*******w
发帖数: 687 | 21 bless！
1.写一段程序比较两棵树是否一样。
常见题。
2.有一个奇怪的linkedlist,除了next pointer还有一个random pointer指向一个随机
的节点。问如何实现clone函数。
最近版上刚讨论过。先creat big linkedlist然后split。
3.写一段代码，给一个字符串，例如"30*(5+10)",输出计算结果。
经典算法。两个stack，一个操作数，一个操作符。写代码其实不简单，要定义操作符
的优先级。
4.写一段代码，输入一个数组和一个数字，找两个数组元素和为给定数字。
经典题。允许O(n)空间，hashtable。否则先sort，一前一后两个指针往中间找。
5.输入一个linkedlist和一个数字例如：9->7->8->6->1->2 和 3，输出还是一个
linkedlist但是每三个数reverse一下，例如8->7->9->2->1->6。
版上最近刚讨论过。递归或者iterative都有。
6.如何求一个树的mirror（将所有节点的children节点反序排列）.
跟遍历similar。
7.下面这道题目是吃饭的时候问的... 阅读全帖 |
|
O******i
发帖数: 269 | 22 看去比二分法求方程的根要难些，应该怎么做？为什么最坏情况是O(n)呢？
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
现在给定一个函数，f(x),x在某值之前是非递减的，在某值之后是非递增的。设计一个
算法快速查找这个值。我给出的算法最坏情况是O(n)，平均是O(lgn)，但是很可惜没来
得及写完。
跟离散的情况类似。比如一个数组n个元素，先递增然后递减。找到最大元素的index。
类似binary search。版上最近也讨论过。连续的情况，给定x0，要比较f(x0)与f(x0+
delta)的大小然后binary search。delta是答案的精度。 |
|
t****a
发帖数: 1212 | 23 不好意思，之前的方程写错了，少写了开根号
修改以后的偏导数方程组为
-sum((x[i]-x0)/f(x0,y0)) = 0
-sum((y[i]-y0)/f(x0,y0)) = 0
我不会求解这个方程组，无法给出解析解
只好用梯度下降方法来迭代求解数值解。
这一题好像是某一年google code jam某一轮列出的题目之一。我记得当时那道题目也
只要求精度若干的数值解。 |
|
B***i
发帖数: 724 | 24 安慰一下
1. 点 (x0, y0, z0), norm (a, b, c),
(x - x0, y - y0, z - z0) * ( a, b, c ) = a ( x - x0) + b ( y - y0) + c (
z - z0) = 0
=> ax + by + cz + ( -ax0 - by0 - cz0) = 0 |
|
P**********k
发帖数: 1629 | 25 你再好好想一想吧
完全不用计算点到直线的距离。
假设(x0, y0)到直线L1的距离是最短的，那么假如与点(0, y0)距离最短的直线不是L1
，而是L２，那说明L1和L2必定在(0, x0)这个区间有交点。
但是题目给的条件是所有直线在这个区间都不相交，所以到(0, y0)距离最短的直线和
到(x0, y0)距离最短的直线是一条。
距离 |
|
x***1
发帖数: 999 | 26 happy moon festival!
V是t的函数，V0 是温度等于0时的V, V0=V(t=0).
那么V'(t=0)=dV/dt(t=0).也就是求导以后在那点的值.对不对？
你看近似：f(x0+dx)=f'(x0)*dx+f(x0) |
|
x***1
发帖数: 999 | 27 happy moon festival!
V是t的函数，V0 是温度等于0时的V, V0=V(t=0).
那么V'(t=0)=dV/dt(t=0).也就是求导以后在那点的值.对不对？
你看近似：f(x0+dx)=f'(x0)*dx+f(x0) |
|
r*****y
发帖数: 53800 | 28 新浪体育讯 亚洲杯是亚洲足坛最重要的赛事。由于世界杯亚洲区预选赛的分组将
按照国际足联的FIFA积分排名来确定种子队。因此，亚洲杯更是亚洲诸强事关世界杯预
选赛能否成为种子队的“重要抢分期”。
中国队在小组赛3战全胜的表现创造了亚洲杯历史来的最佳战绩。其实际意义也直
接体现在了抢分成果上。按照亚洲杯小组赛结束后的赛果计算，本届赛事前亚洲排名第
9位的国足，排名积分已大幅超过身前的卡塔尔、阿曼以及约旦；在4月1日的2018世界
杯预选赛亚洲区分组抽签(40强赛阶段)时，国足必将以亚洲前8的身份成为“种子队”
。这样，国足在亚洲区40强赛阶段突围进军12强赛的希望将大增。
2018世界杯亚洲区预选赛将分为三阶段，国足将直接参加第二阶段，也就是40强赛
阶段。40支球队分为8个小组，每组头名和4个成绩最好的小组第二名将晋级最终阶段，
也就是12强赛阶段。
FIFA世界男足积分排名的计算规则非常复杂。具体公式是——
某队积分=该队(12个月内比赛平均积分×1)+(12个月前至24个月内比赛平均积分X0
.5)+(24个月前至36个月内比赛平均积分X0.3)+(36个月前至48个月内比赛平均积分... 阅读全帖 |
|
n********r
发帖数: 4558 | 29 这些数据遵循指数衰减
y=Aexp(-x/x0)+y0
拟合结果（保留到小数点后一位）：
A=330.6
x0=2.8
y0=48.4
拟合参数的物理意义：
A代表某排名有可能达到的极限最长时间，330.6周，正常人类不可能突破。
x0代表各档次区间，排名每隔2.8位下降一档。
y0代表排名的一般保持时间，48.4周，接近52周，可见目前实行的52周更新的
排名系统乃是精确计算的结果。
Q.E.D. |
|
b*****e
发帖数: 762 | 30 好像应该假设x=vt+x0.
所以 t=(x-x0)/v
Y(x,t)=Asin(kx)cos[k(x-x0)]
这个还是个sin形状。所以一定要考虑谐波。假设弦是线性的。那么高频的谐波的k1,k2
都是k的整数被。不妨设k1=2k,k2=4k.
所以三项叠加。我用mathemacica画了一下 跟那个还有点像。当然具体实际上cmos扫描
速度很快，可能有好几个周期，并且可能有个死时间，这个应该能解释为什么它很多时
候振幅正负不对称。
看图 |
|
E*******e
发帖数: 1309 | 31 你理解错误了。
图中的时刻水还在留，所以x3小于x0.水流完的一刻x3会大于x0，然后最后等于x0.
自己做个试验，在称上的时候，让自己失重，称的读书马上减小。然后把身体定住，马
上进入超重状态，称的度数大于体重，然后恢复到体重。
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~你忽略了，空气
中有一部分水，它离开了上面的容器，确还没有碰到下面的容器 |
|
a*********a
发帖数: 413 | 32 高中时代全国物理一等奖告诉你们这是x0=x1=x3>x2
受力分析：
x0=x3 参见物理竞赛题里面一个经典的卖大米的问题：麻袋放在秤上，大米以恒定流量
流出闸门,经H高度自由落体，秤示数为W时闸门关闭，问袋子里有多少米？答案是W，因
为空中的米重量刚好跟由于米粒撞击带来的额外示数相等。
x1和x2 首先考虑液体内部的压强朝任何方向都是一样的，所以虽然有能量的流失，但
是提供水射出的力可以使任何方向的力。对漏水孔简单的分析告诉我们：△m*v=F*△t
，而且这个速度和射出的水受到的合力是一样的方向。既然这个速度不能忽略那么这个
力就不能忽略。x1水平射出，这个合力最终来自地面的摩擦力。x2垂直向下射出，那么
，水受到自重和秤的支撑反应出来的合力向下，支撑力小于重力，x1>x2
对于不能理解x0=x1的同学可以对容器底部进行受力分析：容器底部液压处处相等，受
到的来自液体的总的压力为压强乘以面积，两者是一样的。 |
|
p********n
发帖数: 98 | 33 x0=x3>x1>x2
period
一个简单的假设,就是认为喷水已经很长时间了,x1只有一半水,x2也只有一半水,还有反
冲.x0满水
x3满水,所以x0=x3>x1>x2
有必要像这么多人讨论的那么复杂嘛?奇怪 |
|
l*3
发帖数: 2279 | 34 呵呵, 我觉得你可能是不会证明这个显而易见的事实:
"刚体动能可以分解为随质心平动的动能+相对质心转动的动能" 这两部分.
我还是给你证明一下吧:
设初始时刻, 质心位置是r0, 刚体上任意一点位置r, 对应的密度是rho(r)
(注意, 这里r是指矢量, 是一个坐标组, 其实就相当于n维向量, 用r 是为了表示其状
态是 "初始状态", 和下面要提到的x做区分)
x0=x0(r0,t) 表示质心的运动轨迹 (依赖于初始值r0, 时间t)
x=x(r,t) 表示刚体上初始时刻位置为r的点的运动轨迹 (依赖于r, t, t=0时x=r)
则质心的定义意味着下式成立:
int (x-x0) rho(r) dr = 0
两端对 t求微分得到:
int (v-v0) rho(r) dr = 0
其中v=v(r,t) 表示陀螺上初始位置在r的点在t时刻的速度. v0=(r0,t) 表示质心在t时
刻的速度.
而陀螺总动能的表达式是什么呢?
是1/2( int v^2 rho(r) dr )
上式 =1/2( int 1/2 (v-v0 + v0)^2 rho(r) dr)
= 1/2( int ... 阅读全帖 |
|
t******n
发帖数: 2939 | 35 ☆─────────────────────────────────────☆
l63 (l63) 于 (Sun Jun 30 11:44:29 2013, 美东) 提到:
具体是这样: 如图, 陀螺理想化模型为一个零质量的杆顶着一个均匀质量的圆盘, 圆盘
重心为C, 陀螺支撑点为O, 固定O点 (就是只是固定O点的位置, 陀螺杆本身是可以往任
意方向摆动的), 初始状态时, 陀螺绕杆的角速度为w0, 杆与竖直方向有一个小的夹角
a0, 初始状态时陀螺重心速度为0. 杆长为D.
计算方法: 建立坐标系, 以柱坐标描述重心的位置C=C(r,theta,z), 由于杆长固定, 故
z=sqrt(D^2-r^2), 相当于重心是两个自由度: r, theta (也就是可以用极坐标描述),
t=0时, theta=0, theta'=0, r=D*cos a0, r'=0; 陀螺还有一个自由度, 是绕杆的角速
度w, 这个w在t=0时的值为w0.
考虑陀螺的拉格朗日量, 即动能-势能.
具体表达式很长, 我会略去一些不必要的部分:
首先考虑陀螺的重心速度对应的平动动能, 重心速度由重... 阅读全帖 |
|
H********g
发帖数: 43926 | 36 题目是不是有漏洞啊。ab显然得改变嘛，水都跑了还能不改变么。假设漏水速度一样，
那么在特定时间，x0＝x3>x1>x2。而且x1-x2的值跟水柱的初始方向相关，换句话说跟
孔所在的平面和地平面的夹角余弦成正比。但是假如孔无限小，那这个差值显然也是无
限小。
而且实际上x3应该是略小于x0的，跟水柱的体积有关，但是既然说孔大小可以忽略，水
柱长度有限，那体积趋近于0所以也忽略掉了。假如水柱相当长，那x0>x3。 |
|
x*********s
发帖数: 5554 | 37 我买t410s的时候那个插槽没有选读卡器，现在想利用expresscard 34的插槽安装一个多
合一读卡器，读不知道什么尺寸合适
newegg.com上面有几个，但尺寸分别是
2.3"x1.9"x0.6"
2.95"x1.29"x0.19"
2.95'x1.34"x0.19"
不知道哪个厚度以及长度正好。
请教一下，谢谢 |
|
c****7
发帖数: 1245 | 38 你这不是说东扯西吗
我贴了两个主板做比较，说其中一个太烂不值得买，你要辩的话那就说说那个烂的怎么好
你贴出第三个主板想说明啥？ 有图为证，证个啥？
另外,
你这 MSI z87 g45 主板是有三个PCIE，但是很蛋痛
3xPCIe 3.0 x16 slots
- The operating mode will be x16/x0/x0, x8/x8/x0, or x8/x4/x4. Please always
install expansion card into the PCI_E2 "First".
三个都能用上的mode只能x8/x4/x4, 3.0又怎样, x16的插x4上的话就很影响带宽了
那张技嘉X58顶级板可是x16/x8/x8, 虽然是2.0, 但是给你贴过的那个link里测了
2.0vs3.0几乎无差别, 不多说了. |
|
b***i
发帖数: 3043 | 39 GraphicsContext gc = canvas.getGraphicsContext2D();
以Brick为模板，在initTimeline()里面
KeyFrame kf = new KeyFrame(Config.ANIMATION_TIME, new EventHandler<
ActionEvent>() {
public void handle(ActionEvent event) {
加入
GraphicsContext gc = canvas.getGraphicsContext2D();
int x0=random.nextInt(Config.SCREEN_WIDTH);
int y0=random.nextInt(Config.SCREEN_HEIGHT);
int x1=random.nextInt(Config.SCREEN_WIDTH);
int y1=random.nextInt(Config.SCREEN_HEIGHT);
gc.setStroke(C... 阅读全帖 |
|
f********i
发帖数: 563 | 40 现有一个问题，算是集合领域的一个应用：
Set1 = { u | u = a0*x0 + a1*x1 + a2*x2 + a3*x3 }
Set2 = { v | v = b0*y0 + b1*y1 + b2*y2 + b3*y3 }
其中 ai 和 bi是常系数，已知。 而 x 和 y 是未知数（整数），但范围已知。想知道
Set2是否是Set1的子集。
注： 未知数x的个数可以是arbitary的，但是x和y的个数必须一致，如上面的4个。也
可以是2个，如下面的例子，但在解题的时候，其个数以及各自的系数，都是已知的。
举个例子：
Set1 = { u | u = 8*x0 + 1*x1 } 其中： 0 <= x0 <= 3, 0 <= x1 <= 3
Set2 = { v | v = 1*y0 + 8*y1} 其中： 0 <= y0 <= 2, 0<= y1 <= 1
经计算： Set1 = { 0, 1, 2, 3, 8, 9, 10, 11, 16, 17, 18, 19, 24, 25, 26, 27 }
而 Set2 = { 0, 1, 2, 8, 9, 10 }
故而 S |
|
b***y
发帖数: 2799 | 41 ☆─────────────────────────────────────☆
observer (笑看人生) 于 (Tue Jun 12 13:04:21 2007) 提到:
N * N的 点阵，点有黑有白，找一个矩形，四角点全黑，面积最大
好象还没人给解法，继续攒人品。
1. Scan 矩阵，找黑点，记录坐标(xi,yi)，存到array， O(N^2)
得到 M 个黑点
2.a 土办法是scan 黑点array, 对任意2个点(x1, y1), (x2, y2)
如果(x1, y2), (x2, y1)也是黑点，（直接查原点阵）
那么是一个合格矩形,
找出max area O(M^2)
2.b 另外一个办法，根据x,y坐标，sort 黑点array，
得到2个sorted array, AX, AY, O(MlgM)
scan AX, 对2个x座标相同的黑点，(x0, y1), (x0, y2)
查 AY找y1, y2, 在y1, y2的黑点中找x座标相同的点,
(x3,y1), (x3,y2), x3 != x0, 那么就是一个合格的矩形
找到max area
一般情 |
|
b***y
发帖数: 2799 | 42 ☆─────────────────────────────────────☆
foolboylei (老高卖菜) 于 (Thu Jun 19 13:18:01 2008) 提到:
现有一个问题，算是集合领域的一个应用：
Set1 = { u | u = a0*x0 + a1*x1 + a2*x2 + a3*x3 }
Set2 = { v | v = b0*y0 + b1*y1 + b2*y2 + b3*y3 }
其中 ai 和 bi是常系数，已知。 而 x 和 y 是未知数（整数），但范围已知。想知道
Set2是否是Set1的子集。
注： 未知数x的个数可以是arbitary的，但是x和y的个数必须一致，如上面的4个。也
可以是2个，如下面的例子，但在解题的时候，其个数以及各自的系数，都是已知的。
举个例子：
Set1 = { u | u = 8*x0 + 1*x1 } 其中： 0 <= x0 <= 3, 0 <= x1 <= 3
Set2 = { v | v = 1*y0 + 8*y1} 其中： 0 <= y0 <= 2, 0<= y1 <= |
|
O*******d
发帖数: 20343 | 43 说一下我怎么做的antialiase。 先把图形计算好，比最后的图像多出宽三列高三行。
图形是一个浮点数的点阵。 两个相邻的点之间的距离是一个单位。 然后用三次方程插
值。每个插值要用图形中4X4的值来计算。 在四个点组成的方格中均匀插入5X5=25个点
（单核机）或7X7=49个点（多核机），给每个插值赋予颜色，最后把一个方格中所有的
颜色平均，就是一个像素的颜色。
由于一组插值用的是同一套三次方程系数，所以写了一个class来做。 其实是一个
template。 做插值时，有两个套着的循环，分别在X方向移动和Y方向移动。 内循环是
在Y方向从上往下移动，这样每次移动一步，只需要计算最后一列的三次方程系数。下
边是我写的Cubic Interpolator. 是一个recursive template。For bicubic
interpolation, N = 2. 所有的重复计算都尽可能的避免了。 这个templete的计算
速度，在N==2时，是单独bicubic interpolation的速度的三倍。
template
class Cu... 阅读全帖 |
|
d****i
发帖数: 4809 | 44 Again, 你没有得到它
At t=t0, x = x0 (asm: MOV R0, #x0)
At t=t0+delta_t, x = x0+1 (asm: ADD R0, 1)
So in real computer programming languages such as C/C++/Java, it is
expressed as
x=x+1 (or simply x++)
哪里有什么问题？ |
|
g******r
发帖数: 139 | 45 谢谢你的建议。我知道qPCR的efficiency, 但确实没有测过各个primer的efficiency到
底怎么样。但理论上effiency对Absolute Quantification影响很大，而对ddCT方法的
影响有限。而我就只用ddCT方法比较各个gene在不同条件下的表达，没有去比较基因间
的表达。至于ddCt的方法有人要argue有缺陷，不理想，那是另外的话题 (体外话，我
发现没有测primer效率的人不在少数，当然绝大部分人也是用ddCt)。
以一个简单的knockdown为例，设目标gene表达量为X0，在Ct时的PCR产物为Xt0，效率
为EX, reference gene的表达量为R0，在Ct时的PCR产物为Rt0，效率为ER；在
knockdown情况下目标gene表达量为X1，在Ct时的PCR产物为Xt1, reference gene的表
达量不变，仍然是R0，在Ct时的PCR产物为Rt1 (两次的Ct可能一样，也可能不一样)。
Xt0=X0 (1+EX)CtX0, Rt0= R0 (1+ER)CtR0 ;Xt1=X1 (1+EX)CtX1, Rt1... 阅读全帖 |
|
w**l
发帖数: 5 | 46 I am surprised that these CS guys are so ignorant of Math. If my understanding
is correct, this is almost a high shool level problem. If his ellipse means
the quadratic surface, then just solve the joint equations:
1. (x-x0)'A(x-x0)=1
2. x=x1+t*(x2-x1)
where A is a square matrix, x0 is the center of the ellipse, x1 and x2 are two
points on the line. Plug 2 into 1, it remain to check whether the following
eqaution about t has solution:
[x1+t*(x2-x1)]'A[x1+t*(x2-x1)]=1
which is a quadratic functio |
|
f********r
发帖数: 50 | 47 问题解决了，多谢rossby的提醒
解法如下（matlab)
function x=normalize(A,B)
x0=eye(size(A,1));
function F=myfun(x)
F=A*x*x'*A'+B*x*x'*B'-eye(size(A,1));
end
x=fsolve(@myfun,x0);
end
只是似乎效率不是特别高，算一个4阶的矩阵有时需要0.5秒
有没有直接求解的方法呢？
还有，似乎不是所有的矩阵A,B都有解，有时解又依赖于初始的x0,
是不是不是所有的矩阵A,B都存在解呢？ |
|
c********i
发帖数: 55 | 48 好冷清……
大概一想，估计个a的初始值，从 ydot = ( y(x1) - y(x0) ) / ( x1 - x0 ),
x1 = exp(-8), x0 = 0
我不知道MATLAB的ODE45到底怎么用，不过楼主可以自己整一个RK4的程序（比较简单）
，然后自己从0积分到1，比较积分误差，给a一个纠正值，再从头积，…… 到满意为止
当然，要是能对a有一个估计最好 要不然可能是超级的慢喽 |
|
j**u
发帖数: 6059 | 49 not sure what you need exactly. the following script generates a figure like
what you showed above.
x0=-10:.1:10;
[X,Y]=meshgrid(x0,x0);
Z=exp(-X.^2-Y.^2);
Z(Z>0.3)=0;Z(Z<0.1)=0;
contour(Z,8); |
|
