Monthly Archives: September 2011

镜先生谈国人技术升级

现在我连“观察”这个词都觉得恶心了,原因就是看到有的人自称“互联网观察者”,然后所发的文章和微博无非就是什么Android啦Apple啦之类的平媒IT版水平的东西。所以我不能说我“观察”了好多年的国内形势,更不能说我“观察”了好多年的国外形势。那我到底“XX”了国内形势呢?难道说“啥也没观察!活在中国这么多年白活了!”——也不至于吧。但假如说自己是“对中外政经科教文卫形势均有精辟独到的见解”——这个就简真是自暴自弃了,好像反正甩不掉“恶心”的嫌疑,那就干脆恶心到极致的意思。我想到一种说法“对中外政经科教文卫形势的个人臆测”,好像比较中肯。

所以,根据本人对中外政经科教文卫形势的个人臆测,我觉得在中国基础研究跟技术市场的结合还不可能。中国能卖得到国外去的东西都没啥高科技;而光靠卖国内,马上赔死。现在国内的所谓“生物科技”还停留在生产“xx素口服液”的水平,不然为什么曹谊林做科研是做组织工程,搞公司不得不去卖“玻尿酸”呢?连学名“透明质酸”都不敢叫。中国老百姓以前一直是农民,种种地,不需要什么高科技。现在也是西方列强告诉我们要高科技的,那当然是西方列强说什么是高科技那就什么是高科技了呗?要是任由中国人搞,最后保证还是回到阴阳气功本草纲目的年代。

今天跟徐老师讨论到这个问题,怎么说他怎么不明白。可见在国外留学多年确实会对中国的政经科教文卫形势失去精辟而独到的见解的。

恰好又看到科学网李维老师的短文,很值得复制粘贴:

可是世界第一不是吹出来的,需要花力气、动脑筋。要靠这个投入带动国人企业产品和技术的升级。 这就是“折寿”的过程呢。很可能是钱都化了,结果是个“平庸”。因为“带动国人企业产品和技术的升级”这样的事情没有什么“世界第一”的说法,只是个良心和认识的结果。在洋人那里,科技的投入推进企业产品和技术的升级是个当然的事情,不是个问题。而在国人这里,就是个问题了。一个设备,要不要国产?国产的话,要有风险。买洋人的产品,对同样的失败,国人的评价也有所不同。洋人产品坏了的话,国人很富有“同情心”,而同样的事情发生在国人的产品上,那说法就多了。作为当事人,如何选择呢?这是个问题。

买洋人的设备、用洋人的试剂、发论文在洋人的杂志上,可能的话再雇几个洋学生、洋教员。然后就可以很“自豪”地说,俺们也达到世界一流了。写到这里,镜某突然想到了一个曾经发生过的鹿鸣馆的故事。镜某怎么看“清华大学-北京大学生命科学联合中心”怎么象“鹿鸣”生命科学联合中心。现实与历史上的事情往往就是如此地相似。几十年后研究课题的名字都准备好了:21世纪中国组建北清生命科学联合中心与19世纪日本建造鹿鸣馆之间思维脉络的比较。

这才是真正地“对中外政经科教文卫形势均有精辟独到的见解”。

科学的“边角料”是啥?

看到有人在豆瓣上提到一本书叫做《科学的边角余料》,非常直白地把“科普”的本质作为了一本科普书的标题。我之前也试用过一句话是“所谓‘科普’跟‘科’一点关系都没有,它只是力图从科学中挖掘一些人文意味进行发酵”。但这句话太文纠纠了,还不如“边角料”这个词来得纯真粗犷。也许科普教主们会反驳说:

其实人民群众也不求知道高深的科学原理,你以为他们笨到以为这些就是科学本身么?想不到你读了这么多书学到的竟是这种想法。

当然,并不是真的有什么科普教主跟我说过这翻话。我只是猜他们是这么想这么看的。与“人民群众”拉近关系的同时攻击离“人民群众”遥远的其他角色,其实不仅“科普”人懂。反正人民群众也不知道富人真正是啥样,想啥。看也只能看到表面。随便你说。“读过很多书”的人也一样。就跟中小学校曾经告诉他们“科学家很伟大”他们只能信一样,现在的科普人告诉他们“科学很好玩”,他们也只能信。你说这就是“边角料”、那些就不是“边角料”,也是随你说的。你的书就不是用来给能够评判你的受众看的,你爱说啥说啥。

但我也做过相反的事情——提出了“科学史比科学更有趣”这种观点,呼吁这些不得不靠搭乘“科普”快车赚大钱的伪公共知识分子们,反正也是忽悠,搞点“科学旧闻”搭配一下,讲讲故人故事,不也很“人文意味”很“边角料”么?关键是,我估计就算有人这么做了,我看了还是会不满意。因为也许问题不在于“普啥”,而在于今时今日什么形式能“普”什么形式则“普”不成。古典音乐会曾经何时跟现在电影院一样,现在电影院就都差不多变成古典音乐会了……现在就算它讲点儿故人故事吧,都必须要伪装成今人今事;明明记录在案的事情,非要伪装成爆料八卦狗仔队交差作品。其实就算“科学旧闻”,也别指望它诚实——无非还是从中找到一些“边角料”而已,有所言有所不言。

在MATLAB做Fourier变换的相关代码

每次遇到Fourier变换都忘记的基本代码。

如何获得“右半边”的Fourier变换谱?

所谓“右半边”的Fourier变换谱其实就是频率大于零那半边。对FFT来说其实是截取“左半边”。


function [X,freq]=lhsFFT(x,Fs)
% Cuts the left-hand side of FFT
% x is the signal sample to be FTT-ed, Fs is the sampling frequency of the signal
% X is the FFT result, freq is the frequency axis

N = length(x);

freq_bin = 0:N-1; % the frequency bin
freq_res = Fs / N; % the frequency resolution
freq = freq_bin * freq_res; %create the frequency range
X = 2 * fft(x) / N; % normalize the data
cutoff = ceil(N/2); %length of the left-hand side

freq = freq(1:cutOff);
X = X(1:cutOff);

获取某频率的幅值

问题其实就是决定某频率值落在第几个frequency bin,然后四舍五入取最接近的那个频率。


N = length(x);

freq_Res = Fs / N; % the frequency resolution
[X, freq] = lhsFFT(x, Fs) % get the left-hand side FFT

f = 1.234 % The queried frequency is 1.234 Hz for example

freq_bin_prime = round(f/freq_res); % finds the frequency bin closest to f
freq_prime = freq_bin * freq_res; % the resolved frequency closest to f
X_prime = X(freq_bin_prime + 1); % the amplitude of the resolved frequency (vector index starts from 1)