极差中有负数运算法则要怎么运算啊?比如17,1...

一个曾经体质极差的人怒答。这不是如何健身的回答,只是一个快要把自己玩死的边缘,把自己拉回来的综合路程。一个吸烟9年,喝酒,作息不规律,熬夜,爱吃重口味油盐等等..的人,一次改变的路程。曾我真的以为我还能呼吸就真的能多活一天似的。大学毕业后的数年里,持续不断的,一个不拉的,坚持不懈的干了如下事:(真特么想抽死自己)
相信很多毕业之后身材走样的人没少干以下的事,其实都是自己的错,怪不得别人。抽烟 :一天至少一包半,前后抽了9年。以至于,后来只要稍微跑动,肺部整个一片灼烧感,有时候有即使不跑动也是灼烧甚至有血味(感觉这尼玛要死了差不多了)..喝酒:有一定程度借酒消愁的行为...作息不规律:4,5年没好好规律生活过了。熬夜:这不是家常便饭麽?顺便抽着烟喝着酒,吃着地沟油烧烤。打两盘DOTA...不极度疲劳不睡..持之以恒坚持不懈的干这事。重油重盐辛辣油腻,大口吃肥肉:一天没一斤肉,加上肥肉..怎么可能。而且全是油腻腻的。而且都是在晚上吃.. 吃很多。坚持不懈持之以恒的吃。 吃完了还加餐出去吃的。还来薯条,喝各种饮料...还特么不运动:我觉得很长一段时间我最大的运动就是下楼走一圈。过度的XXOO (捂脸跑走...):那会年轻没见过女人,跟种马似的....就算体质很差了,还是乐此不疲,战个痛快。不知不觉..人还是会被抽干的。大概就是上个楼都喘,各种失眠盗汗,一天到晚犯困,萎靡。要么昏睡不醒,然后醒来2小时就犯困。虚成这样子!以上都要我一个不拉占全了.. 坚持了好多年...(羞愧难当.我这么糟蹋自己)你们谁能同时把以上全做到了,而且坚持不懈很多年还没挂了...我现在回想起来,那时候的体质虚弱的和快死了没啥区别。却浑然不自知!!后来3件事刺激到我了。第一件:EX嫌弃我20多岁人,30多岁发福的身体。感觉奇耻大辱。我也不知道为毛觉得奇耻大辱。体脂图,哥们我的身材当时就冲着 这个20-22%的造型去了,还是没肌肉的那种。再这么下去就是30%。(我想我没病,很肥的时候不会玩自拍的...我拍的下去麽我...不然真是可以给你们看看)咳咳再强调一次 ,这不是如何健答..只是一个综合过程.. 作为给题主的参考。第二件事:被EX甩了之后,我的身体好虚弱,几天没合眼,睡不着啊。整个精神神经处于应激状态,关在一个小屋子里闷着抽烟,一包接一包。(想象一下那种乌烟瘴气里被熏几天的感觉吧)然后摸鼠标键盘的时候,肌肉都在抖...后来才垂死的睡着,醒来的时候汗把床被子全打湿了。我当时要是挂了。我一点都不意外。此时还发生了一件事,一个快递大哥来送快递,她老婆正好也在,她用一种惊诧的眼睛瞪圆了看着我。 (和快递大哥认识住一个院子,自此之后比较熟了,这大约是昏死睡之后一周的事,还不是当天..)之后我才知道她在惊讶,怎么一年轻小伙子,憔悴成这个样子...第三件事:等我缓过来之后,终于垂死的睡过去。一身湿透的盗汗醒来,有了一点力气,挣扎着爬起来,面对这一切,我觉得我需要锻炼了。 (这已经是那次盗汗昏睡 2周以后了.. 觉得有点力气了才去..)在跑步机上,跑了2公里,很慢的速度,和小姑娘差不多,将近20分钟。下跑步机的那一刻:天旋地转!差点没晕过去,肺里全是血的味道。眼前发黑。(现在想起来都后怕)我当时要是挂了。我还是一点都不意外。在这三件事的综合刺激下,我开始决定改变一些什么:只是为了活下去,对好好的活下去。题目主可以做个参考。做做修改,以适应你自己的方法。我详细的记录过程是想告诉LZ,一切的改变是环环相扣的。不是一步到位的。我也是无意之间环环相扣影响着做了一个综合改变。希望给你参考借鉴。我做的第一件事就是跑步:从2公里发晕到跑完5公里无压力。这个过程我用了2个月。一开始多么的缓慢。真的很累,我每次去健身房的时间,大约是7点,要工作学习啊。然后锻炼要各种等,跑步机满。然后9点才能练完洗完。 还没吃的,要走半小时到地铁站再坐一小时地铁回家。回家后再做吃的睡觉。可能那会神经真的受刺激了,居然要我坚持下来了,差不多一周3-4次。如果你觉得跑步很无聊,我觉得你需要一个很好的耳机,我就靠着这个撑过很多无聊的时候。别再扯着耳机线了。晃啊晃啊你不难受麽!真的很蛋疼好麽。拥有这个:你就会觉得自己飞起来: 跑步很欢乐了!音乐自选,就激情一点吧:这两款运动耳机我都有,而且我都用了很久。... 综合下来推荐索尼的。这时候跑步终于变得有趣多了!不用烦耳机线了!重要的是,你是听歌还是听英语还是听电影随便你,然后一切就在你的 ‘似魔鬼的步伐’中 ,一步两步...摩擦摩擦....咳咳(神曲移步:)这个时候气色开始恢复,以前都是面色发黑的,后来熟人说,最近你气色好多了,脸红润了,终于有色泽了,之前你那憔悴的样子。(看来每个人都看得出,就我自己不知道- -)第二件事就是戒烟戒酒:我觉得跑步这件事简直是一切的基础:希腊:如果你想强壮,跑步吧!如果你想健美,跑步吧!如果你想聪明,跑步吧!无意之间,锻炼的习惯与跑步要我打下了戒烟戒酒的基础。虽然坚持锻炼,但是烟酒不断也不是什么好体质。终于在第三个月的时候,我把烟戒了。 在猛抽了3包烟之后,掐灭了我人生最后的一根烟。然后无痛苦的戒掉了,一点念想头没。 9年的烟民,仅仅4-5天我就把体瘾收拾干净了。如果要给戒烟的人什么建议。我可以另开一贴了。但是第一步一定是锻炼,第二步破心瘾,好好想清楚吸烟你得到了什么。我巧合的完成了第一步,在其他书籍的帮助下猛抽烟看完想清楚。书籍: 这本书能让你戒烟豆瓣评价:(真想个给作者32个赞)如果LZ有什么不良嗜好抽烟喝酒。先从跑步开始,再怒攻心戒之。 跑步会为戒烟打下很好基础。我一个抽的肺部成天灼热,一嘴的血丝味道的人,都能开心 ‘无痛苦’的戒了。对戒的每一秒我都是快乐的,舒服的。所有的戒断反应都是挣脱枷锁的赏赐。之后戒烟至今,一点念想都没!最重要的是:我很开心!戒烟很爽!一点也不痛苦!包括过程!相信我,一根接一根的烟接着抽,简直没有比这个更无聊和愚蠢的事了。完成了戒烟,戒酒就更简单了。这个我本来就没瘾的。第三件事改变作息:当我完成跑步与戒烟这种习惯之后,心肺有了一定恢复,起码我能睡着不大肆出虚汗了。但是这还远远不够。我还有个大问题!作息不规律,到点睡不着.... 这也是健康的大硬伤吧!如果这个不改,其他一切都是白搭。我无意之间看到一句话:你之所以即使疲劳了,也不肯睡觉,是因为你没有勇气结束这一天,也没有勇气迎接新的一天。参考 第二高赞。非常感谢答主。真的帮到了我。一下我脑洞打开,我颇花费了一些日子,告诉自己这一天无论自己做的如何,有没有拖拖拉拉。有没有完成这件事。我都要结束它了,如果想要自己睡的安心,那么白天就积极点。但是!还是有问题的,虽然我努力转变了观点,但是到点还是多少有点睡不着。不过这又开启了第四件事。第四件事打坐,冥想:我觉得是白天的自制力不够,不能好好完成事情,才让自己晚上安心睡着。于是我转而研究了一下自制力,确定它是一种‘肌肉’ 是一种能量,你需要刻意锻炼它,也需要它休息。一旦它疲劳了。你的行为就崩溃了。 所以,我们会发现当我们疲劳的时候,会更加容易暴饮暴食,更加容易的发怒等。参考书籍: (我一开始以为是励志鸡汤文,一直嗤之以鼻,后来无意看到,才发现这本书是斯坦福一门课程。有理有据,我其中借鉴了不少有效的经验...干嘛取个如此鸡汤的名字,叫什么 “前额皮层操作术” 多牛掰.)随着我的实践一路走来。其实对于以下观点不能同意更多:很多问题,不是‘意志力’解决的,其实就是个体力问题。我们中国人,太注重想通了,好像想通了就什么都解决了。其实不然,道理那么多,我们还是过不好一生。别说一生,连个觉都睡不好其实。 很多都是身体问题其实,我觉得这点西方就很务实。锻炼它的第一步就是冥想,打坐。专注于呼吸之间。一开始专注于呼吸之间,坐3分钟我都坐不住。但是我知道自制力是需要锻炼的。专注于呼吸之间的冥想就是最好的打基础,不但可以提高注意力也能提高自制力。于是我每天睡前都进行练习,带着耳机。听很清淡的曲子,但是尽力不跟着调子跑了。慢慢放松,专注呼吸之间。过去谁能想象,我这个坐着就要狂抖腿的人,能像现在一样安安静静的坐着打坐。随着时间的推移与累积。你会在中间学会很多。很重要的一个附带就是,你会意识到你的很多情绪,观点,一股一股往脑中冒的想法都是那么的多,不过你将学会的是,不用刻意压抑什么,任由他们生,也任由他们灭。你只要专注于你的呼吸之间,哪怕你意识各种想法无论是惊涛骇浪 还是 香艳无比 你只要默默的专注于你的呼吸之间,看着他们起看着他们灭。由于这件事,我无意之间加上第三件事的观点:只是你没有勇气结束这一天。有了一个很好的结合。 冥想之后,我觉得这一天我给它画了一个句号。于是我能安安静静睡着了。即使偶尔睡不着,我也会告诉自己,这一天结束了,即使我没睡着,我也休息着。特别坦然。结果是..必须睡的着啊!于是我的作息能够保证了,这还是不够把身体变好,只是叫做正常,叫做恢复。于是水到渠成的第五步有了第五件事,健身:一段日子后,我觉得这样算是稳定了,起码和正常人差不多生活了。再也没那么虚弱的要死了。也没有肺里的血气,但是里强壮和强大差的远呢。于是有了个新的想法;我要强壮。当一个男人发自内心的有了我要强壮的念头之后,一切都好办了。这里请参考:知乎,硬派健身系列。 我都跟哪学的。不过我个人从此不往健身房跑了, 没那时间,也不喜欢对着一堆金属吭哧吭哧的。我个人比较喜欢无器械,囚徒健身等。因为相对自由,个人喜好。 殊途同归。不要争好坏。都挺好。书籍:用失传的技艺练就强大的生存实力优酷: 动作参考。很多训练的动作和基础:可以在知乎上的硬派健身多做了解: (那些要减肥,健身的看这里好吗!你们想要的一切都在这.我可不专业噢, 非常感谢
你的文我全部拜读.受益匪浅)再此我不得不再次强调,这不是一个健身贴,这是一个从颓废走向健康的心路贴。专业知识,这里水平不够没法提供- -,也不做任何误人子弟的指点。OK继续,所以后来慢慢的从跑步,转向了无氧与HIIT锻炼。insanity08
地址:从一开始10分钟跪地上,到现在无氧完了还能跳30分钟。 如果只是做做核心肌肉的话,我能跳完全成。 这是一个对我而言很好的进步。尤其是你早上起来之后,发现,额,我胳膊上怎么这里突起一块肌肉,额过几天又发现, 额原来这个部分是有肌肉的,原来这一块是几块肌肉组成的... 你看着自己身材的变化觉得惊喜不已。是我目前最大的快乐之一。相信也会给你带来快乐。我主要是在家里练的,这一路要我意志力有了巨大的磨练。 无论刮风下雨,多晚回家。我肯定锻炼妥妥的。最晚一次我记得练到夜里1点,虽然不科学,但是只是偶尔,有时候坚持的意义比偶尔的不科学重要。随着体质渐渐的增强。慢慢就开始注重减肥与长肉了。一个曾经毫无自制力的,虚弱的走上发福道路的人都能这么练。50个仰卧起坐就嚷嚷腹肌为何不来...还是退下吧。为了继续加强锻炼,我开始用一个旧书包装了一点铁,注意,这并科学,不要随便模仿,但是一定要做到自己瞧得起自己,尽自己最大努力。然后背上它:完成200个俯卧撑。然后脱掉....然后100多个卷腹,加上10-30分钟HIIT。 这是某一天内容。每天我的的运动内容不同。无氧+HIIT。(9.16 更新:时间我已经修改方案与强度缩短至1小时出头。)这特么才叫减肥...别50个仰卧起坐就嚷嚷腹肌怎么还不来。你们觉得已经走到今天的我,还是个虚弱的胖子吗... 咳咳,不过我是暂时比不了腹肌男们。第一我条件不够。第二我基础薄弱。TAT臣妾也不是神啊...还不得脚踏实地一点点去坚持去做吗 TAT(9.16更:自重与无器械的弊端:本人已经在重新考虑过段日子搬家,到健身房附近)加油吧,痛苦只在一瞬间。它顶多一分钟,一小时,甚至顶多一年。我用这个和志同道合者打气:PS:送上一些我锻炼时候,累的不行的时候用来打气的视频:第六件事,饮食: 当前5件事养成习惯后。这件事就水到渠成的来了。作为一个重油重盐,吃辣,无肉不欢的人。我也不认为我的饮食也可以做一个转变。我也干过半夜突然受不了馋死了跑出去大吃一通。但是如今,蛋奶蔬菜粗粮为主。早餐:香蕉脱脂牛奶燕麦葡萄干 是铁定的晚餐同早餐。午餐:绿豆粥黑米粥各种粗粮蔬菜。 肉:.当然还是吃的,但是无油无盐(不放油盐,白水煮,或者放点萝卜土豆吧,其实挺香单纯的肉味啊,土味啊) 定量,在健身后,吃100-200g (还不是每次,大约2-3天才吃一次);想起以前晚上一顿吃2斤油腻腻的肉就恐怖。喝:杜绝了饮料, 只喝牛奶和水,绝大部分时候。 顶多偶尔锻炼完会喝一点功能性饮料。水是过滤过后的水,和自来水直接烧开了喝伟大很有差别。 要尽量注重水的品质。现在有很多过滤饮水机,LZ买的这个:感觉比自来水烧开强多了。EX可以见证我家冰箱简直是各大饮料厂展览。只喝饮料不喝水的。我曾经一度觉得,不喝饮料怎么活的下去。悲剧的EX,我又把你拎粗来了。提起她,在无亲无故的这座城市里,独自经历了这一场之后。如今,我只能坦诚的说,爱过。看客们请淡定,咳咳。我找到了真爱:自爱.
TAT饮食调整的经验就是:1.如果你突然特别渴望高热量的食物,比如KFC啊 大肉 油腻腻的食物。 那就是很可能你体内需求热量不足了。你可千万别委屈自己,忍着不吃啊,那会让你吃的更多。该吃就吃!但是更好的做法是:摄入蛋白质与有碳水化物的食物:我一般就是 香蕉 葡萄干 牛奶 蛋白 等。 只要一小会之后,你就不会对这些油腻腻的高热食物非常渴望了。甚至吃不吃无所谓了。 因为你成功的把大脑给‘忽悠’过去了。它获得满足了之后就不会折磨你了。你只要给他足够的蛋白质和碳水。(果然讲科学 讲人性才是硬道理)高蛋白,低脂肪,高碳水化物,高纤维素 是你更好的选择:食物热量查询:参考:2.你会渴望吃肉的,但是你知道知道,你需要的是肉,而不是加满佐料各种奇怪的食物。(这不科学,有调料才是美味,我只是怕自己重口味复发,刻意的尽量淡。 完全不用学我,不过少油少盐对身体很好)我大量锻炼后的土豆牛肉:就是土豆+牛肉+一点料酒去腥味。 生炖。 无油无盐。你这时候才会知道,你其实只是需要肉,不是一堆油腻腻的配料。(大概就是,你饿了3天,干馒头你都觉得香,所以..对我而言这时候没什么油盐我也吃的很香)一个重口味,食肉动物主食现在的样子:脱脂奶粉+桂格燕麦+葡萄干: (其他若干香蕉, 黑米,绿豆,黄豆等...就不拍了..总之几乎就素为主。素中我又偏好以 黄豆绿豆黑米土豆萝卜白菜空心菜西兰花这些...尤其是豆类 蛋白质很高,碳水化物丰富,热量并不高,纤维素丰富,锻炼之后会加餐肉食为增肌,个人喜好,非专业推荐。仅供参考试试)若是加餐嘛:比如牛肉炖土豆: (就直接炖..加了点料酒去腥味,我自己反正是觉得很满足了)(补充:很多人误以为我完全油盐不进..怎么可能,臣妾做不到啊...我还活不活了。土豆和香蕉有大量的K+, 其次作为一个人类,肯定下意识的缺什么会去补充什么的。很多食物里都有一定的盐和脂肪含量的。)如今,你要问我痛苦吗?一点都不。一开始有点不习惯,但是慢慢你会觉得浑身清爽。你要永远记得,你要改变体质,也同时意味着要与过去的生活方式做一个诀别。以上。如今的我每天吃的好,睡的好。精神很好,肺里每天都很舒畅,每天醒来浑身都充满着能量。 我在呼吸,我的心脏在有力的跳动。健康很重要,我很幸福。10.7 更新:(内些个人气旺盛的微博,在转载文章的时候,请注明一下来自知乎什么的...谢谢咯;身高180 BMI:21.9 最近正在努力增重,长肌肉中;微博:
;)PS:调侃一句,人家星座都是秀的恩爱,摩羯都在那秀坚强。 TAT
今日更新牙齿部分,正文在最后。比较短小精悍,但十分重要。-----------------------------------------------------分割线-------------------------------------------------------------------------关于体质极差这点,幼年的我真是相当符合。跟现在排名第一的@梦尧不同,身体不好是天生的,而且我是个女生,所以也算作一个小小的补充供大家参考。先细数一下我从小到大的疾病史:1.
过敏体质,对几乎所有抗生素过敏,发烧了就只能喝热水。2.
长期生病。小学基本有一半时间都不在上课,要么在家里,要么在医院。3.
哮喘,经常发烧,偶尔还带血。三年级的时候一边喘一边40几度烧了整整一个星期,又因为对抗生素过敏,光靠葡萄糖抗过一周。4.
不能剧烈运动,走动超过15分钟就胸闷难过,在大太阳底下晒5分钟就晕。再说下现在的情况:1.
慢跑10000米基本没什么问题。2.
平板支撑最长7分半,杠铃深蹲35公斤。(其他数据就不列了)3.
跆拳道黑带。4.
有腹肌,腰臀比一直维持在0.68~0.7之间。5.
基本不生病了,偶尔小感冒喝点热水注意休息3天左右就好。6.
每天睡眠6~7小时也很精神。说一下开始锻炼的契机吧,大概是在初预的时候。因为孤独。第一,
总是卧病在床,从小陪伴在我身边的就只有家人和书,不是说不好,而是说不够好。我很想交朋友,也很想跑跑跳跳,但是由于羞涩和体力差,总是和同学玩不到一起。因为体质差,我的童年就有些单调和无趣了。第二,
喜欢上一个男生(开窍开的比较早,哈哈)。他体育很棒,平常爱打羽毛球。我深知自己配不上他,但我总想如果能站在他身边,就算是做朋友也很好啊。没有结果的暗恋,却让我变得更好,很感激他。第三,
舒婷的《致橡树》。幼年也爱读诗,尤爱这首。那时虽然年幼,却也期待两人相互独立,尊重钦慕的爱情。如果身体糟糕,又如何能做到人格独立,不依附于他人呢?现在想想也挺好玩的,长大之后成了个女汉子,小时候却是这么多思柔情,真是退步了。下面是从开始到现在的锻炼历程。体质差的人锻炼起来最怕的就是操之过急,大学城里晨跑跑死人的事也层出不穷。所以我建议刚开始锻炼时先从最轻量的开始。以下的话只是我个人的锻炼过程,仅供参考。一开始是走路,这个或许大家可以直接略过,因为这太平常了。每天放学之后不做公交车,慢悠悠地走回家,回家的路的话正常速度半个小时左右。第一天开始尝试花了90分钟,把我妈急的。后来老妈知道我要开始锻炼身体了,每天就骑着自行车,把我的书包放在前面的筐子里,陪着我慢慢走回家。也在这里谢谢我妈,毕竟刚开始的时候最痛苦,有人在身边陪伴我鼓励我也给了我很大的力量。这里还有一点小小的建议,大家平常坐地铁坐公交的时候可以不坐电梯,多走走也是好的。从慢走15分钟就晕到正常速度走30分钟无大碍大概花了3个月吧。之后就开始打羽毛球,心怀不轨地跟着爱慕的男生学,趁机凑到他身边去,哈哈哈。这个时候体质也开始改善了,至少不用一个月看几次医生了。羽毛球对于女生来说真的是很不错的运动,可以将身体练得修长,日常打不用耗费太多体力,也好上手。这个时候觉得身体开始变得轻松,我开始意识到之前没有生病的时候自己也是不健康的,并且开始单纯地渴望健康的。之后开始我这一生最重要的爱好——跆拳道,它给了我健康的体质,自信心,礼节与克制。今年20,跆拳道练了有10年了。当然,一开始就给我一个下马威。我本来觉得身体已经好的跟一般人差不多了,结果还是跪。体能比不上同龄人,韧带比不过小朋友,平衡感差,手脚也不协调,综合能力倒数第一(手动撇嘴)。我教练(又是一大恩人)也没嫌弃我,就算我跟不上,他也会默默地安排我做些力气花的少点的活动,但他从未叫我停下,因为他知道我想练,就算很苦也有人嘲笑。每天训练结束之后,我都会留下来跟教练请教一些问题。如果体质差,就更不能瞎练。我就是个活生生的例子,因为刚开始锻炼的时候没有注意膝盖,现在右膝盖永久性损伤,泡冷水或者运动过量都会很疼。这个时候有人在旁边帮你一把,认真地告诉你怎么做最好,这很重要,毕竟在我小时候还没有知乎么。不过现在你可以看@高科教练和@李浩轩教练的回答,知乎上的一些健身建议都被我用在日常锻炼中,这里就不赘述了。每周6个小时的训练,而寒暑假,基本上就是泡在道馆了。这也说明不把锻炼当成负担,真正爱上规律的锻炼对体质提升大有益处啊。到了高中,也就是跆拳道练了4年之后,我终于可以自信地大喊一声:老子的体质跟正常人差不多了!(体能早就超过了,这里指生病频率低了)本来小时候的初衷是变成普通人的体质,少生病就够了。但我开始不满足,因为我想试试看,如果体质再上一层次会不会让身体更加舒适,更少生病,身材也更好。于是就去打网球了。结果真的是血泪史啊。我勒个去,姑娘千千万万不要打网球啊。虽然身体素质提升地超快,但是这右胳膊比左胳膊粗一圈我实在是不能接受啊。男生倒是可以尝试一下,贵族的运动嘛,单手反拍超帅。(好像偏题了不好意思)在练跆拳道打网球的同时,我也一直在打篮球。虽然身高比较矮(只有160),但是有幸交到校队的朋友,当了校队的陪练。因为打篮球需要更多的跑动和团队能力,这也算是提升了另一层次的体质了吧,也交了不少好朋友。这个时候半场比赛打下来会累,但是睡一觉就好了。先写到这里,放个预告。上面是提升体质篇,之后也会补充更详细的内容,明天写上大学之后的健身篇(进阶)。放两张图在这里,求赞╭(╯3╰)╮(不要害怕啊,女生一般不会练成这样的腹肌的,我稍微吸了下肚子= =)宿舍乱,不要在意这些细节。。。。。。觉得对于题主来说,或许方法并不是那么重要,开始练起来才是重要的。每个人都是独一无二的,看看别人的经验,在思考的基础上不断修正自己的锻炼方法比较好。写下这些文字,或许没有多少人会从头至尾地看完。但我依旧期待这个题下的回答能给那些还没开始的人一些鼓励,一个契机。这是我的初衷。8.2 update关于牙齿我的牙齿还是比较整齐的,但是和一般人不同的是,由于体质弱,貌似到2岁才开始长牙(一般是4~10个月),然后到18岁才彻底换完牙= =。小的时候不注意护齿,不喜欢刷牙,所以有两颗牙齿蛀掉,高中的时候因为一次事故,右门牙整颗掉落,现在装了颗假牙。不过也因为这次事故,让我打开了新世界的大门,医生认为我的口腔环境有点糟糕,给我提出了一些建议,我也好好遵守了。现在每天起床后刷牙+漱口水,晚上睡觉之前牙线+刷牙+漱口水。一年一次全面检查,半年洗一次牙去除结石。周围的人可能会不解或者觉得我有些强迫症,还有些朋友时常提起牙线会使牙缝变大或者洗牙对牙齿不好这类不科学言论,基本我就一笑而过。因为当时掉了一颗牙就给我的生活带来极其糟糕的影响,所以我就想把牙齿掉光的年龄再往后推迟再推迟。反正那些反对我的人牙齿都没我白。哼。(0)(2)(1)(15)(2)(4)(1)(1)(2)
模型不能算或不收敛,都需要去monitor,msg文件查看原因,如何分析这些信息呢?这个需要具体问题具体分析,但是也存在一些共性。这里只是尝试做一个一般性的大概的总结。如果你看见此贴就认为你的warning以为迎刃而解了,那恐怕令你失望了。不收敛的问题千奇万状,往往需要头疼医脚。接触、单元类型、边界条件、网格质量以及它们的组合能产生许多千奇百怪的警告信息。企图凭一个警告信息就知道问题所在,那就只有神仙有这个本事了。一个warning出现十次能有一回参考这个汇总而得到解决了,我们就颇为欣慰了。
&& & Fixed time
is too large
&&&Too many
attamps have been made
&&&THE SOLUTION
APPEARS TO BE DIVERGING.
&&&CONVERGENCE
ISJUDGED UNLIKELY.
increment required is less than the minimum specified
这样的信息几乎是无用信息(除了告诉你的模型分析失败以外,没有告诉你任何有用的东西)。宜再查找别的信息来考察。根据经验,改小增量步也不一定能收敛,虽然也有人报告过改好的先例,我是从来没有遇到过,也从来没有那个奢望。所以我一般从模型的设置入手。
必须说明的是:Error和warning的性质是完全不同的。Error意味着运算失败,but出现warning可能还能算,而且有些运算必定会出现warning(比如接触分析必定出“负特征值”,下有详述)。很多警告只是通知性质的,或者只是说明一下而已,不一定都是模型有问题。比如以下warning完全可以忽略:
&&&xxxxx will
(not) printed,这种只是通知你一声,某些玩意儿不输出了。还有:
& & The parameter frequency
cannot be used with the parameter field. It will be
ignored(都说某某被ignored了).
如果模型能算,且结果合理,那么大部分警告信息可以不管。但是以下除外:
&numerical&&sigularity(数值奇异):刚体位移(欠约束)
problem. numerical&&sigularity
when processing node105 instance
&&&pile D.O.F.
1&&ratio=1.735e13
pivot(零主元):& &&
&过约束或者欠约束。
这2个问题一般都意味着模型约束存在问题。1)、2)都会伴随着产生大量负特征值。解决方案当然第一步是检查约束了。
有一些直接导致计算aborted,那就得仔细分析了,比如:
&&&xxxxx is not
a valid in ABAQUS/Standard(告诉你这种计算standard不支持了,换别的)
&&&missing
&在perperty步检查材料属性是不是都加上了。如果有梁单元,看看梁法向定义对了没有。
&&&Detected lock
file Job-1.lck. Please confirm that no other applications are
attempting to write to the output database associated with
& this job before removing the lock file and
resubmitting.
&删除.lck文件就可以了,它是一个自动生成的文件。你也可以另存为(另取名),再运算。
&&&The rigid
part xx is missing a refernce point
&刚体(or刚体约束)都必须通过stools--reference
point给它定义一个参考点(RP),载荷都加在这个RP上。
&&&The area of
54 elements is zero, small, or negative. Check coordinates or node
numbering, or modify the mesh seed.
& The elements
8&&have been identified in
element set ErrElemAreaSmallNegZero
& 这个一般是节点编号不对的问题。必须是逆时针方向。
&&&The value of
256 MB that has been specified for standard_memory is too small to
run the analysis and must be increased.
& minimum possible value for standard_memory is
ABA的问题:集合和面的几何的名称最好不要用特殊符号和数值(特别是从hm转过来),全部用英文字母是最安全的。
&&&令很多人抓狂的error
&& & 1)使用了子程序,
子程序有问题, (例如数组定义跟实际赋值不一致,这个我也遇到过)
&& & 2 )模型有问题,
通常模型很大,很复杂(这个我没遇到过)
&& & Please make
sure that the mesh density of the slave surface in the tie
&& & pair(
mbly__pickedset37_cns_,assembly__pickedsurf36) is finer than the
master surface.
analysis may run slower, may yield inaccurate results, and may
require more memory if this is not the case
)硬盘没空间了(这个其实不会引起error code 5,但是出错是肯定的了), 或者是内存太小.或者产生的文件太大.
)关闭杀毒软件试试(特别是卡巴)
& 5) 有人认为边界条件不正确,也会引起这个错误.&
&system error code 29539
&& &关闭杀毒软件and
&&&**ERROR:
Issue cannot be deleted Not all data Released
在windows中,单击“控制面板”--“系统”--“高级”--
“性能设置”--“数据执行保护”命令,
&&&把pre.exe和standard.exe添加进去。重起动后尝试
&Surfaces associated with analytical rigid part
MANDREL may have their orientation flipped
刚体相连的接触面方向定义反了,在接触定义的地方edit--flip
& CONTACT PAIR
(ASSEMBLY_BLANKBOT,ASSEMBLY_TIE-1_DIEDURF)&&
BLANK-1.5 IS&&OVERCLOSED
&&&BY 0.0512228
WHICH IS TOO SEVERE
&&&这往往是因为接触面的法线方向定义反了。定义刚体和shell的surface时,
要注意选择外侧
&123456 elements are distorted。Excessive
distortion of element number 5 of instance PART-1-1
&&&如果有子程序,一般不是材料设置有问题,就是边界条件的问题
&& &XML parsing
failure for job 1.&&Shutting down
socket and terminating all further
messages.&&Please check
&& & the .log,
.dat, .sta, or .msg files for information about the status of the
& The number of history output requests in this
ABAQUS analysis (&5000)
&& &may cause
SIGNIFICANT performance problems during analysis and
postprocessing
&输出项太多,恐硬件资源不够。要是你确保硬件够,这条也不怕了。一般的,应该减少History中的输出项,尽量输出你最感兴趣的内容。
& Value for parameter nset will be truncated to 80
characters
&nset名字取太长了,80字符限制
& compilation - ifort.exe 问题
&&&Problem
during compilation - ifort.exe not found in PATH.
&&&安装的时候没有装好或是二次开发版本冲突。检查环境变量的设置;然后
&&&verify一下,看看是子程序功能否能通过?
如上所说,有很多warning并一定意味着你的模型存在问题。常被问起的有:
1&&负特征值问题
& & THE SYSTEM MATRIX HAS 8
NEGATIVE EIGENVALUES.
&负特征值是非线性分析的必然产物。所以不必大惊小怪,甚至久而久之,对于你熟悉的问题,你都会视而不见了。若出了问题,可先检查下有没有伴随的
&&&numerical&&sigularity(数值奇异)和&
&Zero pivot(零主元)产生。如果没有,可以参考这几个方面:
&1).刚体位移
&2).单元异常,过度变形、过度扭曲等
&3).应力应变关系有负斜率
如果有流体的话,在容器发生形变的话,也可能出现negative eigenvalue
的情况,不过不会出现警告,这是被允许的
&&&The ratio of
deformation speed to wave speed exceeds 1.0000
&&&这个警告是指单元形变速度V(单元最大形变率/特征尺寸)和膨胀波速C
&(通过材料本构关系求得)的比例超过1。
&&&解决这个问题的方案有以下几种:
&(1)检查单位是否封闭(参数设置有数量级的错误),此错误新手常犯;
&检查网格质量 ;
&检查加载速度,如果条件允许的话就降低速度,该方法也很有效,但在很多
情况下无法降低速度;
&调整STEP中的TIME SCALING FACTOR;调整STEP中的
SCALING FACTOR;
&(5)& &加*SECTION
CONTROLS,NAME=SC,DISTORTION CONTROL,
&& & LENGTH
RATION=0.1
或者YES也可以,加在MATERIAL 前面;或加* DIAGNOSTICS,
&& & DEFORMATION
SPEED CHECK=OFF;
或者加*DIAGNOSTICS, CUTOFF RATIO=RATIO(具体数值),在其他
&方法修改后还有问题的的情况下使用增加
&关键字的方法见(17楼)
force/ZERO MOMENT问题
&&&THERE IS ZERO
MOMENT EVERYWHERE IN THE MODEL BASED ON THE DEFAULT
&&&CRITERION.
PLEASE CHECK THE VALUE OF THE AVERAGE MOMENT DURING THE
&&&CURRENT
ITERATION TO VERIFY THAT THE MOMENT IS SMALL ENOUGH TO BE
&&&TREATED AS
ZERO. IF NOT, PLEASE USE THE SOLUTION CONTROLS TO RESET
&&&THE CRITERION
FOR ZERO MOMENT.
&这个警告是告诉你模型中没有弯矩,没问题的,可以继续计算。
&如果提示中出现特征值奇异的时候才是计算有可能出现不收敛的问题。
&& & Degree of
freedom 4 is not active in this model and can not be
restrained
&&&有限元软件计算对于实体步考虑转动自由度,所以你在边界条件中限制了456的自由度后,软件会忽略的啊.
&& &The option
*boundary,type=displa check status file
between steps for warnings on any jumps
&& &prescribed
across the steps in displacement values of translational dof. For
rotational dof make sure that there are
&& &no such
jumps.All jumps in displacements across steps are ignored.
&你采用了位移边界条件,但在平动自由度上,可能在不同的分析步骤里面有突变(你可以从sta文件里面查看),
&并且应保证转动自由度无突变。
通知性质的warning,一般是因为你采用位移加载方式,都出这个。
&& & The strain
increment has exceeded fifty times the strain to cause first yield
at 377 points
检查下约束够不够,约束够了就不用管了,这只是通知你,你的模型塑性应变很大,一般没多大问题。
&&&123 nodes are
used more than once as a slave node in *TIE keyword.One of the *TIE
constraints at each of these
&&&nodes have
been identified in node set WarNodeOverconTieSlave
&&&定义接触的时候,公共节点重复定义了好几次,这样可能会出现过约束问题(只是可能影响)..
&&&There are 2
unconnected regions in the model.
&&&可能是接触面由空隙,最好在接触属性中定义一个容差范围。一般各个parts之间定义接触,aba都会这样通知用户的,只要接触设置对了,
一般没事。
&& &Boundary
conditions are specified on inactive dof of 124 nodes.
&&&The nodes
have been identified in node set WarnNodeBCIactiveDof
&&&边界条件定义的有问题:在124个节点的非自由度上有边界加载
plasticity/creep/connector friction algorithm did not converg
&一般是塑性应变太大,单元扭曲导致的。可以先改为弹性模型看看是否收敛;
&&&The ratio of
the maximum incremental adjustment to the average
characteristic
&& & length is
1.82846e-02 at node 10868 instance jiti1 on the surface pair
(assembly_jq22,assembly_q22).
可以通过调大预设值消除该提示and检查网格质量。
&&&ELEMENT 42
INSTANCE SOIL3-1 IS DISTORTING SO MUCH THAT IT TURNS
应改进单元质量
& 650 nodes are either missing intersection with
their respective master surface or outside the adjust zone.
&&&改改tie里的tolarance试试
&Dependent part instances cannot be edited or
assigned mesh attributes
模型树--assembly-打击part 右键--make
independent。也可以到模型树part步展开点mesh。
& The aspet ratio for nnn elements exceeds 100 to
&&&单元划分网格长宽比不合适。如果这些单元在不重要的区域(对结果肯定有些影响,
&&&影响大小取决于这三个单元的位置,在模型中的作用等),而且能计算,那就没问题了
&&&123&&elements
are distorted
&&&存在单元扭曲,如果这些单元在不重要的区域(对结果肯定有些影响,
&&&影响大小取决于这三个单元的位置,在模型中的作用等),而且能计算,那就没问题了(同15)
& ***WARNING: DEGREE OF FREEDOM 1 IS NOT ACTIVE ON
NODE 6 - THIS BOUNDARY
&& &CONDITION IS
&约束了单元没有得自由度对求解没有影响,可以查看下
&&&热分析时出现了这样的警告“THERE
IS ZERO HEAT FLUX EVERYWHERE
&& &There is
zero HEAT FLUX everywhere in the model based on the default
criterion. please check the value of the
&&&average&&HEAT
FLUX during the current iteration to verify that the HEAT FLUX is
small enough to be treated as zero.
&& & if not,
please use the solution controls to reset the criterion for zero
HEAT FLUX.
&& & 试试:
(1)是不是热源定义的问题,错误信息是说热源量几乎为零。
(2)定义热源的子程序调用命令流应该为*HEAT GENERATION,在材料模块中定义,子程序为HETVAL。
&& &The elements
in the element set WarnElemSurfaceIntersect-Step1 are
&&&surface
intersections. Refer to the status and message file for further
&&&检查一下你单元集合的定义以及面的定义,看是否出现了相交或重复定义的情况
&& &Boundary
conditions are specified on inactive dof of 36 nodes. The nodes
have been identified in node set
WarnNodeBCInactiveDof.
&& &Integration
and section point output variables will not be output for
deformable elements that are declared as rigid using
&& & *rigid body
这个仅是通知性质的(在interaction步设置为rigid
body,不输出应力应变),你在interaction步定义了刚体约束的话,都会出这个警告。&
&& & For a self
contact surface, the facets of the elements in element set
WarnElemFacetThickPt63d-Step1 are thicker than 0.6 times
&&&an edge or
diagonal length
&& &of the
facets. Use the MAXRATIO parameter on *SURFACE DEFINITION to allow
automatic rescaling of the contact
&&&thicknesses
where necessary for this surface.
&& &Refer to the
status file for further details.
&& &NO VALID
RADIATION OUTPUT REQUESTS HAVE BEEN GENERATED. THIS MAY BE DUE TO
EARLIER INPUT ERRORS OR SPECIFICATION
NONEXISTANT CAVITY OR SURFACE NAME
&&&检查一下你的output設定裡是不是有些set或surface沒有設定到
&123&&nodes may
have incorrect normal definitions. The nodes have been identified
in node set WarnNodeIncorrectNormal。
&&&先用看看WarnNodeIncorrectNormal在哪儿。这个不一定是致命的警告,有时候可以忽略。如果模型不收敛,可以检查下是否有过约束,
&&&在接触上存在边界条件or加载。
1&&上文已经说过,类似于
----------------------------------------------------------------------------------------------------------
&&&ERROR:Too
many attamps have been made
&& & Too many
attamps have been made....
SOLUTION APPEARS TO BE DIVERGING.
&& & CONVERGENCE
ISJUDGED UNLIKELY.
INCREMENT REQUIRED IS LESS THAN THE MINIMUM SPECIFIED
&& &Analysis
Input File Processor exited with an error.
-----------------------------------------------------------------------------------------------------------
&这样的信息几乎是无用信息(除了告诉你的模型分析失败以外,没有告诉你任何有用的东西)。宜再查找别的信息来考察。比如:
)&&Numerical&&sigularity&
& solver problem.
numerical&&sigularity
processing node105
instance&&表示:数值奇异:刚体位移(欠约束)&
&Zero pivot 表示:过约束
这样的信息(当然不仅仅是这些),才是比较有价值的。
&对于TIME INCREMENT REQUIRED IS LESS THAN THE
MINIMUM SPECIFIED
&& &Too many
attamps have been made
3& &对于“网格扭曲”的警告: excessively
distorted elements 前面有提到。
第一步:采用二楼底下的方法用display查看“ ErrElemExcessDistortion-Step1
”在模型的哪些部位,做到心中有数。
第二步:检查模型的网格质量: mesh步---verify----Analysis
Check选取模型。这种情况,一开始计算即出现“distorted
element”的信息。Besides:很多其他问题也会网格扭曲警告。比如,几何模型导入有误需要修补、单元类型选取错误、边界条件有误、材料属性错误、接触设置不合理、子程序错误等。
第三步:即使你的网格划分很好,如果变形过大,也会导致网格扭曲。然后修改网格划分,怒要出现红色,关键区域不要出现黄色。(当然最好是所有的网格都用structure划分,且都没有红色、黄色出现。网格质量就比较好。这种情况,警告信息往往是在计算到一定步骤之后才出现“distorted
element”。 这种情况建议采用ALE等方式,此不详述,搜索论坛“网格重划”“ALE”等技术。
1.强化、等向强化、随动强化定义
2.等向强化、随动强化理解
,qf*(ij)=0
f(ij)=f*(ij)C(q)=0
f(ij)=f*(ijij)=0
3.等向强化、随动强化应用范围
1、双线性随动强化(BKIN)
,所以有两个斜率,弹性斜率和塑性斜率,由于随动强化的Von
高达六条不同温度下的曲线。
3.表来定义非线性特性
2、双线性等向强化(BIS0)
Mises BKIN
MP,EX,1,200e9
MP,NUXY,1,0.25
MP,GXY,1,150e9
TBDATA,1,300e6,2000e6
3、多线性随动强化(MKIN)
随动强化效应,这个选项使用Von
Mises BKIN
MPTEMP1070
MPDATAEX330ES25ES
TBTEMPSTRA2N
TBDATA0.010.050.1
TBDATA300003700038000
TBDATA2250003100033000
4、多线性等向强化(MISO)
ANTYPE,STATIC
ET,1,solid45
MP,EX,1,2.1E11&&&&
MP,NUXY,1,0
TB,BKIN,1,1&&&&&&&&&
!TB,Biso,1,1&&&&&&&&&
TBDATA,1,250e6,8.6e9
,0.01,0.01,0.05
OUTRES,,1&&&&&&&&&&
WRITE SOLUTION ON RESULTS FILE FOR EVERY SUBSTEP
NSUBST,120&&&&&&&&&&&&&&&&&
! BEGIN WITH 120 SUBSTEPS
SFA,2,1,PRES,-&&&&
sfadele,2,,all
SFA,2,1,PRES,
sfadele,2,,all
ANSOL,2,84,EPEL,Z,EPELZ_2
ANSOL,3,84,EPPL,Z,EPPLZ_4
ANSOL,4,84,S,z,SZ_4
ADD,5,2,3, , , , ,1,1,1,
/AXLAB,X, DEFLECTION
/AXLAB,Y, Stress
模量”可以理解为是一种标准量或指标。材料的“模量”一般前面要加说明语,如弹性模量、压缩模量、剪切模量、截面模量等。这些都是与变形有关的一种指标。
杨氏模量(Young's Modulus):
杨氏模量就是弹性模量,这是材料力学里的一个概念。对于线弹性材料有公式σ(正应力)=Eε(正应变)成立,式中σ为正应力,ε为正应变,E为弹性模量,是与材料有关的常数,与材料本身的性质有关。杨(ThomasYoung)在材料力学方面,研究了剪形变,认为剪应力是一种弹性形变。
1807年,提出弹性模量的定义,为此后人称弹性模量为杨氏模量。(有点类似虎克定律^_^)
弹性模量(Elastic Modulus)E:
弹性模量E是指材料在弹性变形范围内(即在比例极限内),作用于材料上的纵向应力与纵向应变的比例常数。也常指材料所受应力如拉伸,压缩,弯曲,扭曲,剪切等)与材料产生的相应应变之比。
弹性模量是表征晶体中原子间结合力强弱的物理量,故是组织结构不敏感参数。在工程上,弹性模量则是材料刚度的度量,是物体变形难易程度的表征。
弹性模量E在比例极限内,应力与材料相应的应变之比。对于有些材料在弹性范围内应力-应变曲线不符合直线关系的,则可根据需要可以取切线弹性模量、割线弹性模量等人为定义的办法来代替它的弹性模量值。根据不同的受力情况,分别有相应的拉伸弹性模量modulus
of elasticity for tension (杨氏模量)、剪切弹性模量shear modulus of elasticity
(刚性模量)、体积弹性模量、压缩弹性模量等。
一个弹性常数,它等于应变(或应变分量)对应力(或应力分量)之比。对一个完善的弹性材料来说,它是弹性模量的倒数,即材料每单位应力的变形率。常见的实验测定的柔量有拉伸柔量、剪切柔量、蠕变柔量等。
剪切模量G(Shear Modulus):
剪切模量是指剪切应力与剪切应变之比。
剪切模数G=剪切弹性模量G=切变弹性模量G
切变弹性模量G,材料的基本物理特性参数之一,与杨氏(压缩、拉伸)弹性模量E、泊桑比ν并列为材料的三项基本物理特性参数,在材料力学、弹性力学中有广泛的应用。
其定义为:G=τ/γ,
其中G(Mpa)为切变弹性模量;
τ为剪切应力(Mpa);
γ为剪切应变(弧度)。
体积模量K(Bulk Modulus):
体积模量可描述均质各向同性固体的弹性,可表示为单位面积的力,表示不可压缩性。公式如下K=E/(3&(1-2&v)),其中E为弹性模量,v为泊松比。具体可参考大学里的任一本弹性力学书。
性质:物体在p0的压力下体积为V0;若压力增加(p0→p0+dP),则体积减小为(V0-dV)。则被称为该物体的体积模量(modulus
elasticity)。如在弹性范围内,则专称为体积弹性模量。体积模量是一个比较稳定的材料常数。因为在各向均压下材料的体积总是变小的,故K值永为正值,单位MPa。体积模量的倒数称为体积柔量。体积模量和拉伸模量、泊松比之间有关系:E=3K(1-2μ)。
压缩模量(Compression Modulus):
压缩模量指压应力与压缩应变之比。
储能模量E':
储能模量E'实质为杨氏模量,表述材料存储弹性变形能量的能力。储能模量表征的是材料变形后回弹的指标。
储能模量E'是指粘弹性材料在交变应力作用下一个周期内储存能量的能力,通常指弹性;
耗能模量E'':
耗能模量E''是模量中应力与变形异步的组元;表征材料耗散变形能量的能力, 体现了材料的粘性本质。
耗能模量E''指的是在一个变化周期内所消耗能量的能力。通常指粘性
切线模量(Tangent Modulus):
切线模量就是塑性阶段,屈服极限和强度极限之间的曲线斜率。是应力应变曲线上应力对应变的一阶导数。其大小与应力水平有关,并非一定值。切线模量一般用于增量有限元计算。切线模量和屈服应力的单位都是N/m2
截面模量:
截面模量是构件截面的一个力学特性。是表示构件截面抵抗某种变形能力的指标,如抗弯截面模量、抗扭截面模量等。它只与截面的形状及中和轴的位置有关,而与材料本身的性质无关。在有些书上,截面模量又称为截面系数或截面抵抗矩等。
强度是指某种材料抵抗破坏的能力,即材料抵抗变形(弹性塑性)和断列的能力(应力)。一般只是针对材料而言的。它的大小与材料本身的性质及受力形式有关。可分为:屈服强度、抗拉强度、抗压强度、抗弯强度、抗剪强度等。
如某种材料的抗拉强度、抗剪强度是指这种材料在单位面积上能承受的最大拉力、剪力,与材料的形状无关。
例如拉伸强度和拉伸模量的比较:他们的单位都是MPa或GPa。拉伸强度是指材料在拉伸过程中最大可以承受的应力,而拉伸模量是指材料在拉伸时的弹性。对于钢材,例如45号钢,拉伸模量在100MPa的量级,一般有200-500MPa,而拉伸模量在100GPa量级,一般是180-210Gpa。
刚度(即硬度)指某种构件或结构抵抗变形的能力,是衡量材料产生弹性变形难易程度的指标,主要指引起单位变形时所需要的应力。一般是针对构件或结构而言的。它的大小不仅与材料本身的性质有关,而且与构件或结构的截面和形状有关。
一般用弹性模量的大小E来表示.而E的大小一般仅与原子间作用力有关,与组织状态关系不大。通常钢和铸铁的弹性模量差别很小,即它们的刚性几乎一样,但它们的强度差别却很大。
最近瞅了一些关于时频分析工具箱的matlab函数使用方法,总结一下吧.
我使用的是2011a的matlab,貌似没有自带的时频分析工具箱,可以到网上下载,google一搜就能搜到,安装后就可以使用了(所谓安装就是把工具箱的目录包含到matlab工作目录中即可).
下面说一些时频工具箱函数的用法(由于我下的工具箱没有html版的帮助,所以只能使用help
funname的方法查看帮助信息).
Generate gaussian amplitude modulation.
Y=AMGAUSS(N,T0,T) generates a gaussian amplitude modulation
centered on a time T0, and with a spread proportional to T.
This modulation is scaled such that Y(T0)=1
and Y(T0+T/2) and Y(T0-T/2) are approximately equal to 0.5 .
产生高斯幅值调制信号,其以信号点数为N,中心为T0,传播0.5的比例到T,即Y(T0)=1,Y(T0+T/2) =
Y(T0-T/2) = 0.5。 T0默认为N/2,T默认为2*sqrt(N)。
如z=amgauss(160,90,40); plot(z);
Signal with constant frequency modulation.
[Y,IFLAW] = FMCONST(N,FNORM,T0) generates a frequency
modulation
with a constant frequency fnorm.
The phase of this modulation is such that y(t0)=1.
产生一个固定频率的频率调制信号。N为产生信号的点数,FNORM为标准化频率(默认为0.25),T0表示此时刻为正幅值(默认为round(N/2)),相当于规定了相位。IFLAW为设置的频率变化情况,这里就是一个常数。
如z=fmconst(128,0.05,50); plot(real(z));
Signal with linear frequency modulation.
[Y,IFLAW]=FMLIN(N,FNORMI,FNORMF,T0) generates a linear
modulation.
The phase of this modulation is such that Y(T0)=1.
产生一个线性调频信号。N为信号的点数,FNORMI为开始标准频率(默认为0.0),FNORMF为结束标准频率(默认为0.5),相位由T0规定,即Y(T0)=1。IFLAW为设置的频率变化情况,这里就是一个线性变化的直线。
如[z, f]=fmlin(128,0.05,0.3,50); plot(real(z));figure,
以上是比较常用的产生信号的函数。
下面说一些时频分析常用的分析方法,主要有短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD等,其他还有很多分析方法,都有相关函数,具体可以查看帮助。
短时傅里叶变换STFT
功能:计算时间序列的短时离散傅里叶变换,得到瞬时频率。
[tfr, t, f] = tfrstft(x) %
计算时间序列x的短时傅里叶变换,参数tfr为短时傅里叶变换系数,t为系数tfr对应的时刻,f为归一化频率向量
[tfr, t, f] = tfrstft(x, t) % 计算对应时刻t的短时傅里叶变换
[tfr, t, f] = tfrstft(x, t, n) % 计算n点对应时刻t的短时傅里叶变换
[tfr, t, f] = tfrstft(x, t, n, h) % 参数h为归一化频率平滑窗
[tfr, t, f] = tfrstft(x t, n, h, trace) % trace显示算法进程
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
h--频率滑窗,h归一化为单位能量(缺省值为hamming(n/4))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
sig=[fmconst(128,0.2);
fmconst(128,0.4)]; % 产生由两个常值调频信号(即正弦信号)的组合信号
tfr=tfrstft(sig);
subplot(211); imagesc(abs(tfr)); % 瞬时频率
subplot(212); imagesc(angle(tfr)); % 瞬时相位
Wigner-Ville时频分布图
功能:计算时间序列的Wigner-Ville时频分布图,得到瞬时频率
[tfr, t, f] = tfrwv(x)
[tfr, t, f] = tfrwv(x, t)
[tfr, t, f] = tfrwv(x, t, n)
[tfr, t, f] = tfrwv(x, t, n, trace)
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
sig=fmlin(128,0.1,0.4);
tfrwv(sig);
伪Wigner-Ville时频分布图
功能:计算时间序列的伪Wigner-Ville时频分布图,得到瞬时频率
[tfr, t, f] = tfrpwv(x)
[tfr, t, f] = tfrpwv(x, t)
[tfr, t, f] = tfrpwv(x, t, n)
[tfr, t, f] = tfrpwv(x, t, n, trace)
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
sig=fmlin(128,0.1,0.4);
tfrpwv(sig);
还可以设置成三维显示,如下
其他时频分析
计算时间序列的快速梅林变换
计算快速梅林逆变换。
生成时间序列的Bord-Jondan时频分布图,得到瞬时频率。
生成时间序列的Butterworth时频分布图,得到瞬时频率。
生成时间序列的Choi-Williams时频分布图,得到瞬时频率。
生成时间序列的广义矩形时频分布图,得到瞬时频率。
生成时间序列的Margenau-Hill时频分布图,得到瞬时频率。
生成时间序列的Margenau-HillP谱图时频分布图,得到瞬时频率。
计算时间序列的频谱图分布,得到瞬时频率。
计算时间序列的尺度图时频表示,即连续小波变换的幅值的平方。
计算时间序列的时频分布和它的重排形式,得到瞬时频率。
计算时间序列的瞬时频率密度。
从时间序列的重排时频分布中提取瘠。
计算给定时间序列的瞬时频率规律的理想时频表示。
以上给了很多函数,但常用的基本就是短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD和提取瘠ridges、得到理想时频表示tfrideal。
我们对信号进行频率分析时常用到的一个评价指标就是信号的功率谱密度函数,那么我们知道功率谱密度函数的自变量可以以频率f来表示,亦可以用圆频率w来表示,那么两种表示方法之间有什么量的关系呢?我们有时候容易弄混,本文笔者以最简单的限带白噪声信号为例,说明两者之间的转换关系。
设该白噪声的频率范围为0Hz-1Hz,即w的范围为0-2Pi,且设该白噪声信号为位移信号的功率谱密度函数。那么该信号可以用两个函数来描述,一个是以f为自变量的S1(f),一个是以w为自变量的S2(w)。取其中一个频率点w0和f0,它们代表同一个频率,即w0=2*pi*f0;大家注意:最容易犯的错误出现了:S1(f0)=S2(w0)。
你觉得同一频率处对应的功率谱密度值应该相同,其实这是错误的。要知道在功率谱密度的表达式中你计算出来的仅仅是函数值,并没有体现量纲。没错,同一频率处的功率谱密度相同,但在不同的量纲体系中,它们的值是不相等的。那我们就来看下两个函数的量纲分别是什么,看量纲的话我们从积分的角度看容易一些,我们知道功率谱密度图上曲线与横轴包围的面积相当于该信号在该带宽上的功率,即所测物理量的平方。对于S1(f),其函数值的量纲为m^2/Hz,而对于S2(w),其函数值的量纲为m^2/(rad/s),而我们知道1Hz=(2*pi)rad/s。所以在某频率点正确的表示为:S1(f0)m^2/Hz=S2(w0)m^2/(rad/s)=2*pi*S2(w0)m^2/Hz,即:S1(f0)=2*pi*S2(w0)。
结束了!最后做个小测验:某信号的功率谱密度函数为S(w)=w^2,那么其功率谱密度函数写成用f的表达式该怎么写?
答案:S2(f)=2*pi*S(2*pi*f)=8*pi^3*f^2.
希望大家都做对了!
大家都知道,在最普通的傅里叶变换公式中,自变量是圆频率与时间之间的转换,关系式如下:
而有时我们要得到的关系是以频率f为自变量的表达式,这时我们有时会常常就变换式前面的系数是否有1/2pi
而纠结,现在我就要说明这个问题。但是在讲这个问题之前,我先纠正一个大部分人都在犯的错误,不知道你在不在这些人中,那就是关于函数名的意义,即F(w)中的F代表的含义。函数名是函数的代号它代表的是由自变量到函数值之间的映射关系(规律),举例来说;f(x)=3x+1,中的f代表将自变量乘以3之后再加1的一种运算。
而f(x)=6x+1中的f则代表将自变量乘以6之后再加1的运算,由此我们可以知道,直接在那里写什么F(f)和F(w)之间的关系的文章都犯了这个错误。所以本文中进行了区分。
按上面的公式进行傅里叶变换之后得到一个自变量为w的函数,其函数值代表该信号在圆频率为w处的频率分量,显然圆频率w=w0处的函数值,应该与频率f=w0/(2*pi)处的函数值相等,因为在实际的物理意义上它们代表着同一种频率。由此可得下面的推导。
笔者在第一时间就更新的微软的windows8.1系统,但是使用了一段时间之后发现了各种问题。
第一,IE11极差的兼容性,新版的浏览器再界面和功能上看不出与IE10有什么区别,但在第三方扩展的兼容性方面却直线下降,现存软件插件几乎都与之不兼容,从而导致迅雷下载无法使用(期待新版迅雷赶紧开发),ie关闭时崩溃等现象,只得在高级选项中,关闭“启用第三方浏览器扩展”这一选项。
第二,新系统虽然集成了skydrive这一应用,并在电脑的硬盘上同时划出一部分空间用以脱机使用,这样的功能确实使得同步文件变得简单,无需通过上传操作,只需把文件移动或复制到skydrive文件夹就可以了。然而,笔者之前使用win8时常用的一个功能是skydrive中集成了onenote笔记本,在skydrive中可以直接调用onenote程序打开,而新版的系统是只在skydrive中保留了一个笔记本的快捷方式,并且根本无法调用程序打开。
第三,新系统的一个变化是在“这台电脑”下新增了若干文件夹,如视频,音乐,图片等,笔者到现在也是不知其所以然,不知道这东西是干什么用的,到是有点像xp系统中我的文档,我的音乐,我的图片那种感觉。
第四,新系统的微软中文输入法在使用很多程序时都会发生较为明显的卡顿情况,使得我不得不转向第三方输入法,而原来使用win8系统时,微软输入法已经能满足我的需求了。
首先声明本人使用的版本为Maple13,系统为win8专业版,最近装完Maple之后发现无法正常运行,总是出现“与内核之间连接失败”,后来在网上发现是系统不兼容,我可不想重装系统,看完网上各方英雄的回答,最终有效的一种是这样的:
把D:\Program Files (x86)\Maple 15\java 下把类似(可别手贱全给删了)mapleResources.jar
的选项全部删除,仅保留mapleResources.jar和mapleResourcesen.jar。
最后虽然变成英文的软件了,但可以正常使用,也就知足了!
用过Matlab的人都知道,Matlab的命令总是能一石三鸟,通过改变输入参数的注释项即可实现不同功能,如今xcorr命令的难点就在于其有四个注释项,这些注释项使得计算的结果各有不同,本文将详细介绍对应每个注释项Matlab是如何计算的,当然本文考虑输入的是一个简单一维序列x=[1,2,3],序列中数据对应的序号依次为1,2,3(请读者在阅读下文时,不要把序号和数据值弄混,这里只是个特例),其他情况读者可以轻松扩展得到。首先确定一下,该序列的均值为2,方差值为0.6667.这个应该不用说怎么算的吧。然后读者需要了解的是该命令意在计算序列中间隔不同距离的数据之间的关系。
第一:缺省注释项,[a,b]=xcorr(x),通过该命令计算的结果为:a=3 8 14 8
3;b=-2 -1 0 1 2.
下面介绍一下,该过程计算机是如何计算的,首先讲b的计算,设一维序列的长度为N,则序列中任意两个数据序号相减,最小值为1-N,最大值为N-1,且能取遍两者之间的所有整数,将这些数从小到大排列得到的就是b;然后讲a的计算,在缺省注释项的情况下,a的计算是这样的,a的每一项是对应b的每一项的
1、当b(1)=-2时,计算a(1)时只用到一组数据——(3,1)注意顺序,只有这两个数据的序号相减(后面数据的序号减去前面的)满足b=-2,因此a(1)的计算公式为:3*1=3
2、当b(2)=-1时,计算a(2)时用到两组数据——(2,1)和(3,2),这两组数据的序号相减(后面数据的序号减去前面的)满足b=-1,因此a(2)的计算公式为:2*1+3*2=8
3、当b(3)=0时,计算a(3)时用到三组数据——(1,1)、(2,2)、(3,3),这三组数据的序号相减(后面数据的序号减去前面的)满足b=0,因此a(3)的计算公式为:1*1+2*2+3*3=14
4、当b(4)=1时,计算a(4)时用到两组数据——(1,2)和(2,3),(读者请对比和情况2的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=1,因此a(4)的计算公式为:1*2+2*3=8
5、当b(5)=2时,计算a(4)时用到一组数据——(1,3),(读者请对比和情况1的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=2,因此a(4)的计算公式为:1*3=3
第二:注释项为‘unbiased’,[a,b]=xcorr(x,‘unbiased’),通过该命令计算的结果为:a=3
4 4.;b=-2 -1 0 1
2.下面介绍计算机如何计算该过程,b的计算在四种注释项的情况下是相同的,就不再讲述了。a的计算仍是和b的每一项相对应的。
1、当b(1)=-2时,计算a(1)时只用到一组数据(记N=1)——(3,1)注意顺序,只有这两个数据的序号相减(后面数据的序号减去前面的)满足b=-2,因此a(1)的计算公式为:(3*1)/1=3
2、当b(2)=-1时,计算a(2)时用到两组数据(记N=2)——(2,1)和(3,2),这两组数据的序号相减(后面数据的序号减去前面的)满足b=-1,因此a(2)的计算公式为:(2*1+3*2)/2=4
3、当b(3)=0时,计算a(3)时用到三组数据(记N=3)——(1,1)、(2,2)、(3,3),这三组数据的序号相减(后面数据的序号减去前面的)满足b=0,因此a(3)的计算公式为:(1*1+2*2+3*3)/3=4.6667
4、当b(4)=1时,计算a(4)时用到两组数据(记N=2)——(1,2)和(2,3),(读者请对比和情况2的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=1,因此a(4)的计算公式为:(1*2+2*3)/2=4
5、当b(5)=2时,计算a(4)时用到一组数据(记N=1)——(1,3),(读者请对比和情况1的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=2,因此a(4)的计算公式为:(1*3)/1=3
第三:注释项为‘biased’,[a,b]=xcorr(x,‘biased’),通过该命令计算的结果为:a=&
1.0000,b=-2 -1 0 1 2。下面介绍计算机如何计算该过程,注意到本次计算用到的序列x的长度为3,记为M=3。
1、当b(1)=-2时,计算a(1)时只用到一组数据——(3,1)注意顺序,只有这两个数据的序号相减(后面数据的序号减去前面的)满足b=-2,因此a(1)的计算公式为:(3*1)/M=1
2、当b(2)=-1时,计算a(2)时用到两组数据——(2,1)和(3,2),这两组数据的序号相减(后面数据的序号减去前面的)满足b=-1,因此a(2)的计算公式为:(2*1+3*2)/M=2.6667
3、当b(3)=0时,计算a(3)时用到三组数据——(1,1)、(2,2)、(3,3),这三组数据的序号相减(后面数据的序号减去前面的)满足b=0,因此a(3)的计算公式为:(1*1+2*2+3*3)/M=4.6667
4、当b(4)=1时,计算a(4)时用到两组数据——(1,2)和(2,3),(读者请对比和情况2的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=1,因此a(4)的计算公式为:(1*2+2*3)/M=2.6667
5、当b(5)=2时,计算a(4)时用到一组数据——(1,3),(读者请对比和情况1的区别),这两组数据的序号相减(后面数据的序号减去前面的)满足b=2,因此a(4)的计算公式为:(1*3)/M=1
第四:注释项为‘coeff’,[a,b]=xcorr(x,‘coeff’),通过该命令计算的结果为:a=
0.2143,b=-2 -1 0 1
2,下面介绍计算机如何计算该过程,这种情况实际是将第三种情况下得到的结果进行归一化,使得b=0时对应的值为1,a(1)=1/4.3;a(2)=2.7=0.5714,a(3)=4.7=1,a(4)=2.7=0.5714,a(5)=1/4.3
另:xcorr命令在工程上的应用通常是对时间上的采样数据序列x进行处理,当数据点采完之后交给Matlab处理时,Matlab是不知道你的采样时间间隔的,它仅仅根据上文所述的计算过程对输入的数据序列x进行计算,但我们可以自己定义时间间隔,例如dt=0.01,此时t=dt*b即代表相关性计算中的时间延迟,前半部分是超前,后半部分是滞后,若R=xcorr(x,‘unbiased’),则通过命令:plot(t,R)即可得到该时域信号的自相关函数曲线。
前几日,去上海车展看了一下,发现各种活动都要求扫二维码来参加,于是拿出手机,用Uc浏览器自带的二维码功能想扫一下,可是却蹦出这么一句话:“安卓相机出现问题,请重启设备”,重启就重启呗,谁怕谁?关键是重启不管用啊,于是上网查,发现网上各种说法都有,又是什么刷机造成的了,又是删除了邮件应用了,还有的说软件不兼容的,说换一个软件就好了。但这些说法,我都不以为然,于是照经验先看安全软件,打开之后发现安全软件在防录音、照相的模块中,把我的所有软件都设为了禁止,怪不得呢,于是手动设为允许,问题解决了!当然这也可能只是其中一种原因.

我要回帖

更多关于 负数运算法则 的文章

 

随机推荐