o****n
发帖数: 9475 | 1 突然想起这位室友，因为刚毕业时一直认为，到09年就会有她的消息，虽然两年前就有
了消息，最近还是老想起她
我们学校每年都招中保生，就是从中专报送上来的，我们那年有四个，室友是其中之一
。她是初中中专，读了四年还是五年的那种，她也是四个中保生里唯一一个成绩好的，
另外三个每学期都补考。中保生除了和我们一样的学费，还有另外付培养费，这钱他们
原来的学校出，他们和学校也有合同，毕业后必须回学校任教多少多少年，要不要付几
万块的违约金。
我们班七个女生，六个北方的，我们住一个宿舍，唯一一个南方的住隔壁，她叫这位室
友阿如，我们就都跟着这样叫。阿如老家在中原，父母支边，她上的是边区的中专。按
规定，毕业后回边区。
大一刚开学那会儿，我们其他几个玩疯，那有活动往那噌，除了上课，其他时间都在玩
，还好没逃课。阿如从来不和我们一起，她没上过高中，那时她不知道什么是摩尔，什
么是动量守恒，什么是矢量、标量，高数，大学物理，普通化学这些课对她无疑是天书
，还有英语，她说的，中专时也上英语课，基本是玩的。我们疯的时候她在自学高中的
数理化，补英语。前两年我们基本没和她一起活动过，上课以外见她一面都不容易。平 |
|
s**********l
发帖数: 8966 | 2 [心得]图形和图像
发信站: 水木社区 (Thu Dec 18 21:45:24 2014), 站内
既然这个板块是讲这个两个领域，恰巧学在图形（Graph），干在图像（Image），聊聊
两者的异同。一家之言，想到哪儿，写到哪儿，还望指正。
图形是以点线面等3D的元素来描述这个世界，大多数的时候是以矢量的形式来呈现。而
图像是是以2D的一个一个像素来真实的呈现这个世界，多数时候是标量的单个像素的矩
阵表示。图像加上时间维度就成了视频。图像是人眼对于世界最原始的还原，而图形则
是人类对世界的抽象。通俗点，如果要描述一张桌子，拿着相机咔嚓一下就是一张图像
，而图形则要用点线面来构建一个模型，然后加上材质光照来通过一定的模型来渲染而
变成一个图像。
图形的主要内容是建模（Modeling）和绘制（Rendering），建模的代表软件是建筑类
的CAD类的软件，譬如AutoCAD，比较通用的3DMax，Maya等。绘制最著名的软件是
RenderMan，专门用来制作电影的。图形学现在的发展方向是模型编辑，大场景建模，
真实感绘制，非真实感绘制，群体动画，运动捕获，基于物理的仿真等... 阅读全帖 |
|
w***g
发帖数: 5958 | 3 用内积算比较直观。要判断x和Y=[a,b,...]是否独立只要把x逐个投影到Y的每个元素，
并把投影的那部分从x中减去。如果一圈下来x还有剩的那就是线性独立了:
foreach y in Y {
x -= (x . y) / (y . y) y
------- -------
内积 内积
----------------- ---
标量 矢量
}
然后判断x是否为零就行。 |
|
t****t
发帖数: 6806 | 4 matlab的矩阵计算核心应该本来就是公共的库, blas, lapack之类的, 不过matlab都给
你包装好了, 优化都做好, 比如说在intel的CPU上就会用MKL, 所以矩阵规模越大越合
算, 如果都是标量就慢了. |
|
c*m
发帖数: 1599 | 5 申威基于当年的alpha处理器，还行。
申威5效率和功耗都比NVidia Titan GPU强，还是CPU，已经不错了。
只是申威的单核性能确实不济（虽然性能每瓦不错），做个人使用的桌面电脑确实不行。
放君搞得几款处理器都还行，虽然借（tou）鉴（dao）了别人的设计，但是搞出来的
CPU不算差。例如郭芳科大搞得Sparc指令集CPU和ARMv8指令集CPU,都还可以。
郭芳科大最早是逆向山寨Intel的IA64体系处理器安腾，只是安腾CPU在市场上举步维艰
，这个逆向工程虽然成功了却不得不放弃（立项时迷信intel，没想到intel也有失误的
时候）。然后因为太阳微系统免费开放了Sparc的设计代码，于是郭芳科大就在Sun的设
计基础上改改，弄出了飞腾1500，用在天河2号巨星机的通信上。 （Sparc的特点是通
信很牛但是浮点性能不足。）天河2的主要计算节点还是靠intel的CPU。因为Sparc指令
集市小，所以放君又瞄向了ARM指令集。郭芳科大ARMv8指令集，四发射乱序执行，拥有
多达64个核心，主频达到2GHZ的服务器CPU的理论浮点计算能力是512G DP Flops，采... 阅读全帖 |
|
T********i
发帖数: 2416 | 6 超级小型机。
俺爹70年代初期参与的海军项目，每秒百万次，指令数据cache，超标量长流水。最后
也做出来了。
架构上，早有人想出来了。20年后，Intel的80486也是这些东西。关键的形成工业的能
力。 |
|
发帖数: 1 | 7 这里面讲了，为什么goroutine太多，在多路众核处理器上会大量造成cache-line miss
和ping-pong效应。我们目前最强的机器是2路200+核心16个内存控制器，还没有测试，
估计golang性能因为scheduling会很难看。做处理器的最怕跨路访问内存，一个load应
该有200~300ns延时，再加上超标量，相当于等待成百上千条指令啊，性能下降10倍算
少的。
The main concepts in Go runtime are:
M - OS thread (Machine).
P - Logical CPU (Processor), there are exactly GOMAXPROCS P's.
G - Goroutine.
Netpoll - network poller (epoll descriptor).
RunQ - scheduler queue with runnable G's.
MHeap - global memory allocator state (contains central caches and free
spans)... 阅读全帖 |
|
发帖数: 1 | 8 这里面讲了，为什么goroutine太多，在多路众核处理器上会大量造成cache-line miss
和ping-pong效应。我们目前最强的机器是2路200+核心16个内存控制器，还没有测试，
估计golang性能因为scheduling会很难看。做处理器的最怕跨路访问内存，一个load应
该有200~300ns延时，再加上超标量，相当于等待成百上千条指令啊，性能下降10倍算
少的。
The main concepts in Go runtime are:
M - OS thread (Machine).
P - Logical CPU (Processor), there are exactly GOMAXPROCS P's.
G - Goroutine.
Netpoll - network poller (epoll descriptor).
RunQ - scheduler queue with runnable G's.
MHeap - global memory allocator state (contains central caches and free
spans)... 阅读全帖 |
|
b******d
发帖数: 28 | 9 【 以下文字转载自 Programming 讨论区 】
【 原文由 Bluewind 所发表 】
1. 在emacs下编译c程序时,我的emacs总是只从~/下
找makefile,如何将其设成自动从当前buffer中文件所在的路径下找? 如:
输入M-x compile, emacs提示
compile command: make -k
回车后出现compilation buffer,里面提示
cd /home/myaccount/mydirectory/
make -k
make: *** No targets specified and no makefile found. Stop.
这时makefile就在mydirectory下,可是emacs找不到,
把makefile易到我的~/下就找到了
该怎么解决?多谢
2. perl中的正则表达式
现有一标量$a,我要替换其中的某个词,然后保存到一个新的变量$b,
但我还需要原来的$a,我可以这样做
# for example, $a = "hello, world";
$b = $a;
$b =~ s/hello/hell/; |
|
f*******d
发帖数: 339 | 10 译的不太准确， scalar field 应译为标量场， dark energy 一般译为暗能量。
唉， 两年前我就知道这个塌缩的结果， 还和当时的老板讨论过，
不过没想到要深入研究下去写成文章， 因为当时觉得没太大意义，而且以前也有人写过文章。
现在想想亏了， 后来Steinhardt, Ovrut, Turok, Linde 这些大腕往string， Supergravity 上一扯，
硬是把这玩意炒火了。 |
|
c****u
发帖数: 584 | 11 你真能唐啊。
这些球差，慧差，场曲。。。不是 Abbe limit考虑的东西了。严格的来讲，
aberration的衍射过程要上有限元才能得矢量解的。标量解对大NA系统事很不准，更别
说只考虑折射率了。
当然这些末节不会入真正学物理的法眼的。 |
|
c****u
发帖数: 584 | 12 俺设计过物镜，一但NA 靠近全反射，Ray Tracing 就歇菜了。CodeV不知道， Zemax用
标量仰射＋Monte Carlo来模拟残不忍睹. |
|
c****u
发帖数: 584 | 13 你真能唐啊。
这些球差，慧差，场曲。。。不是 Abbe limit考虑的东西了。严格的来讲，
aberration的衍射过程要上有限元才能得矢量解的。标量解对大NA系统事很不准，更别
说只考虑折射率了。
当然这些末节不会入真正学物理的法眼的。 |
|
c****u
发帖数: 584 | 14 俺设计过物镜，一但NA 靠近全反射，Ray Tracing 就歇菜了。CodeV不知道， Zemax用
标量仰射＋Monte Carlo来模拟残不忍睹. |
|
n******e
发帖数: 110 | 15 研究一个蛋白在正常组织和癌症组织的表达量。找公司买的病人匹配的样品。问题来了
，总的上样蛋白量根据公司的浓度是一样的，考马斯染显示两个样品不同条带的浓度都
不一样，但总量估计是差不多，符合公司的标量。但是用actin做，正常样品估计只有
肿瘤样品actin的十分之一。问题来了，怎么比较这个蛋白的含量？
很多文献里都用一个house keeping基因的蛋白做loading control。我可以用蛋白总量
做control吗？被接受程度都差不多？
谢谢！ |
|
n******e
发帖数: 110 | 16 研究一个蛋白在正常组织和癌症组织的表达量。找公司买的病人匹配的样品。问题来了
，总的上样蛋白量根据公司的浓度是一样的，考马斯染显示两个样品不同条带的浓度都
不一样，但总量估计是差不多，符合公司的标量。但是用actin做，正常样品估计只有
肿瘤样品actin的十分之一。问题来了，怎么比较这个蛋白的含量？
很多文献里都用一个house keeping基因的蛋白做loading control。我可以用蛋白总量
做control吗？被接受程度都差不多？
谢谢！ |
|
c******x
发帖数: 67 | 17 单参考方法经常被看作是黑箱工具。对于比较轻的主族元素构成的体系，这通常不会有
什么问题。但是对于重元素体系，即便是最简单的标量相对论计算，情况就复杂得多。
研究工作者至少对所用的量化模型，相对论哈密顿量，基组，电子组态等有所了解，否
则连自己做了什么都不知道。
这里，我分析一篇Chem. Phys. Lett.杂志今年年初发表的文章：
Tao Su, Chuan-Lu Yang, Xin-Qiang Wang, Feng-Juan Bai, and Mei-Shan Wang,
Theoretical characters of the ground states of YbX (X = F, Cl, Br, I, At),
Chemical Physics Letters, 467(4-6), 265-269, 2009.
在这篇文章中，作者犯了很多错误。
1. 对所用基组不熟悉。例如，作者称氟原子用的基组是aug-cc-pVQZ，但基组收缩形式
却是cc-pVQZ；Yb原子所用的ECP使用了错误的名称（ECP28MDF）和错误的文献（Ref.
24-25）。文献中对Yb原子没有发表 |
|
w****j
发帖数: 237 | 18 【 以下文字转载自 Mathematics 讨论区 】
发信人: winglj (coolking), 信区: Mathematics
标 题: 一个简单的矩阵函数问题
发信站: BBS 未名空间站 (Thu Feb 22 18:59:53 2007)
想请教各位一个简单的矩阵函数问题：
一个 n-by-n 矩阵,
其本征值（eigen value)为 e1,e2,...,en,
本征量矩阵E=diag(e1,e2,...,2n),
本征矢量矩阵为 T
A=T*E*T^(-1)
一个标量函数f(x),例如 exp(x)
那矩阵函数 f(A) 是不是等于
如果本征值不同:
f（A)=T*f(E)*T^(-1)=T*diag(f(e1),f(e2),...,f(en))*T^(-1)
如果有本征值相同:
f(a)=T*[f(e1),f'(e1)/1!,f''(e1)/2!,....;0,f(e1),...]*T^(-1)
好像隐隐约约记得是这样定义的，手头没有工具书，现谢谢各位了 |
|
l*2
发帖数: 21 | 19 请教MATLAB高手们：怎么样改写下面的代码，把循环体去掉。或者有什么其他优化的方
法。感觉这个循环太慢了。
其中 A,B,C,D都是固定的二维矩阵, f是标量。一个二维矩阵点乘另一个矢量的各个元
素，有没有办法不用循环体？
for f_idx = 1:length(freq)
f = freq(f_idx);
K = eye(5) + j*2*pi*f*A;
M = inv(K);
N = det(K);
S = exp( j*2*pi*f * D - 2*pi* f^2 * C*M*B ) / sqrt(N));
end
多谢！ |
|
l*2
发帖数: 21 | 20 这个实现什么无所谓，我想问的就是像这样的二维矩阵乘标量能不能避免使用循环体。
像下面这个简单的程序：
for f_idx = 1:length(freq)
f = freq(f_idx); % f是矢量freq的元素
S ＝ A*f; % A是二维矩阵。至于S怎么处理不用关心
end |
|
s******h
发帖数: 54 | 21 看上去应该没有区别的，但是在程序里跑出来结果就是不一样。。。。
这些数都出现在矩阵元里面的，然后矩阵又经过若干的操作，最后的结果是一个标量，
类似于求了矩阵的迹（trace）。。。
有没有达人知道这究竟是怎么回事啊。。。。。 |
|
h**********c
发帖数: 4120 | 22 用计算几何二分法，
input两个有向线段，
output 二维点坐标
固定一个，
0.初始一个标量t
1.判断另外一个线段上一点，在上一个线段的左侧还是右侧（一般的cg书上都有介绍）
2。沿左或右相反方向延长加倍或减半t，算出新点。
重复1，2直至找到收敛点
线性问题不解方程最好的办法就是二分。 |
|
I**********8
发帖数: 292 | 23 x^HA_1x=a_1
x^HA_2x=a_2
...........
x^HA_Nx=a_N
其中A_i,a_i是常数矩阵（非正定）和标量。 x是未知向量。
多谢 |
|
|
|
l***x
发帖数: 21 | 26 it's going to be a Jacobian matrix. please check matrix calculus. |
|
s*****r
发帖数: 69 | 27 我觉得还是矢量，先做分解，在各分量上分别求导，然后合成 |
|
t**s
发帖数: 52 | 28 d(row vector)/d(column vector) = matrix
d(dolumn vector)/d(row vector) = matrix
d(row vector)/d(row vector) = row vector
d(column vector)/d(column vector) = column vector
详细解释参见此书 A Matrix Handbook for Statisticians, by G. Seber. |
|
d*******d
发帖数: 3382 | 29 哪里来的问题？
如果rho是一个标量数字，难道不是A=B*rho? |
|
x*****d
发帖数: 427 | 30
弦论通俗演义
李 淼
第七章 先声
(第四节)
我们前面介绍的T对偶，既可以用在玻色弦理论中，也可以用在超
弦理论中。用于玻色弦时，情况很简单，无非由一个玻色弦得到另
一个玻色弦；用于超弦时，情况稍复杂，在T对偶下，IIA理论变成
IIB理论，反之亦然。这个现象有一个简单的世界面上的解释。在世
界上，当我们做T对偶时，是将动量模与绕量模互换，这个互换可
以通过改变世界面上对应的标量场 (即紧化的那个空间) 的左手模的
符号达到。由於要保持世界面上的超对称，对应的世界面上的费米
子的左手模也要改变符号。我们知道，时空中的费米子 |
|
x*****d
发帖数: 427 | 31 弦论通俗演义
第十章 第二次革命：D膜
第三节 非阿贝尔对称性
李淼
膜物理学带来最有意思的物理是当有若干个膜同时存在的时候，特
别是若干个平行的维度相等的D膜。威顿在95年十月份泡耳钦斯基的
文章出现在网络的两周后，写出了一篇研究多D膜和D膜束缚态的文
章。这篇文章第一次指出当多D膜存在时，会出现非阿贝尔对称性。
前面已经看到，当两个D膜同时存在时，除了两个端点都搭在同一个
D膜上的开弦，还有两个端点搭在不同D膜上的开弦。前者组成两个D
膜上的低能激发态，包括每个D膜上的无质量规范场和标量场。在超
弦理论中，它们形成两个阿贝尔矢量超多重态。现在的问题是，两个
搭在在不同D膜上的开弦是什么样的激发态？
其实，这个问题在弦论的早期就有了回答。现在我们用1来标志搭在
第一个D膜上的端点，用2来标志搭在第二个D膜上的端点。如果开弦
是可定向的（在型II弦论中必须是可定向的，原因是两个开弦可以合
并成一个闭弦离开D膜，从而是型II理论中的一个闭弦，必须是可定
向的），用1、2来标志端点还简单了一点，因为定向性本身要求在一
个D膜上有两种不同的端点，一种对应于弦的一个定向，其实相对于
该D膜上的 |
|
x*****d
发帖数: 427 | 32
第十二章 黑洞
第四节
李淼
黑洞的膜的理论到了马德西纳和斯特劳明格的灰体谱计算可以说到了高
潮。
我们在上一节提到达斯和马图关于D膜辐射的计算，这个计算假定黑洞
离临界黑洞不远，这样除了标准黑体谱中应有的普朗克公式外，还有一
个和吸收截面成正比的因子。这个因子的存在是因为被辐射出的无质量
粒子要克服黑洞的引力，幅射谱不再是简单的黑体谱。在黑洞的半经典
计算中，计算这个因子比较复杂。而D膜的图象中，我们只要计算两个
开弦合并成一个闭弦的几率，计算相对简单。这样简单的计算所给出的
结果与半经典计算完全一样，不能不令人惊讶。在离临界黑洞不远时，
对于一个被辐射出的标准标量粒子来说，当能量低于霍金温度时，这个
吸收因子与黑洞视界的面积成正比，也就是说，对于长波来说，黑洞看
起来象一个黑盘。
而马德西纳和斯特劳明格则考虑黑洞离临界极限比较远的情况，同样考
虑低能极限。在这种情况下，吸收因子是能量的一个复杂函数，不仅仅
与黑洞视界面积有关。这个因子破坏了标准的黑体谱，所以这个因子被
称为灰体因子。同样，吸收因子在黑洞背景下的半经典计算很复杂，因
为黑洞的几何比较复杂。但是，他们的计算结果表明 |
|
w****j
发帖数: 237 | 33 想请教各位一个简单的矩阵函数问题：
一个 n-by-n 矩阵,
其本征值（eigen value)为 e1,e2,...,en,
本征量矩阵E=diag(e1,e2,...,2n),
本征矢量矩阵为 T
A=T*E*T^(-1)
一个标量函数f(x),例如 exp(x)
那矩阵函数 f(A) 是不是等于
如果本征值不同:
f（A)=T*f(E)*T^(-1)=T*diag(f(e1),f(e2),...,f(en))*T^(-1)
如果有本征值相同:
f(a)=T*[f(e1),f'(e1)/1!,f''(e1)/2!,....;0,f(e1),...]*T^(-1)
好像隐隐约约记得是这样定义的，手头没有工具书，现谢谢各位了 |
|
x******i
发帖数: 3022 | 34
1
1。x和dx是啥关系？
2。x和dx是标量还是D维向量？ |
|
g******s
发帖数: 733 | 35 不好意思,因为随机数均值为零, dx=x. x是标量, 在给定T的情况下,每次产生一个符合
g(x)分布的随机数,共产生D=32个随机数. |
|
g****t
发帖数: 31659 | 36 Bourbaki的书里面指数函数是微分方程定义的吧。
x'=x定义出exp函数。
标量微分方程解x'=x存在是用Picard iteration方法证明的。
所以推广到x'=Ax我想是比较自然的想法。
但往矩阵推广的严格证明其实还是挺麻烦的，
一般的微分方程课本上很多证明都有漏掉环节。
矩阵是有限维的线性算子，这个幂级数的定义已经
是最本质的了。幂级数还不够 “elementary”吗？
更加不简单的理解也有很多，可以查 “谱分解”，
“算子半群” 或者 “李群指数映射”
spectral decomposition
semigroup of operators
exponential map |
|
D*******a
发帖数: 3688 | 37 我记得当时老师说exp(At)就是用微分方程定义的，是个符号
只不过最后“碰巧”跟标量函数exp(x)的特性一致 |
|
I**********8
发帖数: 292 | 38 【 以下文字转载自 EE 讨论区 】
发信人: IlIlIlIl2008 (kuaichezou), 信区: EE
标 题: 请问这组方程怎么解
发信站: BBS 未名空间站 (Mon May 4 10:05:10 2009)
x^HA_1x=a_1
x^HA_2x=a_2
...........
x^HA_Nx=a_N
其中A_i,a_i是常数矩阵（非正定）和标量。 x是未知向量。
多谢 |
|
x*****d
发帖数: 427 | 39 你这不是胡说八道吗？ 力怎么就跟坐标有关了？ 应力就跟坐标无关了？
1阶张量怎么就是标量了？ |
|
h******g
发帖数: 121 | 40 看书说，2*2*2 rank-3 tensor 和一个三维向量相乘会得到一个标量，不明白怎么乘的
，谢谢。 |
|
