| 查看: 1032 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
MATLAB求解方程后没有报错，但无法绘制出曲线，求高手帮忙解决... 已有1人参与
|
||
|
最近利用MATLAB求解方程后没有报错，但无法绘制出曲线，求高手帮忙解决... 方程不是很复杂，具体方程如下： clear;clc; P=29.6; D=8245; R=0.021; LL=0.021; x0=0.010; pp=8960; a=0.002; b=0.002; erfa=30*pi/180; angle=atan(LL/R); n=10; h=a/n; N=x0/h; %步长总数 k=1.2; t=1/(2*k); Y=zeros(N,1); L=zeros(N,1); H=zeros(N,1); L(1)=R*tan(angle); H(1)=0; x(1)=h; y(1)=h*tan(erfa); for i=1:1:N-1; x(i+1)=(i+1)*h; %离散后xi横坐标 y(i+1)=x(i+1)*tan(erfa); %离散后yi纵坐标 psaiI=atan((LL+(i+1)*h)/R); z=fsolve(@(z)tan(psaiI)*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))-1/((2*z-1)^t)-sqrt((2*z-1)/z^2),1); Dm=D*sqrt(z^2/(2*z-1)); xx=fsolve(@(x)cos(x)/sin(psaiI-x)-sqrt(z^2/(2*z-1)),0); X=[x,xx]; L(i+1)=(L(i)*tan(X(i))*tan(psaiI)+tan(psaiI)*(R-H(i)))/(1+tan(X(i))*tan(psaiI)); H(i+1)=(R*tan(X(i))*tan(psaiI)-L(i)*tan(X(i))+H(i))/(1+tan(X(i))*tan(psaiI)); if H(i)>Y(i) Pm=z*P; else Pm=P; end if i==1 vv(1)=(1/3)*pi*(h*tan(erfa))^2*h; m(1)=pp*vv(1); v=Pm/m(1); elseif i>n vv(i)=(1/3)*pi*(i*h*tan(erfa))^2*i*h-vv(n)-(1/3)*pi*((i-n)*h*tan(erfa))^2*(i-n)*h; m(i)=pp*vv(i); v(i)=Pm/m(i); else vv(i)=(1/3)*pi*(i*h*tan(erfa))^2*i*h-vv(i-1); m(i)=pp*vv(i); v(i)=Pm/m(i); end end plot(x(i),v(i)) 想绘制横坐标为x, 纵坐标为v的曲线。请熟悉MATLAB的高手给点建议，谢谢... |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有16人回复
-
FileCode能看出啥？
已经有27人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有12人回复
-
静等基金结果
已经有19人回复
-
阿姨
已经有3人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
4楼2016-08-15 09:56:39
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-13 20:47:14
|
谢谢您的回复。现在我将程序拆分出来分别计算求得x，z值; H值；Pm值；和微元质量m。最后绘制plot(x,v)。发现还是不行！ 发现：这样可以得到离散的x(i)值，y(i)值，X(i)值。但是离散的z(i)值得不到！所以最后无法得到Pm=z*p值。请问给出意见如何修改？谢谢... clear;clc; P=29.6; D=8245; R=0.021; LL=0.021; x0=0.010; pp=8960; a=0.002; b=0.002; erfa=30*pi/180; angle=atan(LL/R); n=10; h=a/n; N=x0/h; %步长总数 k=1.2; t=1/(2*k); %求X值和z值 for i=1:N x(i)=i*h; y(i)=x(i)*tan(erfa); psaiI=atan((LL+i*h)/R); z=fsolve(@(z)tan(psaiI)*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))-1/((2*z-1)^t)-sqrt((2*z-1)/z^2),1); xx=fsolve(@(x)cos(x)/sin(psaiI-x)-sqrt(z^2/(2*z-1)),0); X(i)=xx; end %求H值 L=zeros(N,1); H=zeros(N,1); L(1)=R*tan(angle); H(1)=0; for i=1:N-1 L(i+1)=(L(i)*tan(X(i))*((LL+i*h)/R)+((LL+i*h)/R)*(R-H(i)))/(1+tan(X(i))*((LL+i*h)/R)); H(i+1)=(R*tan(X(i))*((LL+i*h)/R)-L(i)*tan(X(i))+H(i))/(1+tan(X(i))*((LL+i*h)/R)); end %比较上面计算出的H和y值的大小后，求Pm值 if H(i)>y(i) Pm=z*P; else Pm=P; end %求各个微元的质量 for i=1:N if i==1 vv(1)=(1/3)*pi*(h*tan(erfa))^2*h; %第一个微元体积 m(1)=pp*vv(1); v(1)=Pm/m(1); elseif i>n vv(i)=(1/3)*pi*((i)*h*tan(erfa))^2*(i)*h-vv(n)-(1/3)*pi*((i-n)*h*tan(erfa))^2*(i-n)*h; m(i)=pp*vv(i); v(i)=Pm/m(i); else vv(i)=(1/3)*pi*((i)*h*tan(erfa))^2*(i)*h-vv(i); m(i)=pp*vv(i); v(i)=Pm/m(i); end end plot(x,v) |
3楼2016-08-15 06:44:45
5楼2016-08-16 03:36:42
