| 查看: 583 | 回复: 5 | ||
[求助]
虫友们，看一下这个程序，运行一下 已有1人参与
|
| 这个程序我就想将输出的四个系数值，全部编辑到coefficient的矩阵中，然后一起输出来，不想一次输出一个来，这样的话，后面紧接着画图，感觉不能直接画在一幅图上，还有就是为什么这个程序最后会输出个ans=。。。。，求助谢谢 |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
| 不好意思，附件没传上 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：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
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
谢谢，图画出来了，只是有一个问题，为什么最后出现个ans，你是怎么做到的？ 6楼2015-11-17 20:12:08