k5vin
2011-09-26, 10:12
本人新手,最近有一方程不得解,往高手相助:
clear all, close all
syms x;
y0=0.12/0.8; y1=1-y0;
hg=2/20;
c1=1.8624*hg^2-3.95*hg^3+16.375*hg^4-37.226*hg^5+76.81*hg^6 ...
-126.9*hg^7+172*hg^8-143.97*hg^9+66.56*hg^10;
h=0.02; L=0.8;
c0=5.346*h*c1;
T=[0 1 0 1 0 0 0 0
1 0 1 0 0 0 0 0
0 0 0 0 0 -1 0 1
0 0 0 0 -1 0 1 0
sin(x*y0) cos(x*y0) sinh(x*y0) cosh(x*y0) -sin(x*y1) -cos(x*y1) -sinh(x*y1) -cosh(x*y1)
cos(x*y0)-c0*x*sin(x*y0) -sin(x*y0)-c0*x*cos(x*y0) cosh(x*y0)+c0*x*sinh(x*y0) sinh(x*y0)+c0*x*cosh(x*y0) cos(x*y1) -sin(x*y1) cosh(x*y1) sinh(x*y1)
-sin(x*y0) -cos(x*y0) sinh(x*y0) cosh(x*y0) sin(x*y1) cos(x*y1) -sinh(x*y1) -cosh(x*y1)
-cos(x*y0) sin(x*y0) cosh(x*y0) sinh(x*y0) -cos(x*y1) sin(x*y1) cosh(x*y1) sinh(x*y1)];
f=det(T);
现在必须用到fzero函数求解f。求解
clear all, close all
syms x;
y0=0.12/0.8; y1=1-y0;
hg=2/20;
c1=1.8624*hg^2-3.95*hg^3+16.375*hg^4-37.226*hg^5+76.81*hg^6 ...
-126.9*hg^7+172*hg^8-143.97*hg^9+66.56*hg^10;
h=0.02; L=0.8;
c0=5.346*h*c1;
T=[0 1 0 1 0 0 0 0
1 0 1 0 0 0 0 0
0 0 0 0 0 -1 0 1
0 0 0 0 -1 0 1 0
sin(x*y0) cos(x*y0) sinh(x*y0) cosh(x*y0) -sin(x*y1) -cos(x*y1) -sinh(x*y1) -cosh(x*y1)
cos(x*y0)-c0*x*sin(x*y0) -sin(x*y0)-c0*x*cos(x*y0) cosh(x*y0)+c0*x*sinh(x*y0) sinh(x*y0)+c0*x*cosh(x*y0) cos(x*y1) -sin(x*y1) cosh(x*y1) sinh(x*y1)
-sin(x*y0) -cos(x*y0) sinh(x*y0) cosh(x*y0) sin(x*y1) cos(x*y1) -sinh(x*y1) -cosh(x*y1)
-cos(x*y0) sin(x*y0) cosh(x*y0) sinh(x*y0) -cos(x*y1) sin(x*y1) cosh(x*y1) sinh(x*y1)];
f=det(T);
现在必须用到fzero函数求解f。求解