kcjf.net
当前位置:首页 >> MATLAB求教 >>

MATLAB求教

clc;clear f1=@(x)[x(1)-0.6*sin(x(1))-0.3*cos(x(2)),x(2)-0.6*cos(x(1))-0.3*sin(x(2))]; f2=@(x)[x(1)^2-x(2)^2+x(3)-9,3*x(1)+5*x(2)+6*x(3),x(1)-3*x(2)-6*x(3)-1]; x=fsolve(f1,[0,0]) x2=fsolve(f2,[-1,1,-1]) 结果:Equation solved.fsolve completed because the

你看看我编的这个合适吗? r=2;for i=0:1:5; x1=r*cos(2*pi*i/5); y1=r*sin(2*pi*i/5); x2=0; y2=0; x3=r*cos(2*pi*(i+1)/5); y3=r*sin(2*pi*(i+1)/5); x=[x1 x2 x3]; y=[y1 y2 y3]; mypatch=patch(x,y,[i*0.1/5 i*0.5/5 i*0.4/5]); %(x, y, color) hold on axis equalend

你只要把数据排列成一个断面矩阵就行:以深度和断面宽度作为横纵坐标,每个格点为数据,然后contour或contourf或imagesc就可以.你仔细想想其实是很简单的.

原式可变为:dy/dx=1/3/(x+1)y(x=0)=0;由于表达式较简单,可以直接用dsolvedsolve('Dy=1/3/(x+1)','y(0)=0','x')

好像没有这个函数rref(A) 是初等行变换化行简化梯矩阵,但它用了交换两行看来你只能自己写函数了

做个测试,希望有所帮助. 代码 % By lyqmathclc; clear all; close all;syms x tft = limit(t*(1+1/x)^(2*t*x), x, inf);diff(ft, t) 结果 ans = exp(2*t) + 2*t*exp(2*t) >>

<p>x0是初值,那迭代后变化的应该是x,那就是画x100次迭代的值的图.</p> <p>把function [x,val,k]=grad(fun,gfun,x0)几个地方改下</p> <p>k=0; epsilon=1e-5;</p> <p></p> <p>x=zeros(2,maxk);</p> <p></p> <p> k=k+1;</p> <p>x(:,k)=x0;</p>

L是个m*n矩阵L(k,j)是一个元素的值你这里L(k,j)作为area向量下标时能够保证不越界吗?

>>x=[520,280,156.6,78,39.62,3.1,0];>>x1=[0,3.1,39.62,78,156.6,280,520];>>y=[0,-30,-36,-35,-28.44,-9.4,0];>>y1=[0,9.4,28.44,35,36,30,0];>>x2=520:-10.4:0;>>y2=spline(x,y,x2);>>x3=0:10.4:520;>>y3=spline(x1,y1,x3);>>plot(x2,y2,'-',x3,y3,'-') 其实要加入更多的点的话,你只需要把自变量x的范围按要求变化下就好了.但是要注意的就是样条函数的自变量x要求是递增或是递减,所以分成两段来解.

网站首页 | 网站地图
All rights reserved Powered by www.kcjf.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com