d**g
发帖数: 1031 | 1 网上看来的, 我擦, 都不太会做.已经没法子跳槽啦.
1、简单描述一个单片机系统的主要组成模块，并说明各模块之间的数据流流向和控制
流
流向。简述单片机应用系统的设计原则。（仕兰微面试题目）
2、画出8031与2716（2K*8ROM）的连线图，要求采用三-八译码器，8031的P2.5,P2.4和
P2.3参加译码，基本地址范围为3000H-3FFFH。该2716有没有重叠地址？根据是什么？
若
有，则写出每片2716的重叠地址范围。（仕兰微面试题目）
3、用8051设计一个带一个8*16键盘加驱动八个数码管（共阳）的原理图。（仕兰微面
试
题目）
4、PCI总线的含义是什么？PCI总线的主要特点是什么？ （仕兰微面试题目）
5、中断的概念？简述中断的过程。（仕兰微面试题目）
6、如单片机中断几个/类型，编中断程序注意什么问题；（未知）
7、要用一个开环脉冲调速系统来控制直流电动机的转速，程序由8051完成。简单原理
如
下：由P3.4输出脉冲的占空比来控制转速，占空比越大，转速越快；而占空比由K7-K0
八
个开关来设置，直接与P1口相连（开关拨到下方时为"... 阅读全帖 |
|
h****a
发帖数: 580 | 2 向各位牛人请教一个问题：
现有一套数据（x[m],y[n],z[m,n])，其中x[m]和y[n]是参数，z[m,n]是（x[m],y[n]）
相对应的值，m从1到M，n从1到N，一共有M*N个数据点.
根据我的观测和初步检验，我认为下面的方程可以很好地拟合(fit)这些数据：
z=f(x,y)=a11*x^2*y^2+a12*x^2*y+a13*x^2+a21*x*y^2+a22*x*y+a23*x+a31*y^2+a32*y
+a33
现在的问题是，如何准确地得到方程中的系数coeffients（a11,a12,a13,a21,a22,a23,
a31,a32,a33)?
据我所知，对于二维数据（x[i],y[i]), 很容易用方程y=f(x)来拟合，有很多工具（
gnuplot，xmgrace）可以做到。但是对于三维数据（x[m],y[n],z[m,n]),请问有什么现
成的工具可以用来拟合？若是没有现成的工具，我需要自己写程序，该如何着手？
我在网上搜索“least-squares fitting”，找到一些说明和子程序，但全都是用y=f(x
)来拟合数据(x,y)，没有用z= |
|
|
t*s
发帖数: 1504 | 4 用那个啥Rayleigh...忘了
1000000*
最 |
|
a**a
发帖数: 416 | 5 en, 好像有稀疏矩阵的特征值迭代解法。去netlib上翻翻故纸堆吧。算法是GRMES类的。 |
|
h****a
发帖数: 580 | 6 向各位牛人请教一个问题：
现有一套数据（x[m],y[n],z[m,n])，其中x[m]和y[n]是参数，z[m,n]是（x[m],y[n]）
相对应的值，m从1到M，n从1到N，一共有M*N个数据点.
根据我的观测和初步检验，我认为下面的方程可以很好地拟合(fit)这些数据：
z=f(x,y)=a11*x^2*y^2+a12*x^2*y+a13*x^2+a21*x*y^2+a22*x*y+a23*x+a31*y^2+a32*y
+a33
现在的问题是，如何准确地得到方程中的系数coeffients（a11,a12,a13,a21,a22,a23,
a31,a32,a33)?
据我所知，对于二维数据（x[i],y[i]), 很容易用方程y=f(x)来拟合，有很多工具（
gnuplot，xmgrace）可以做到。但是对于三维数据（x[m],y[n],z[m,n]),请问有什么现
成的工具可以用来拟合？若是没有现成的工具，我需要自己写程序，该如何着手？
我在网上搜索“least-squares fitting”，找到一些说明和子程序，但全都是用y=f(x
)来拟合数据(x,y)，没有用z= |
|
s*i
发帖数: 5025 | 7 【 以下文字转载自 THU 讨论区,原文如下 】
发信人: XiaoP (小屁), 信区: THU
标 题: [转载] 菜鸟问个有限元的弱问题
发信站: Unknown Space - 未名空间 (Mon Sep 22 14:18:30 2003), 转信
【 以下文字转载自 Computation 讨论区 】
【 原文由 XiaoP 所发表 】
我们现在的project想用有限元模拟材料neck以后的样子,我自己做怎么也不收敛,如果不
加入塑性性能的话,abaqus就纯粹把样品当弹性棒处理不出现necking
但我知道肯定是能做的,我看70年代的paper人家就做出来了,他们那时是遍的程序,我想
30多年了,应该有现成的子程序可以用,但不知道去哪里找.
呵呵,不好意思,比较弱,abaqus也是刚学了没多久,老板又逼得急,各位大侠们请帮帮忙,
千恩万谢 |
|
v****e
发帖数: 41 | 8 我查了一下，Nastran可以算。但是我们这儿只有ABAQUS。 用Mises没有问题，主要是要
求用Tresca准则。
UMAT是ABAQUS的用户材料子程序。 |
|
f**********g
发帖数: 117 | 9 怎样用命令来编译多个源程序.f（一个主程序+多个子程序）
好久不用，都忘了。
谢谢。
= |
|
s***h
发帖数: 592 | 10 FEA做app的越来越多，有CAD化的趋势，做基础理论的还好，比如本构之类的
做本构的大牛一个子程序的版权就能卖出天价 |
|
s***h
发帖数: 592 | 11 各人有各人的理解，你有什么想法就说出来，
交流用不着拿门外门内的来压人，要比资历就更没劲了，
非要搬黄历的话，我的umat和hypela2子程序三年以前就卖给公司了，
手。 |
|
m****n
发帖数: 293 | 12 不是大牛。 只说自己的理解：
作FEA用到C++的地方还挺多的， 有C++写的FEA， 也有跟商业FEA软件结合的子程序。
我就用C++ call ABAQUS API 来实现自己的算法。
当然如果就是应用商业FEA软件的话， 会它们的操作和命令就可以了。
如果有误， 望大牛指正。 |
|
h********2
发帖数: 33 | 13 大家好,小弟是中西部一所公立学校的ME博士毕业。背景主要是有限元，平时商业软件
ABAQUS, ANSYS, HYPERMESH, LS-DYNA用的多一点，博士期间主要做的是复合材料，主
要是偏向于工业应用，编程的东西比较少，有项目也写一点FORTRAN子程序,对CFD也稍
微了解一点。现在是拿的OPT EAD, 今年八月份OPT第一年到期 。
现在美国的任何地方都愿意relocate，求建议或ME的工作机会或内推。可以的话请站内
信联系我，我会把简历发过去。
十分感谢！ |
|
发帖数: 1 | 14 求工作实习或博士后机会，本人美国ME专业PHD， 计划今年12月份毕业，
研究方向是有限元模拟(高速撞击，齿轮拟合，柔性结构振动，复合材料固化反应)，断
裂损伤力学，复合材料弹塑性变形，声固，热力耦合分析。有多年使用Matlab和
Fortran写有限元代码的经验，对Abaqus(包括用户子程序:UMAT,VUMAT,UMATHT,USDFLD
等)和Solidwrok比较熟悉。
详细简历可通过站内联系，或者邮箱：[email protected]，联系后我会将简历发到您的
邮箱，谢谢。 |
|
发帖数: 1 | 15 大家好,我是中部一所公立学校的ME博士毕业。背景主要是有限元，力学试验，平时用
到的商业软件
ABAQUS, ANSYS, HYPERMESH, COMSOL, SOLIDWORKS 博士期间主要做的是复合材料本构
建模，力学测试以及有限元分析，编程主要用FORTRAN做ABAQUS用户子程序开发以及
MATLAB代码。
现在已经拿到OPT, 随时可以开始工作。
美国的任何地方都愿意relocate，direct hire和contractor的工作都可以，求ME的工
作机会或内推。
可以的话请站内信联系我，我会把简历发过去。
十分感谢！ |
|
h****a
发帖数: 580 | 16 向各位牛人请教一个问题：
现有一套数据（x[m],y[n],z[m,n])，其中x[m]和y[n]是参数，z[m,n]是（x[m],y[n]）
相对应的值，m从1到M，n从1到N，一共有M*N个数据点.
根据我的观测和初步检验，我认为下面的方程可以很好地拟合(fit)这些数据：
z=f(x,y)=a11*x^2*y^2+a12*x^2*y+a13*x^2+a21*x*y^2+a22*x*y+a23*x+a31*y^2+a32*y
+a33
现在的问题是，如何准确地得到方程中的系数coeffients（a11,a12,a13,a21,a22,a23,
a31,a32,a33)?
据我所知，对于二维数据（x[i],y[i]), 很容易用方程y=f(x)来拟合，有很多工具（
gnuplot，xmgrace）可以做到。但是对于三维数据（x[m],y[n],z[m,n]),请问有什么现
成的工具可以用来拟合？若是没有现成的工具，我需要自己写程序，该如何着手？
我在网上搜索“least-squares fitting”，找到一些说明和子程序，但全都是用y=f(x
)来拟合数据(x,y)，没有用z= |
|
k***i
发帖数: 662 | 17 Matlab的最大优点是调试界面很友好，向量处理的子程序丰富，而且图形功能强大。
缺点是不适于大规模计算。
无论如何，我的观点是，能及时完成任务就行，管他什么工具呢。 |
|
b***k
发帖数: 2673 | 18 也许这个程序用c实现比C++来实现要容易的多。呵呵，我是c出身的。
这个话题似乎转到了如何学习C++的问题。
以前litaihei在板上说过一个方法，我觉得很好。
就是写个例子程序来练习C++的各种简单功能
比如写个matrix class，包括
constructor (default, copy)
destructor
copy assignment
operator overloading(including +,-,*,transpose, inverse, [],etc)
member function (eigenvalue)
output/input overloading
inheritance (bandmatrix, sparse matrix)
template(generize the matrix class)
这个matrix class写完善了后，可以再编写一点小main程序来测试它。
这个过程中也可以理解const，inline，extern等等的keyword
通过这个实践，可以初步理解C++中一些最基本的但是也是最重要的概念。 |
|
S*****H
发帖数: 90 | 19 公元某年某月某日，某著名投行某部门某组某个头要我去跟他谈一下。大概觉得我大老
远去的，只有一个人见我不太好意思，该头谈完了又找了个挨踢部门的日本裔帅哥考我
技术。
该日本帅哥先考我数学，题目是1+q+…+q^n。我用0.1秒弄清了题目，0.1秒把答案写给
帅哥。帅哥感受到了我的鄙夷不屑，说道“还是有人不知道的”。
帅哥再考我编程，让我写一个排序子程序。我想这难不倒我，快速排序烂熟于胸呢。岂
料该日裔帅哥再三申明，任何排序都可以，那我就只好从简了，写了一个三行的冒泡排
序交给帅哥。为了检查我的三行程序是否正确，帅哥写出了大约15个数，然后开始一步
一步地验算。我不禁瞪大了眼睛，三行的程序还要这样检验。有10分钟的时间，我看着
窗外曼哈顿灰蒙蒙的天空。完了帅哥说是对的。就这么结束了。 |
|
s***e
发帖数: 911 | 20 SUBROUTINE rkdumb(vstart,nvar,x1,x2,nstep,derivs)
INTEGER nstep,nvar,NMAX,NSTPMX
PARAMETER (NMAX=50,NSTPMX=200)
REAL x1,x2,vstart(nvar),xx(NSTPMX),y(NMAX,NSTPMX)
EXTERNAL derivs
COMMON /path/ xx,y
CU USES rk4
INTEGER i,k
REAL h,x,dv(NMAX),v(NMAX)
do 11 i=1,nvar
v(i)=vstart(i)
y(i,1)=v(i)
11 continue
xx(1)=x1
x=x1
h=(x2-x1)/nstep
do 13 k=1,nstep
call derivs(x,v,dv)
call rk4(v,dv,nva |
|
j**n
发帖数: 20 | 21 go to www.nr.com....
they have electronic form of Numerical Recipes...and codes written in
both C and Fortran is available....
but the pages are all in pdf format... so you will need to either buy the
source code from the publisher or use some tools to extract the codes from
the pdf files. |
|
a**u
发帖数: 99 | 22 能不能这样算一个多边形的重心?
选取一个中心为原点,如,取各点x坐标的平均,y坐标的平均.
把每个端点和该点相连,就形成乐N个三角形.
然后,调用N次计算三角形重心和重量的子程序.
得到(M_1,C_1),(M_2,C_2)，……，(M_N,C_N)
M_i为重量,C_i=(x_i,y_i)为重心位置.
这样,该多边形的重心为
X=∑（M_i*x_i）/∑M_i
Y=∑（M_i*y_i）/∑M_i
同理，可以得到凹多边形的重心。
整个怪异多边形的重心则易得到。 |
|
r****y
发帖数: 1437 | 23
SUBROUTINE four1(data,nn,isign)
INTEGER isign,nn
REAL data(2*nn)
INTEGER i,istep,j,m,mmax,n
REAL tempi,tempr
DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
n=2*nn
j=1
do 11 i=1,n,2
if(j.gt.i)then
tempr=data(j)
tempi=data(j+1)
data(j)=data(i)
data(j+1)=data(i+1)
data(i)=tempr
data(i+1)=tempi
endif
m=n/2
1 if ((m.ge.2).and.(j.gt.m)) then
j=j-m
|
|
r****y
发帖数: 1437 | 24
此程序four1.for用于追赶法加FFT解椭圆方程会出问题, 需要使用
realft.for如下
SUBROUTINE realft(data,n,isign)
INTEGER isign,n
REAL data(n)
CU USES four1
INTEGER i,i1,i2,i3,i4,n2p3
REAL c1,c2,h1i,h1r,h2i,h2r,wis,wrs
DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
theta=3.141592653589793d0/dble(n/2)
c1=0.5
if (isign.eq.1) then
c2=-0.5
call four1(data,n/2,+1)
else
c2=0.5
theta=-theta
endif
wpr=-2.0d0*sin(0.5d0*theta) |
|
s***e
发帖数: 911 | 25
我用fortran子程序, 你拿去自己转成c的? |
|
m**a
发帖数: 16 | 26 法1: 进行均匀抽样，每隔若干次取平均值，根据中心极限定理，最
终分布趋于GAUSS分布。
法2: 找NUMERICAL RECIPE，有标准子程序。 |
|
|
b*******x
发帖数: 100 | 28 。。什么类型的输出？你是只想要return的value吗？
如果是warning和error信息，可以用 try(exprs,T) 来抑制
如果不是，是程序内的打印输出，他们一般又自己的参数比如trace=F,verbose=F,来抑
制这种输出
如果没有相关参数，试一试sink（），把信息结果送进一个临时文件。
偶知道的就这些了。。。 |
|
b*******x
发帖数: 100 | 29 假设你没有参数能控制，这里有个函数模仿一下，打印一个结果，但是你只想要return
的那个value
> myfun <- function(x){print(x); return(x+1)}
直接运行，又output 1的输出，你不想要他
> z <- myfun(1)
[1] 1
就创造个临时文件，把输出都扔进去
> tf <- tempfile()
> sink(tf)
> z <- myfun(1)＃没有显示那个1
> sink()
> unlink(tf)
> z #但是你有那个返回值
[1] 2 |
|
z**k
发帖数: 378 | 30 f <- function() print("hello cat!")
f()
sink("stdout")
f()
sink()
f() |
|
|
g********r
发帖数: 8017 | 32 TNND. 这个程序不知道是什么底层输出. sink()光把标准输出给除掉了.那些乱七八糟
的灰字一个也没少:( 别人的C程序搞起来太难了.忍着了. |
|
D******n
发帖数: 2836 | 33 try this argument
sink(.........,type = c("output", "message")) |
|
|
D******n
发帖数: 2836 | 35 其实你就让它输出乱码呗。
到时候parse一下输出文件就可以。 |
|
g********r
发帖数: 8017 | 36 嗯也还好.不是乱码,是很多不需要的信息.现在我把需要的sink()出来,连parse也不用:)
就是刷屏厉害,影响速度.mac里面输出是实时的,几毫秒的运算还要等小半秒的显示. |
|
D******n
发帖数: 2836 | 37 。。。。我的意思是用batch mode，眼不见为净
R CMD BATCH xxx.R
用:) |
|
|
s*r
发帖数: 2757 | 39 it has output buffer. it may hold the output for a while and output to
screen in batches
用:) |
|
T****n
发帖数: 2195 | 40 MATLAB数据都是以钜阵形式（或列向量）存储的，SAS得写数组。
MATLAB运行大程序，方便开中间某一个变量的值（主窗口输入变量名）。SAS还得PROC
PRINT。
MATLAB的子程序块比较容易写，也容易修改。SAS MACRO相对来说语句难一些吧。
个人观点 。。。。。以前学过EE，做过点儿MATLAB的东西。
但统计函数，MODEL，TEST什么的，definitely 用 SAS。 |
|
n****u
发帖数: 110 | 41 主程序里面调用了其他几个子程序，
怎样记录各个子程序运行所花的时间？
谢谢 |
|
s********n
发帖数: 26222 | 42 允许包子群发，我指的是不通过wbcenter, id个人可以通过金融中心自助解决 |
|
m****s
发帖数: 18160 | 43 抱歉，这个不会开放，个人转账只能一次转给一个用户。 |
|
h****n
发帖数: 29 | 44 刚才我们说到了复制，是的，这正是我打算讨论的重点：意识的可复制性。
原帖分析了意识的本质，如果同意其中的逻辑，我们可以顺理成章地推导出意识能够复制——
人的意识包括记忆与思维方法。思维方法是可以改变的，不是唯一的，对区分个体来说并不足够。记忆包括：人生成长经历记忆与学习知识的记忆，其中知识的记忆也不是唯一的，只有人生成长经历记忆才是个人唯一的标志，无法更改，是他人没有的。例如：你成长中所经历的所有事情，你对亲人、同学、同事、熟人、朋友的相互关系的记忆，是唯一的。人不断有新的记忆产生，也会忘记许多过去久远的记忆，因此，我们只要有部分重要的记忆，就是“我”了。不同的记忆会促成神经元进化出不同的网络结构。保存住这个网络结构，就等于保存了记忆。
如何保存网络结构？实验表明，不同网络结构的神经元，其信息处理能力不同。信息处理能力指：输出信息/输入信息，或表示为：C/R （C=输出，R=输入）。神经元的C/R就是该神经元的信息值。知道了大脑神经网络中所有神经元的C/R值，并且知道了神经元的连接方式，我们就获得了一个人的所有意识信息，就可以用电子程序仿真一个等效的电子... 阅读全帖 |
|
m***i
发帖数: 2013 | 45 Good point. But theoretically 可复制 is never equal to that you can really
复制it, hehe
能够复制——
体来说并不足够。记忆包括：人生成长经历记忆与学习知识的记忆，其中知识的记忆也
不是唯一的，只有人生成长经历记忆才是个人唯一的标志，无法更改，是他人没有的。
例如：你成长中所经历的所有事情，你对亲人、同学、同事、熟人、朋友的相互关系的
记忆，是唯一的。人不断有新的记忆产生，也会忘记许多过去久远的记忆，因此，我们
只要有部分重要的记忆，就是“我”了。不同的记忆会促成神经元进化出不同的网络结
构。保存住这个网络结构，就等于保存了记忆。
信息处理能力指：输出信息/输入信息，或表示为：C/R （C=输出，R=输入）。神经元
的C/R就是该神经元的信息值。知道了大脑神经网络中所有神经元的C/R值，并且知道了
神经元的连接方式，我们就获得了一个人的所有意识信息，就可以用电子程序仿真一个
等效的电子大脑出来。这就是意识转移进电脑的理论基础。
元，其中每一个又与数千个其他神经细胞相连，如果这么计算的话，大脑里至少有150
万... 阅读全帖 |
|
e*r
发帖数: 1121 | 46 【 以下文字转载自 Programming 讨论区 】
发信人: ear (猪头洗心革面), 信区: Programming
标 题: numerical recipe里的快速傅立叶变换
发信站: BBS 未名空间站 (Tue Jul 1 22:39:16 2014, 美东)
最近需要编个cross correlation 的程序， 调用了numerical recipe里的four1 （快
速傅立叶变换）。用简单的正弦函数做测试，通过。上真实数据，发现傅立叶变换后返
回很多NAN的value。数据点也不是很变态，大概是0.01 到3之间。网上google有些人说
这个four1毛病很多。我不是很明白。
这里有懂的可以讲讲么？另外可以推荐个比较robust 的快速傅立叶变换子程序吗？非
常感谢！我用的是fortran.
xiexie! |
|
c*********r
发帖数: 19468 | 47 纸面可以参考啊，马力、车重、变速器、驱动方式、空动、轮胎、车架、悬挂、各种电
子程序……各种因素都有关系，不见得每次只凭纸面就能准确判断哦 |
|
d******0
发帖数: 1598 | 48 9，感恩节，Naperville，晚8点
没人唱了，都坐沙发上。卡拉OK嘛，我知道是怎么回事：晃常吼的那几首过上一遍，心
里憋那股劲儿发出来，也就歇菜了，跟感冒出一场大汗差不多。
半张墙那么大的电视依旧字幕闪烁。你的脸有几分憔悴，你的眼有残留的泪。背景乐海
潮般起伏，人声部分被单调的电子音代替。李蔷薇和她那杯可乐加冰，你的唇美丽中有
疲惫。这姑娘手里还握着麦克，两只手一上一下紧紧握着，黑黝的柄，灰白的头，幸好
她没有张开嘴——你的唇美丽中有疲惫——对着它唱。我暂停住“垫子”里的“国内十
大禁片”，对她笑了一下。她也对我笑了一下。
我和李蔷薇交往半年了。同年出的国，不同的专业，小小一个大学城，在散发着鱼腥味
的中国店对面相见不相识。直到临毕业才在某某某的生日聚会上通了姓名，留下电话，
去酒吧吃汉堡包，独立日的夜的草坪上铺张毯子对着烟花嚼爆米花，听了半场学生票半
价的交响乐，送她回公寓楼，隔了座拥抱，倚着窗接吻，手牵手上楼，晨雾弥漫的
Wendy’s drive in，两杯咖啡一份墨西哥鸡肉卷，最后直接送她去实验室。哪怕是一
套程序，也可以浪漫起来，计算机科学博士生李蔷薇如是说。这姑娘看起... 阅读全帖 |
|
