c***u
发帖数: 190 | 1 就是5000个linear 线型方程组,I=RxD, R是5000x5000矩阵, 局部会有zero, 知道I和D
(实验上得到的).5000x5000计算对电脑要求比较高吧? 这个有效计算就是仪器采集大量
的数据,然后经过矩阵运算, 转换出我们需要的数据..
哥们能个我一个matlab的简单code吗? 我们可以先用random生成的数试.
谢谢 |
|
c***u
发帖数: 190 | 2 现在做个简单的数值计算, 老板需要知道怎样能很快,有效的计算5000x5000的矩阵. 我
本人不是很懂科学计算. 也不是很懂matlab 和mathematics. 有人做过类似的活吗?
谢谢 |
|
t***s
发帖数: 4666 | 3 what's the dimention of I and D? left divide a 5000x1 vector
by a 5000x5000 matrix in matlab on my machine only takes ~40 seconds. |
|
T**********t
发帖数: 449 | 4 基本假设：
（1）班次信息用点到点方式，5000个火车站，单个信息不到100字节，
5000X5000 = 25 M X 100 = 2500 M = 2.5G，这可以全部装入内存，不用考虑硬盘表现
了。
（2）主要流量来自于对当天出来的票源的抢购。每天新增票源按照每两个连续站就算
一票，来存储
6000 班次 X 30站 /每班次= 18000区间 X 14车厢 X 5种位置 = 1.260M
（3）每小时查询峰值：5000万次
（4）每小时交易峰值：200万次， 每秒交易量： 555次
下面是基本架构：核心数据库、剩余票源发布服务器、中心城市订票端服务器
其中
1、核心数据库：完成真正的交易，更新现有座位数目。每秒预计交易量： 555次。
Capability of 1 typical Oracle Database: 30,000 TPS for credit card
transaction
Capability of 1 IBM Power Database: 50,000 TPS for credit card transaction
结论：数据库不是bottleneck，... 阅读全帖 |
|
a*********s
发帖数: 122 | 5 雖然知道版上萬佛都是用神掃9K ED
但是還是想問問有沒有人知道A&I的掃描品質如何?
所謂的high resolution for 6x6檔案大小為67mb
這樣大概多少pixel呢? 能有將近5000x5000嗎
謝謝 |
|
l*******n
发帖数: 155 | 6 我在word里画的图，准备copy到photoshop里,每次我在photoshop里新建一个文件的时
候，photoshop会自动根据剪贴板里像素来创建画布的大小。我遇到的情况挺奇怪，同
一个图，copy过去后，像素的大小去不一样（分辨率都一样是72）。比如有时是
500x500，有时却是5000x5000。反复搞了好多次，就是搞不懂。哪位高人有经验给说说
，谢谢啦。 |
|
发帖数: 1 | 7 做一回雷锋
如果你想从raw fastq data到结果，大概两个步骤：
1. raw fastq--bwa--mutation/gene expression
主要是read mapping，把GB level的原始数据BWA mapping得到bam file，然后如果你
的目标是找variants那就用GATK这样的主流软件
1a. 对于你来说，这个步骤里你主要需要学会linux environment，bash，学会用HPCC
来submit job，学会tune BWA/GATK的参数，是很容易的。
1b. 同时你要学会面对big data，都是比如300 million行的数据，处理一个全基因组
需要1TB空间，这个是很可怕的。大数据的storage，transfer都要注意。
1c. 什么BWA/GATK这种复杂算法高级原理C/JAVA的你不需要care
你可以理解这一步就是从海量海量的数据里初步filter出你要的东西，但是粗糙的东西
，不是完整产品
2. gene expression---统计分析/美丽的图图
这个主要是RNA-seq，你tophat之类的得到一堆基... 阅读全帖 |
|
t***s
发帖数: 4666 | 8 what do you mean by "compute effectively"? |
|
c******m
发帖数: 599 | 9 5k*5k
如果是稀疏矩阵, matlab就搞定了 |
|
|
|
x****g
发帖数: 6597 | 12 这个得要Matlab， 我搞过6000x6000的矩阵.
其他软件内存不够。 |
|
|
y*****a
发帖数: 580 | 14 If spares matrix, not a problem at all (particularly n-diagonal matrix with
smaller bindwidth.), and I think you can try D=I\R in matlab. The problem
of using a random generator is that, you may need pre-conditioning of your
matrix. If you want to code, get a "Numerical Recipes" book to look at Chapter 2. There also have a code in the book to solve for sparse matrix using conjugate gradient method.
It would be better to identify the structure of your matrix, the size itself is not a problem ( |
|
h******g
发帖数: 69 | 15 At octave 3.2.0 Windows XP
b=rand(5000,1);
A=rand(5000,5000);
tic(); x=A\b; toc()
only takes 30 seconds.
Computer is Xeon CPU 2.8G, 2G RAM. |
|
|
p*****e
发帖数: 310 | 17 matlab就是用lapack算的，不用舍近求远 |
|
m**********o
发帖数: 19 | 18 Use an iterative method such as conjugate gradient or inverse with singular
value decomposition. |
|
k*********g
发帖数: 791 | 19 用迭代法；
用gauss－seidel；
不过，每行得先找到绝对值最大的，再gs； |
|
