| 查看: 737 | 回复: 5 | |||
阿毛虫银虫 (正式写手)
|
[求助]
虫友们，看一下这个程序，运行一下 已有1人参与
|
| 这个程序我就想将输出的四个系数值，全部编辑到coefficient的矩阵中，然后一起输出来，不想一次输出一个来，这样的话，后面紧接着画图，感觉不能直接画在一幅图上，还有就是为什么这个程序最后会输出个ans=。。。。，求助谢谢 |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有8人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
关于Filecode分析方法
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 547.7
- 散金: 81
- 帖子: 757
- 在线: 117.2小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
| 不好意思，附件没传上 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : clamped_coefficient.m
2015-11-17 08:43:57, 1.84 K
2楼2015-11-17 08:44:05
【答案】应助回帖
感谢参与，应助指数 +1
| 改了一下，你可以看看 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : clamped_coefficient.m
2015-11-17 16:03:10, 1.95 K
3楼2015-11-17 16:03:17
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 547.7
- 散金: 81
- 帖子: 757
- 在线: 117.2小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
4楼2015-11-17 16:09:06
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
阿毛虫: 金币+8, ★★★很有帮助, 谢谢 2015-11-17 20:12:20
阿毛虫: 金币+8, ★★★很有帮助, 谢谢 2015-11-17 20:12:20
|
function [orthogonality, coefficient_all]=clamped_coefficient(c) format long clear all G=8.0e10; K=1.2; E=2.06e11; P=7.9e3; l=4; D=0.6; H=0.2; k=1/K; I=D*(H^3)/12; A=D*H; figure; coefficient_all = []; r=sqrt(I/A); s=sqrt(E*I/(k*A*G)); syms x ; syms c; frequency=[405.35308 1093.84137 2087.15447 3341.29751]; for ii=1:length(frequency) b=sqrt(P*A/(E*I))*frequency(ii); alfa=b*sqrt(-(r^2+s^2)+sqrt((r^2-s^2)^2+4./(b^2)))/sqrt(2); beta=b*sqrt(r^2+s^2+sqrt((r^2-s^2)^2+4./(b^2)))/sqrt(2); y=cosh(alfa*x)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sinh(alfa*x)..., -cos(beta*x)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sin(beta*x); y_int=int(y^2,0,l); fai=(alfa^2+b^2*s^2)/alfa*sinh(alfa*x)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2))*(alfa^2+b^2*s^2)/alfa*cosh(alfa*x)..., -(b^2*s^2-beta^2)/beta*sin(beta*x)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2))*(b^2*s^2-beta^2)/beta*cos(beta*x); fai_int=int(fai^2,0,l); orthogonality=c^2*y_int+c^2*r^2*fai_int-1/(P*A); coefficient_result=solve(orthogonality,c) ; coefficient_results=vpa(coefficient_result,6); coefficient = abs(coefficient_results(1)); coefficient_all = [coefficient_all; coefficient]; x_shuzhi=0:0.01:4; y_shuzhi=cosh(alfa.*x_shuzhi)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sinh(alfa.*x_shuzhi)..., -cos(beta.*x_shuzhi)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sin(beta.*x_shuzhi); hold on; h(ii) = plot(x_shuzhi,y_shuzhi); end |
5楼2015-11-17 16:15:45
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 547.7
- 散金: 81
- 帖子: 757
- 在线: 117.2小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
谢谢，图画出来了，只是有一个问题，为什么最后出现个ans，你是怎么做到的？ 6楼2015-11-17 20:12:08
